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Abstract 


Red  Sprites  and  Blue  Jets  are  two  different  types  of  recently  discovered  optical  flashes  ob¬ 
served  above  large  thunderstorm  systems.  Sprites  are  luminous  glow's  occurring  at  altitudes 
typically  ranging  from  ~50  to  90  km.  In  video  they  exhibit  a  red  color  at  their  top  which 
gradually  changes  to  blue  at  lower  altitudes.  Sprites  may  occur  singly  or  in  clusters  of  two 
or  more.  The  lateral  extent  of  “unit”  sprites  is  typically  5-10  km  and  they  endure  for  several 
milliseconds.  Jets  are  upward  moving  (~  1 00  km/s)  highly  collimated  beams  of  luminosity, 
emanating  from  the  tops  of  thunderclouds,  extending  up  to  ~50  km  altitude  and  exhibiting 
a  primarily  blue  color. 

We  propose  that  sprites  result  from  large  electric  field  transients  capable  of  causing 
electron  heating,  breakdown  ionization  and  excitation  of  optical  emissions  at  mesospheric 
altitudes  following  the  removal  of  thundercloud  charge  by  a  cloud-to-ground  discharge.  De¬ 
pending  on  the  history  of  charge  accumulation  and  removal,  and  the  distribution  of  ambient 
atmospheric  conductivity,  the  breakdown  region  may  have  the  shape  of  vertically  oriented 
ionization  column(s).  Results  of  a  two-dimensional  and  selfconsistent  quasi-electrostatic 
(QE)  model  indicate  that  most  of  the  observed  features  of  sprites  can  be  explained  in  terms 
of  the  formation  and  self-driven  propagation  of  streamer  type  channels  of  breakdown  ion¬ 
ization.  Comparison  of  the  optical  emission  intensities  of  the  1st  and  2nd  positive  bands  of 
N2,  Meinel  and  1st  negative  bands  of  NJ,  and  the  1st  negative  band  of  0}  demonstrates  that 
the  1st  positive  band  of  N2  is  the  dominant  optical  emission  in  the  altitude  range  ~50-90 
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km,  which  accounts  for  the  observed  red  color  of  sprites.  Optical  emissions  of  the  1st  and 
2nd  positive  bands  of  N2  occur  in  carrot-like  vertical  structures  with  typical  transverse  di¬ 
mension  ~5-10  km  which  can  span  an  altitude  range  from  ~80  km  to  well  below  ~50  km. 
The  appearance  of  optical  emissions  associated  with  sprites  can  be  delayed  in  time  (~l-20 
ms)  with  respect  to  the  causative  cloud  to  ground  discharge.  Theoretical  model  results  are 
found  to  be  in  good  agreement  with  recent  video,  photometric  and  spectral  measurements 
of  sprites. 

We  propose  that  blue  jets  are  streamer  type  processes  occurring  on  an  atmospheric  spatial 
scale  during  the  charge  accumulation  stage  in  a  thundercloud  prior  to  lightning  discharges. 
The  quasi-electrostatic  model  produces  results  in  general  agreement  with  most  features  of 
jets  as  observed  in  video.  Results  indicate  that  the  velocity  of  upward  propagation  of  jets 
is  ~  100  km/s,  and  that  the  general  shape  of  jets  is  in  the  form  of  upward  expanding  beams 
of  luminosity  with  cone  angles  <30°.  The  model  results  agree  with  video  observations 
showing  that  the  brightness  of  jets  drops  ~2  orders  of  magnitude  between  the  lower  part 
and  the  tip  of  the  jet,  and  the  fact  that  the  luminosity  fades  away  simultaneously  everywhere 
along  the  cone  of  the  jet.  The  blue  color  of  jets  as  observed  in  video  is  naturally  explained 
as  emission  from  the  2nd  positive  band  of  N2  excited  by  electron  impact.  Jets  are  not 
necessarily  associated  with  lightning  discharges  and  may  appear  only  in  relatively  rare 
cases  of  large  (~300-400  coulombs)  thundercloud  charge  accumulation  at  high  (~20  km) 
altitudes.  Jets  occur  during  the  accumulation  phase  of  the  thundercloud  charge,  so  that  it  is 
not  necessary  for  large  amounts  of  charge  to  be  removed  in  order  to  produce  them. 
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In  this  introductory  chapter  we  briefly  survey  the  information  available  concerning  the 
recently  discovered  phenomena  of  optical  flashes  above  thunderstorms,  spanning  the  altitude 
range  from  cloud  tops  and  up  to  the  lower  edge  of  the  ionosphere  (~90  km),  and  referred 
to  as  “sprites”  and  “blue  jets”. 

Sprites  are  luminous  glows  occurring  at  altitudes  typically  ranging  from  ~50  to  90  km.  In 
video,  they  exhibit  a  red  color  at  their  top  which  gradually  changes  to  blue  at  lower  altitudes. 
Sprites  may  occur  singly  or  in  clusters  of  two  or  more  [Sent man  et  al.,  1995].  The  lateral 
extent  of  “unit"  sprites  is  typically  5-10  km  and  they  endure  for  several  milliseconds.  Blue 
jets  or  jets  in  short  are  upward  moving  (~  1 00  km/s)  highly  collimated  beams  of  luminosity, 
emanating  from  the  tops  of  thunderclouds,  extending  up  to  ~50  km  altitude  and  exhibiting 
a  primarily  blue  color  [Wescott  et  al.,  1995a], 

The  intent  of  this  introductory  survey  is  to  give  the  reader  a  general  understanding  of  the 
observed  properties  of  the  phenomenon  of  sprites  and  jets  which  constitute  an  important 
component  of  the  complex  of  physical  processes  which  are  not  yet  fully  understood  and 
which  manifest  upward  electrodynamic  coupling  of  energy  originating  in  tropospheric  thun¬ 
derstorms  to  the  mesosphere  and  lower  ionosphere.  The  scope  of  this  dissertation  work  and 
the  specific  contributions  to  knowledge  are  discussed  at  the  end  of  this  chapter  in  Section  1 .3. 
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1.1  UPWARD  ELECTRODYNAMIC  COUPLING  OF  THUNDERSTORMS 

A  thunderstorm  is  the  final  product  of  the  growth  of  a  series  of  clouds  which  are  formed 
through  the  upward  transfer  of  energy  from  surface  heating  effects  leading  to  buoyant  ver¬ 
tical  air  motions  or  through  the  lifting  effects  associated  with  the  motion  of  atmospheric 
frontal  weather  zones.  The  thunderstorm  electromagnetic  activity  is  spectacular  and  dan¬ 
gerous,  in  spite  of  the  fact  that  it  represents  an  energy  expenditure  of  less  than  1%  of  the 
thermal,  gravitational,  and  kinetic  energy  associated  with  condensation  of  moisture  and  the 
development  of  drafts  and  precipitation  [Bra ham,  1952]. 

The  importance  of  tropospheric  thunderstorms  in  the  coupling  of  electromagnetic  energy 
to  the  upper  atmospheric  regions  on  a  global  scale  is  underscored  by  the  fact  that  the 
thunderstorms  and  associated  lightning  activity  are  the  generators  which  maintain  the  more 
than  5  x  1 05  Coulombs  of  the  earths  charge  against  the  conduction  loss  due  to  the  total  global 
current  of  ~  1 800  A  flowing  through  the  conducting  atmosphere  in  fair  weather  regions  [e.g., 
Stergis ,  1957].  There  are  more  than  2000  thunderstorms  active  at  any  given  time  over  the 
Earth's  surface  and  on  average  lightning  strikes  the  Earth  ~100  times  per  second  [Volland, 
1984].  On  the  order  of  1000  coulombs  of  bound  space  charge  are  separated  and  maintained 
during  the  late  growth  and  active  periods  of  a  thunderstorm  [e.g.,  Marshall  et  al.,  1996] 
with  cloud-to-ground  lightning  discharges  involving  the  transfer  to  the  ground  of  up  to  300 
coulombs  of  charge  in  several  ms  [ Brook  et  al.,  1982],  leading  to  large  quasi-electrostatic 
(QE)  fields  which  exist  in  the  mesosphere  and  lower  ionosphere  regions  over  millisecond 
time  scales.  These  fields,  combined  with  intense  electromagnetic  pulses  (EMP)  of  ~20 
Gigawatts  peak  power  [liman,  1987,  p.  1 15]  generated  by  lightning  currents,  are  thought 
to  intensely  heat  the  ambient  electrons  [Inan  et  al.,  1996b,  and  references  therein]  and 
accelerate  energetic  runaway  electrons  (REL)  [Roussel-Dupre  et  al.,  1994;  Bell  etal.,  1995]. 
producing  ionization  and  optical  emissions. 
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1.1.1  Experimental  Evidence 

The  experimental  evidence  indicating  strong  upward  electrodynamic  coupling  of  thunder¬ 
storms  to  the  mesosphere  and  lower  ionosphere  includes  early/fast  perturbations  of  subiono- 
spherically  propagating  VLF  signals  associated  with  lightning  discharges  in  underlying 
thunderstorms  [Armstrong,  1983;  Inan  et  al.,  1988;  1993;  1995;  1996a;  Dowden  et  al., 
1994].  and  optical  emissions  in  clear  air  above  thunderstorms  associated  with  sprites  [e.g., 
Sentman  et  al.,  1995;  Lyons,  1994,  1995,  1996;  Boeck  et  al.,  1995;  Rairden  and  Mende, 
1 995;  Winckler  et  al.,  1996]  and  blue  jets  [  Wescott  et  al.,  1995a]  and  airglow  enhancements 
[ Boeck  et  al.,  1992].  The  latter  are  believed  to  be  closely  associated  with  recently  discov¬ 
ered  rapid  (<  1  ms)  optical  emissions  at  80-95  km  altitudes  with  lateral  extents  up  to  300 
km  preceding  sprites  and  referred  to  as  “elves”  [Fukunishi  et  al.,  1996;  Lyons,  1996;  Inan 
et  al.,  1996b],  Sprites  are  clearly  associated  both  temporally  and  spatially  with  intense 
lightning  discharges  in  underlying  thunderstorms  [Boccippio  et  al.,  1995]  as  well  as  with 
early/fast  VLF  events  [Inan  et  al.,  1995]  providing  strong  evidence  of  significant  changes 
in  the  mesospheric  electrical  conductivity.  Sprites  are  almost  uniquely  associated  with  pos¬ 
itive  cloud-to-ground  lightning  discharges  [Boccippio  et  al.,  1995;  Winckler  et  al.,  1996; 
Lyons,  1996].  The  first  observations  of  gamma  ray  bursts  of  terrestrial  origin  [Fishman 
et  al.,  1994],  associated  with  positive  cloud-to-ground  lightning  discharges  [Inan  et  al., 
1996c],  and  intense  wideband  VHF  bursts  [e.g.,  Holden  et  al.,  1995]  constitute  additional 
new  examples  of  upward  coupling  of  energy  originating  in  thunderstorms. 

1.1.2  Theoretical  Mechanisms 

The  new  observations  listed  above  have  led  to  the  discovery  of  new  interaction  mechanisms 
which  result  in  troposphere-mesosphere/lower  ionosphere  coupling.  These  mechanisms  are 
based  on  heating  of  the  ambient  electrons  by  electromagnetic  pulses  generated  by  lightning 
discharges  [Inan  et  al.,  1991,  1996b;  Taranenko  et  al.,  1993a,b;  Milikh  et  al.,  1995;  Rowland 
etal.,  1995]  or  by  large  quasi-electrostatic  thundercloud  fields  [Paskoet  al.,  1995;  1996a,b.c; 
Boccippio  et  al.,  1995;  Winckler  et  al.,  1996],  and  on  runaway  electron  processes  [Bell 
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et  al.,  1995;  Winckler  et  al.,  1996;  Roussel-Dupre  and  Gurevich,  1996;  Taranenko  and 
Roussel-Dupre ,  1996].  It  has  been  proposed  by  Inan  et  al.  [1996d]  that  quasi-electrostatic 
thundercloud  fields  are  capable  to  maintain  the  ionospheric  electrons  at  a  persistently  heated 
level  well  above  their  ambient  thermal  energy.  Changes  in  the  thundercloud  charge  (e.g., 
in  lightning  discharges)  lead  to  heating/cooling  above/below  this  quiescent  level,  and  are 
registered  as  early/fast  VLF  events  [Inan  et  al.,  1996d].  The  simultaneous  observations  of 
early/fast  VLF  events  and  sprites  [Inan  et  al.,  1995],  indicate  that  these  two  effects  may  be  a 
manifestation  of  the  same  physical  process  consistent  with  model  predictions  of  both  optical 
emissions,  heating  and  ionization  changes  associated  with  quasi-electrostatic  thundercloud 
fields  [Pasko  et  al.,  1995,  1996a,b;  Inan  et  al.,  1996d].  Some  of  these  physical  processes 
are  illustrated  in  Figure  1.1. 

It  should  be  noted  that  in  existing  scientific  literature  the  possibility  of  breakdown  ioniza¬ 
tion  of  the  upper  atmosphere  by  thundercloud  fields  was  first  mentioned  in  1925  by  C.T.R. 
Wilson.  He  recognized  that  the  relation  between  the  thundercloud  electric  field  which  de¬ 
creases  with  altitude  ;  as  ~  and  the  critical  breakdown  field  which  falls  more  rapidly 
(being  proportional  to  the  exponentially  decreasing  atmospheric  density)  leads  to  the  re¬ 
sult  that  “there  will  be  a  height  above  which  the  electric  force  due  to  the  cloud  exceeds 
the  sparking  limit”  [Wilson,  1925].  Three  decades  later  Wilson  speculated  that  a  discharge 
between  the  top  of  a  cloud  and  the  ionosphere  might  often  be  the  normal  accompaniment 
of  a  lightning  discharge  to  ground  [Wilson,  1956]. 

In  this  dissertation  we  investigate  the  heating,  ionization,  attachment  and  optical  emis¬ 
sion  effects  associated  with  intense  quasi-electrostatic  thundercloud  fields  existing  above 
large  thunderstorms  before  (producing  blue  jets)  and  after  (generating  sprites)  lightning 
discharges.  In  the  following  section  we  provide  a  brief  description  of  the  phenomenology 
of  optical  flashes  above  thunderstorms  on  the  basis  of  selected  results  of  video,  photometric 
and  spectroscopic  measurements. 
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Phenomena 


ANC,  ATC,  Ahv 


RCVR 


Mechanisms 


ANe,  ATC,  Ahv 


Fig.  1.1.  Phenomena  and  mechanisms.  Illustration  of  different  phenomena  (top  panel) 
and  theoretical  mechanisms  (bottom  panel )  of  lightning-ionosphere  interactions  operating 
at  different  altitudes  and  producing  optical  emissions  (Ah')  observed  as  sprites,  blue  jets  and 
elves,  as  well  as  heating  (AT)  and  ionization  changes  (A Ne)  detected  as  very  low  frequency 
(VLF)  signal  changes. 
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1.2  OPTICAL  FLASHES  ABOVE  THUNDERSTORMS 

The  most  spectacular  evidence  of  electrodynamic  coupling  between  tropospheric  lightning 
and  the  overlying  mesosphere/lower  ionosphere  are  luminous  glows  above  thunderstorms 
occurring  at  altitudes  ranging  from  the  cloud-tops  to  ~90  km.  Visual  accounts  of  this 
kind  of  glow  in  the  clear  air  above  thunderstorms  have  appeared  in  the  literature  since  the 
19th  century,  the  most  vivid  accounts  being  those  by  air  transport  pilots  [Vaughan  and 
Vonnegut,  1989].  Although  most  of  these  reports  probably  described  the  known  phenomena 
of  lightning  channels  which  terminate  in  the  clear  air  rather  than  in  the  cloud  [e.g.,  Malan, 
1963,p.  6 \Uman,  1969,  p.  2],  there  were  exceptional  cases  that  differed  from  the  description 
of  a  typical  air  discharge.  In  retrospect,  these  unusual  diffuse  optical  events  are  likely  to 
have  been  associated  with  recently  documented  optical  flashes  above  thunderstorms  called 
sprites  and  jets. 

1.2.1  First  Recorded  Image 

The  first  recorded  image  of  unusual  optical  flashes  occupying  large  volumes  of  space  above 
thunderstorms  was  obtained  serendipitously  on  5  July  1989  during  a  test  of  a  low-light-level 
TV  camera  at  the  O’Brien  Observatory  of  the  University  of  Minnesota  near  Minneapolis. 
This  observation  is  described  by  Franz  et  al.  [1990],  and  a  reproduction  of  the  image  from 
this  paper  is  shown  in  Figure  1 .2. 

1.2.2  Observations  from  Space  Shuttle 

Following  the  original  observation  by  Franz  et  al.  [1990],  Vaughan  et  al.  [1992]  and 
Boeck  et  al.  [1995]  reported  video  observations  from  the  space  shuttle  of  tens  of  similar 
optical  events.  Most  of  the  video  examples  captured  from  space  have  the  appearance  of  thin 
luminous  filaments  of  glowing  gases  stretching  into  the  stratosphere  above  thunderstorm 
[Boeck  et  al.,  1995].  Some  of  the  video  images  obtained  from  the  space  shuttle  [Vaughan 
et  al.,  1992;  Boeck  et  al.,  1995]  are  reproduced  in  Figure  1 .3. 
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Fig.  1.2.  The  first  TV  image  of  an  optical  flash  above  thunderstorms  \Franz  et  al.. 
1990].  Low-lighi-level  TV  image  of  a  luminous  electrical  discharge  over  a  thunderstorm 
~250  km  from  the  observing  site.  The  image,  which  persisted  less  than  one  video  frame 
(<30  ms),  resembled  two  plume  discharges  simultaneously  extending  into  the  clear  sky 
above  thunderstorm.  The  image  is  taken  from  http://wwwghcc.msfc.nasa.gov/  skeets.html 
[Courtesy  of  Dr.  O.  H  Vaughan,  Jr.]. 


Typical  heights  of  optical  flashes  were  estimated  to  be  60  to  75  km  above  the  earth  with 
the  uncertainty  arising  from  the  rather  coarse  resolution  of  the  digitized  image  and  from 
the  lack  of  information  on  the  cloud  top  height  or  the  location  at  which  the  stratospheric 
flash  intersects  the  cloud  top  [Boeck  et  al.,  1 995].  The  lateral  width  of  observed  luminosities 
varied  considerably  between  events,  with  some  examples  showing  thin  or  several  thin  vertical 
filaments  while  others  appearing  as  broad  columns  of  some  kilometers  across.  In  addition, 
many  events  exhibited  a  bulge  of  luminosity  at  the  top  of  the  flash  with  dimensions  of 
order  of  kilometers  [Boeck  et  al..  1995],  The  full  vertical  extent  of  the  flashes  were  only 
observable  within  the  time  resolution  of  one  video  field  (17  ms)  [Boeck  et  al.,  1995].  It 
appears  from  the  space  shuttle  observations  that  the  conditions  for  the  stratospheric  flashes 
may  occur  over  most  regions  of  the  globe  (in  temperate  and  tropical  areas,  over  the  oceans, 
and  over  the  land)  and  that  the  luminosity  is  definitely  associated  with  a  discharge  in  the 
cloud  below  and  occurs  soon  after  the  onset  of  the  discharge  [Boeck  et  al.,  1995]. 
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Fig.  1.3.  Observations  of  optical  flashes  above  thunderstorms  from  the  space  shuttle 
[Boeck  et  ai,  1995].  A  single  captured  image  from  each  stratospheric  flash'  (a)  A  distinct 
bulge  or  blob  of  illumination  at  the  upper  end  of  the  stratospheric  flash.  This  perception  of 
an  upper  bulge  is  further  enhanced  in  this  and  other  cases  by  the  break  in  the  illumination 
found  just  below  the  upper  bulge.  The  length  of  the  vertical  column  was  estimated  to  be 
approximately  31  km  f Vaughan  et  al .,  1992];  (b)  A  single-filament  stratospheric  flash,  with 
a  single  break  in  illumination,  and  a  bulge  of  illumination  at  the  top  observed  on  October 
6,  1990,  at  23:37:06  UTC.  The  flash  extended  above  the  horizon;  (c)  A  stratospheric  flash 
(01 :29:48  UTC)  illustrates  very  clearly  a  fan-shaped  structure  at  the  top  of  the  filament;  (d) 
A  double-filament  stratospheric  flash  observed  two  minutes  following  previous  event  (c). 
The  image  is  taken  from  http://wwwghcc.msfc.nasa.gov/  skeets.html  [Courtesy  of  Dr.  O. 
H.  Vaughan,  Jr.]. 
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1.2.3  Observations  During  Summer  1993 

During  the  summers  of  1993-1994  there  occurred  a  marked  increase  in  number  of  reports 
of  optical  flashes  above  thunderstorms. 

Using  a  low  light  level  all  sky  black  and  white  television  system,  Sentnum  and  Wescott 
[  1 993]  imaged  nineteen  examples  of  optieal  flashes  during  a  single  flight  of  NASA's  DC-8 
airborne  laboratory  over  thunderstorms  in  Iowa,  Nebraska  and  Kansas  in  July,  1993.  They 
estimated  the  terminal  heights  of  the  events  to  be  60  km,  with  error  bars  (due  to  unknown 
range)  extending  to  100  km.  The  duration  of  the  flashes  was  estimated  to  be  16  ms  or  less, 
and  the  brightness  as  calculated  from  comparison  with  stellar  brightness  to  be  25-50  kR. 
roughly  that  of  moderately  bright  aurorae.  The  emission  volume  of  flashes  was  estimated 
to  be  >  1000  km',  with  a  horizontal  extent  of  ~  10-50  km  and  their  occurrence  rate  was 
such  that  an  optical  flash  was  observed  in  association  with  every  200-300th  cloud-to-ground 
stroke. 

In  a  ground-based  program  in  the  high  plains  area  of  Fort  Collins,  Colorado,  Lyons  [  1 994] 
obtained  black  and  white  images  of  high  altitude  flashes  above  storm  systems  in  Nebraska 
and  Kansas  in  July  and  August,  1993.  An  image-intensified,  low-light  video  camera  system¬ 
atically  monitored  the  stratosphere  above  distant  (100-800  km  range)  mesoscale  convectne 
systems  over  the  high  plains  of  the  central  United  States  for  21  nights  between  6  July  and 
27  August  1993.  Complex,  luminous  structures  were  observed  above  large  thunderstorm 
clusters  on  eleven  nights,  with  one  storm  system  (7  July  1993)  yielding  248  events  in  410 
minutes.  The  duration  of  events  ranged  from  33  to  283  ms,  with  an  average  of  98  ms. 
The  luminous  structures,  generally  not  visible  to  the  naked,  dark-adapted  eye.  exhibited  on 
video  a  wide  variety  of  brightness  levels  and  shapes  including  streaks,  aurora-like  curtains, 
smudges,  fountains  and  jets.  The  structures  were  often  more  than  10  km  wide  and  their 
upper  portions  extended  to  above  50  km  altitude.  Some  of  the  shapes  appeared  to  be  “frozen 
in  space"  as  if  “illuminated"  during  the  few  video  fields  in  which  they  were  visible.  The 
video  images,  with  a  time  resolution  of  17  ms  per  video  field,  did  not  provide  any  sense  of 
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Fig.  1.4.  An  optical  flash  observed  from  the  ground  by  L\m?[  1994].  The  evcnl  observed 
at  0607  UTC  on  July  7,  1 993.  The  distant  storm  cloud  top  heights  are  estimated  to  be  between 
10  and  15  km  altitude.  If  so,  the  uppermost  extent  of  the  structures  would  be  from  50  to 
65  km  altitudes  [Lyons.  1994],  The  image  is  taken  from  http://wwwghcc.msfc.nasa.gov/ 
skects.html  [Courtesy  of  Dr.  O.  H.  Vaughan.  Jr.]. 


propagation  either  upward  or  dow  nward.  Figure  1 .4  illustrates  a  typical  image  of  an  optical 
flash  observed  on  July  7.  1993  [Lyons,  1994], 

With  heights  reaching  towards  100  km  and  diameters  of  20  km  or  more,  assuming  cir¬ 
cular  symmetry  suggests  that  the  larger  luminous  structures  can  attain  emission  volumes 
of  thousands  of  cubic  kilometers.  While  the  mesoscale  convective  system  as  a  whole  is 
electrically  active,  mesospheric  optical  flashes  tend  to  occur  in  regions  characterized  by 
low  cloud-to-ground  flash  rates  where  positive  cloud-to-ground  flashes  are  more  common 
[Lvwj.v,  1994]. 

Winckler  [1995]  observed  more  than  150  optical  flashes  above  a  thunderstorm  complex 
moving  south-east  across  the  state  of  Iowa  during  the  night  of  9-10  August  1993.  Black 
and  white  images  were  obtained  through  clear  air  with  an  intensified  CCD  TV  cameras  at 
the  O’Brien  Observatory  of  the  University  of  Minnesota  located  about  60  km  north  east 
of  Minneapolis  and  250-500  km  from  the  storm  center.  The  optical  flashes  consisted  of 
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bright  vertical  striations  extending  from  50-80  km  altitude,  often  covering  tens  of  kilometers 
laterally,  with  tendrils  of  decreasing  intensity  visible  in  the  brighter  events  down  to  the  cloud 
tops  below  20  km  altitude.  All  of  the  more  intense  events  were  coincident  with  a  VLF  sferic 
(i.e.,  powerful  impulsive  radio  signal  radiated  by  intense  return  stroke  currents  involved  in  a 
lightning  discharge)  measured  in  the  300Hz-12  kHz  range,  but  smaller  events  often  did  not 
occur  in  association  with  a  detectable  sferic.  The  intense  part  of  the  events  usually  occurred 
within  one  TV  field  ( 1 6.7  ms).  The  brightness  for  the  more  intense  events  was  estimated  to 
be  50-100  kR. 

Comparison  of  results  obtained  with  those  reported  by  Vaughan  et  al.  [  1992],  Sentman 
and  Wescott  1 1 993],  and  Lyons  [  1 994]  led  Winckler  to  clearly  conclude  that  all  these  authors 
were  indeed  observing  the  same  phenomena  [Winckler,  1995]. 

1.2.4  Terminology  for  Optical  Flashes  Above  Thunderstorms 

The  optical  flashes  above  thunderstorms  have  from  the  beginning  been  variously  referred  to 
in  the  literature  as  “upward”,  “cloud-to-stratosphere”,  or  “cloud-to-ionosphere”  lightning 
discharges.  Use  of  these  terms  suggests  of  unwarranted  parallels  to  normal  tropospheric 
lightning,  and  also  implies  an  upward  direction  of  propagation  which  has  not  yet  been 
determined  [see  discussion  in  Lyons,  1994  and  Sentman  et  al.,  1995].  The  term  “sprite" 
has  been  suggested  by  D.  Sentman  [e.g.,  Sentman  et  al..  1995]  as  non-judgmental  as  to 
the  physics  of  the  phenomenon  and  the  specific  production  mechanisms.  All  of  the  above 
described  optical  flashes,  including  those  observed  from  the  space  shuttle,  might  qualify  for 
inclusion  in  the  broad  category  of  sprites. 
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Fig.  1.5.  An  optical  event  recorded  on  color  video  from  an  aircraft  [Sentman  et  al., 

1995].  A  sprite  event  (one  of  the  largest)  observed  on  4  July  1994  al  0400:20  UT  [Sentman 
et  al..  1995],  The  original  color  image  [Sentman  et  al..  1995]  is  reproduced  here  in  black 
and  white.  The  image  is  taken  from  hllp://wwwghce.msfc. nasa.gov/  skcels.html  [Courtesy 
of  Dr.  O.  H.  Vaughan,  Jr.]. 

1.2.5  Observations  of  Sprites  During  Summer  1994 

The  first  color  imagery  of  optical  flashes  above  thunderstorms  and  their  unambiguously 
triangulated  physical  dimensions  and  heights  were  obtained  using  observations  from  two 
jet  aircraft  [Sentman  et  al.,  1995 \Wescott  etal.,  1995a;  Sentman  and  Wescott,  1994,  1995]. 
Low  light  level  television  images,  in  both  color  and  black  and  white,  show  that  there  are  at 
least  two  distinctly  different  types  of  optical  emissions  spanning  part  or  all  of  the  distance 
between  the  anvil  tops  of  the  thunderclouds  and  the  ionosphere.  The  first  of  these  emissions, 
called  sprites,  are  luminous  structures  of  brief  (typically  <  16  ms  although  some  events 
endure  longer)  duration  with  a  main  body  in  red  color  that  typically  span  the  altitude  range 
of  50-90  km.  and  lateral  dimensions  of  5-30  km.  Faint  bluish  tendrils  often  extended 
downward  from  the  main  body  of  sprites,  occasionally  appearing  to  reach  the  cloud  tops 
near  20  km.  Figure  1 .5  shows  one  of  the  brightest  events  recorded  during  the  campaign. 

The  brightest  red  region  of  a  unit  sprite,  its  “head”,  lies  between  the  characteristic  altitudes 
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of  66  and  74  km.  Above  the  head  there  is  often  a  dark  band  (“hairline”)  above  which  a 
faint,  wispy  red  glow  (“hair”)  is  observed  to  splay  upward  and  outward  toward  terminal 
altitude.  Beneath  the  head  there  sometimes  occurs  a  dark  band  (“collar”)  at  66±4  km 
below  which  faint  tendrils  may  extend  downward  to  altitudes  of  40-50  km,  changing  from 
red  near  the  collar  to  blue  at  their  lowest  extremities.  Sprites  may  occur  singly,  but  more 
typically  occur  in  clusters  of  two  or  more.  Clusters  may  be  tightly  packed  together  into 
large  structures  40  km  or  more  across,  or  loosely  spread  out  into  distended  structures  of 
spatially  separated  sprites.  The  onset  of  luminescence  occurs  simultaneously  (within  video 
resolution  of  ~16  ms)  across  the  cluster  as  a  whole,  coincident  with  the  occurrence  of 
cloud-to-ground  lightning  below.  All  elements  of  the  cluster  typically  decay  in  unison  in 
100-  !60  ms,  but  most  of  this  decay  time  may  be  attributed  to  the  image  lag  of  the  video 
cameras  [Sentman  el  al..  1095].  Sprites  seem  preferentially  to  occur  on  the  backside  of 
frontal  storm  systems,  and  occur  in  regions  of  positive  lightning  ground  strokes. 

Observations  of  sprites  from  the  Yucca  Ridge  Field  Station  29  km  northeast  of  Fort  Collins. 
Colorado,  on  the  night  of  July  11-12,  1 994  were  performed  by  Winckler  et  al.  [  1 996]  using 
both  w ide-fteld  and  telescopic  image-intensified  CCD  TV  cameras,  a  telescopic  photometer 
system  and  a  1  to  50  kHz  broad  band  VLF  sferic  receiver.  Telescopic  TV  images  of  bright 
sprites  had  a  fan  shaped  upper  plume  with  not  well  resolved  fine  features,  as  well  as  upward 
forked  and  vertically  striated  forms  adjacent  to  these  plumes  and  bright  points  of  luminosity 
around  the  plume  shaped  regions.  Many  sprites  consisted  entirely  of  groups  of  vertically 
aligned  striations  which  sometimes  appeared  to  diverge  from  a  common  point  of  origin  at 
cloud  tops.  All  observed  sprites  were  preceded  by  a  cloud-to-ground  stroke  with  a  coincident 
sfcric  and  sky  flash.  All  cloud-to-ground  strokes  associated  with  sprites  were  positive,  and 
most  had  100  kA  or  more  inferred  peak  current.  Using  the  photometer  data,  the  duration 
of  sky  flashes  induced  by  the  cloud-to-ground  strokes  was  determined  to  be  ~3  ms  and  the 
additional  sprite  total  light  curve  also  endured  for  ~3  ms.  An  image  of  a  large  bright  sprite 
event  recorded  at  05:58:45.371  UT  is  shown  in  Figure  1.6.  This  sprite  displays  vertical 
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Fig.  1.6.  A  bright  sprite  event  recorded  on  the  night  of  July  11-12, 1994  by  Winckleret 
al.  [1996]. 

striations  on  each  side  of  a  central  core.  This  sprite  brightened  in  the  next  field  (not  shown) 
and  greatly  faded  one  video  field  later. 

1.2.6  Observations  of  Sprites  During  Summer  1995 

Video  images  of  sprites  were  obtained  by  Rairden  and  Mende  [1995]  with  a  CCD  camera 
during  a  July  1995  ground-based  campaign  near  Fort  Collins,  Colorado.  These  unfiltered 
intensified  images  reveal  detailed  spatial  structure  within  the  sprite  envelope.  The  temporal 
resolution  of  the  standard  interlaced  video  imagery  is  limited  by  the  60  fields  per  second 
acquisition  rate  (16  ms).  The  specific  CCD  used  here,  however,  is  subject  to  bright  events 
leaking  into  the  readout  registers,  allowing  time  resolution  on  the  order  of  the  linescan  rate 
(63  /zs).  Typical  sprite  onset  is  found  to  follow  the  associated  cloud-to-ground  lightning 
by  1.5  to  4  ms.  The  onsets  of  the  individual  sprites  within  a  cluster  are  generally,  but  not 
always,  simultaneous  to  within  1  ms.  Sprites  tend  to  have  a  bright  localized  core,  less  than 
2  km  in  horizontal  dimension,  which  rises  to  peak  intensity  within  0.3  ms  and  maintains 


INTRODUCTION 


15 


Fig.  1.7.  A  collage  of  three  sprite  images  obtained  b>  Rairden  and  Mende  [1995].  The 
image  is  reproduced  from  [Rairden  and  Mende .  1995] 

this  level  for  5  to  10  ms  before  fading  over  an  additional  10  ms.  Three  examples  of  captured 
sprite  images  illustrating  fine  spatial  structure  are  shown  in  Figure  1.7. 

1.2.7  Photometric  Properties  of  Sprites 

It  was  clear  from  the  very  first  video  observations  that  essentially  all  the  temporal  and 
spatial  development  of  the  sprites  occurred  within  one  ~16  ms  video  field  [Winckler  et  al., 
1996],  Thus,  photometric  measurements  are  necessary  to  extract  useful  information  on  the 
temporal  dynamics  of  sprites  on  time  scales  less  than  16  ms. 

Winckler  el  al.  [1996]  used  a  photometric  system  with  high  time  resolution  with  the 
analog  photomultiplier  outputs  sampled  at  100  kHz  rate.  The  sprites  and  associated  cloud 
and  sky  flashes  in  each  event  were  identified  on  the  associated  video  images.  Although  the 
video  provided  only  16.7  ms  timing  accuracy,  for  most  events  millisecond  timing  accuracy 
could  be  obtained  from  the  VLF  sferic-photometer  correlated  responses  [Winckler  et  al., 
1996].  An  example  is  shown  in  Figure  1 .8,  covering  30  ms  of  data  which  included  a  bright 
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sprite.  The  cloud-to-ground  stroke  sferic  is  coincident  with  the  first  rise  of  the  photometer 
within  1  ms  starting  at  05:33:28.788  UT  on  July  12,  1994.  The  photometer  trace  is  double 
peaked.  Winckler  et  al.  [1996]  interpret  the  smaller  initial  3  ms  duration  peak  as  Rayleigh 
scattered  light  in  the  sky  due  to  the  intense  light  of  the  cloud-to-ground  discharge  (termed 
sky  flash  in  Figure  1 .8)  and  the  second  larger  peak  3  ms  later  as  direct  light  of  the  sprite  in 
the  photometer  field  of  view. 
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Fig.  1.8.  Photometer  and  sferic  traces  for  a  sprite  event  [Winckler  et  al.,  1996].  The 
initial  photometer  response  within  1  ms  of  the  sferic  is  followed  by  a  small  maximum  due 
to  the  scattered  light  from  the  cloud-to-ground  stroke,  and  a  large  maximum  due  to  the 
luminosity  of  the  sprite.  The  figure  is  reproduced  from  [Winckler  et  al.,  1996]. 


The  puzzling  feature  that  the  total  duration  of  video  images  of  sprites  was  often  longer  than 
that  indicated  by  the  photometric  data  can  be  understood  in  terms  of  the  different  spectral 
response  of  the  two  instruments,  coupled  with  the  spectral  time  history  of  the  sprite  event. 
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The  photometer  used  by  Winckler  et  al.  [1996]  was  blue  sensitive  with  a  response  peaked 
at  450  nm.  In  contrast,  the  video  photocathodes  were  red  sensitive,  with  broad  maximum 
extending  to  850  nm.  Results  suggest  that  the  sprite  began  with  an  intense,  wideband  optical 
flash  covering  both  the  blue  and  red  portions  of  the  spectrum.  The  video  output,  covering 
primarily  the  red  portion  of  the  spectrum,  persisted  for  up  to  10  video  fields  in  contrast  to 
the  photometer  output,  primarily  covering  the  blue  region,  which  decayed  on  a  time  scale 
of  one  video  field  [Winckler  et  al.,  1996]. 

Observations  of  sprites  and  related  optical  phenomena  were  performed  by  Fukunishi 
et  al.  [1996]  using  a  multichannel  high-speed  (time  resolution  of  15  /rs)  photometers  at 
Yucca  Ridge  Field  Station,  Colorado  in  June-August  1995  .  The  photometers  were  used 
to  simultaneously  observe  two  different  elevation  angles  (0.2°  pointing  accuracy)  with  a 
narrow  field  of  view  of  1 .0°  in  the  vertical  direction  and  9.5°  in  the  horizontal  direction.  A 
typical  example  of  transient  luminous  events  observed  is  shown  in  Figure  1 .9.  The  location 
of  a  positive  discharge  (08:57:53.306  UT,  +326  kA)  coincident  in  time  with  this  event  was 
identified  by  National  Lightning  Detection  Network  (NLDN)  to  be  at  a  distance  of  231  km 
South-East  of  Yucca  Ridge.  The  figure  shows  the  luminosity  variations  measured  by  the 
photometer  over  a  40  ms  interval,  with  one  photometer  channel  (without  filter)  set  at  15° 
elevation,  while  the  remaining  two  channels,  one  without  filter  and  the  other  with  a  red  filter, 
set  at  9°  elevation.  The  luminosity  variations  consist  of  two  parts,  an  initial  short  duration 
burst  coincident  in  time  with  the  recorded  onset  of  a  VLF  sferic  and  a  longer  duration 
event  starting  ~7  ms  later.  The  first  burst  actually  consisted  of  two  parts,  with  an  initial 
portion  interpreted  to  be  due  to  Rayleigh  scattered  light  from  the  precursor  cloud-to-ground 
lightning  flash.  The  second  portion  was  observed  only  in  the  high  elevation  photometer  and 
was  interpreted  to  be  a  large  scale  lower  ionospheric  diffuse  flash  lying  between  ~75  to  90 
km  altitudes.  These  brief  flashes  are  called  “elves”  [Fukunishi  et  al.,  1996]  and  are  believed 
to  be  produced  by  the  heating  of  ambient  ionospheric  electrons  by  intense  electromagnetic 
pulses  (EMPs)  radiated  by  powerful  lightning  discharges  [Inan  et  al.,  1996b],  The  longer 
duration  peak  which  started  ~7  ms  later  corresponds  to  a  sprite  event  [Fukunishi  et  al., 
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1996].  The  time  delay  between  the  onset  of  the  VLF  sferic  and  the  onset  of  the  following 
sprite  is  found  in  the  case  shown  to  be  7  ms;  however,  in  different  cases  observed,  this  time 
delay  ranged  from  several  to  several  tens  of  milliseconds  [Fukunishi  et  al.,  1996]. 

In  summary,  the  most  recent  photometric  measurements  demonstrate  that  sprites  are 
typically  delayed  in  time  ranging  from  ~3  ms  [Winckler  et  al.,  1996]  to  several  tens  of  ms 
[Fukunishi  et  al.,  1996]  with  respect  to  the  onset  of  the  causative  cloud-to-ground  stroke. 
The  duration  of  sprite  events  as  observed  in  photometer  data  usually  decays  on  a  time  scale 
of  one  video  field,  although  this  result  may  be  specific  to  the  spectral  range  covered  by  the 
photometers  used. 

1.2.8  Spectral  Properties  of  Sprites 

Spectra  of  sprites  in  the  stratosphere/mesosphere  above  electrically  active  cumulonimbus 
clouds  were  acquired  on  July  16,  1995  from  an  observation  site  near  Ft.  Collins,  Colorado 
by  Mende  et  al.  [1995].  The  spectra,  resolved  from  ~450-800  nm,  included  four  distinct 
features  in  the  600-760  nm  region  which  have  been  identified  as  the  N2  first  positive  system 
with  A=2,3,  and  4  from  the  v=2,4,5,6  vibrational  levels  of  the  B^FIg  state.  The  spectra  were 
lacking  in  other  features  such  as  the  Nj  Meinel  or  the  N£  first  negative  system  which  led  the 
authors  [Mende  et  al.,  1995]  to  a  conclusion  that  the  electron  energy  causing  the  excitation 
is  quite  low.  However,  with  the  available  signal-to-noise  ratio  of  the  measurement,  in  order 
to  be  detectable  these  other  features  would  have  had  to  be  at  levels  no  less  than  a  factor  of 
ten  of  that  of  the  N2  first  positive  system.  The  result  obtained  by  Mende  et  al.  [1995]  is 
shown  in  Figure  1.10. 

Hampton  et  al.  [  1 996]  used  a  video  slit  spectrograph  to  obtain  optical  spectra  of  sprites. 
Twenty  five  events  were  observed  from  the  Mt.  Evans  Observatory  over  a  thunderstorm 
on  the  border  of  Nebraska  and  Colorado  between  0700  and  0900  UT  on  the  night  of  22 
June  1995.  For  10  of  these  events  optical  spectra  were  measured  in  the  wavelength  range 
from  540  to  840  nm.  After  correcting  for  the  spectrograph  response  function,  digitized 
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Fig.  1.9.  Lightning-induced  transient  luminous  event  [Fukunishi  et  ai,  1996].  Anexam 
pie  of  lightning-induced  transient  luminous  events  observed  at  Yucca  Ridge  Field  Station, 
Colorado:  (a)  luminosity  variations  observed  by  a  four-  channel  photometer  during  the  40- 
ms  interval  between  08h  57m  53.300  sec  and  53.340  sec  UT  on  June  23,  1995;  (b)  enlarged 
display  of  photometer  signals  for  the  initial  3-ms  interval.  The  figure  in  reproduced  from 
[Fukunishi  et  cil.  1996]. 
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Fig.  1.10.  Spectral  scan  of  a  sprite  [Mende  et  al.,  1995).  The  thick  trace  represents  the 
spectra  of  the  sprite  while  the  thin  trace  shows  the  spectra  of  a  calibration  neon  light.  The 
figure  is  reproduced  from  [ Mende  etai,  1995). 

spectrograph  video  images  were  used  to  measure  the  wavelength  of,  and  ratios  between,  the 
various  emissions.  All  emissions  were  found  to  be  those  of  the  first  positive  bands  of  N2. 
There  was  no  evidence  of  the  presence  Meinel  bands  of  N£,  indicating  (in  the  opinion  of 
the  authors)  that  the  mechanism  responsible  for  sprites  produces  little  or  no  ionization  at  70 
km  altitude.  However,  it  is  important  to  note  that,  in  order  to  be  detectable,  the  intensity  of 
the  Meinel  band  would  have  to  be  no  smaller  than  a  factor  of  ~10  below  that  of  the  N2  first 
positive  band.  As  we  show  later  (Sections  3.5  and  4.7.3),  the  predicted  levels  of  Meinel 
band  intensities  in  sprite  producing  processes  which  lead  to  intense  ionization  columns  are 
~10*  lower  than  that  of  N2  first  positive. 

In  summary,  the  red  color  of  sprites  was  experimentally  identified  by  two  independent 
groups  of  authors  as  being  primarily  due  to  excitation  of  the  first  positive  band  of  molecular 
nitrogen  [Mende  et  al.,  1995;  Hampton  et  al.,  1996].  The  measurements  of  both  Hampton  et 


INTRODUCTION 


21 


Fig.  1.11.  Reduced  spectrum  of  a  sprite  [Hampton  et  at.,  1996],  Reduced  spectrum  from 
a  sprite  event  observed  at  07:04:39  UT.  The  wavelengths  for  the  band  heads  for  different 
transitions  of  the  first  positive  band  of  N2  arc  marked  at  top.  The  vertical  lines  show  the 
auroral  intensities  of  these  transitions.  The  figure  is  reproduced  from  [ Hampton  et  ai, 
I996|. 


al.  [  1996]  and  Mende  et  al.  [1995]  are  fully  consistent  with  the  theoretical  results  put  forth 
in  this  thesis  and  are  not  inconsistent  with  intense  ionization  occurring  within  the  luminous 
regions  of  sprites.  Quantitative  interpretation  of  these  results  is  given  in  Sections  3.5  and 
4.7.3  of  this  dissertation. 

1.2.9  Observations  of  Blue  Jets  During  Summer  1994 

During  the  Sprites  1994  dual  aircraft  observational  campaign  conducted  by  the  University 
of  Alaska  a  completely  new  and  distinctly  different  (from  sprites)  optical  phenomena  called 
“blue  jets"  was  captured  by  black  and  white  and  color  video  cameras  [  Wescott  et  al.,  1995a; 
Sentman  and  Wescott,  1994],  Blue  jets  are  narrowly  collimated  beams  of  blue  light  that 
propagate  upwards  from  the  tops  of  thunderclouds.  The  jets  appear  to  propagate  upward 
at  speeds  of  about  100  km/s  and  reach  terminal  altitudes  of  40-50  km.  Fifty  six  examples 
were  recorded  during  a  22  minute  interval  during  a  storm  over  Arkansas.  Two  examples  of 
the  images  recorded  are  illustrated  in  Figure  1.12. 
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Fig.  1.12.  Two  images  of  Blue  Jets  obtained  from  an  aircraft  during  summer  1994 
Sprites  campaign  [Wescott  et  ai.  1995a].  Images  are  provided  through  the  courtesy  of 
E.  Wescott,  University  of  Alaska.  Images  are  taken  from  http://wwwghcc.msfc.nasa.gov/ 
skeets.html  [Courtesy  of  Dr.  O.  H.  Vaughan,  Jr.]. 
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Blue  jets  are  narrowly  collimated  (full  cone  angle  ~  15±7.5°)  with  an  apparent  fan  out 
near  their  tops  (40-50  km).  The  estimated  brightness  ranges  up  to  nearly  1 03  kR  [E.  Wescott, 
personal  communication].  Their  apparent  source  duration  is  ~200  ms  at  the  base  of  the 
jet.  The  overall  brightness  at  the  bottom  is  ~  100  times  that  at  the  top  of  the  jet  and  the 
luminosity  decays  simultaneously  along  the  jet  beginning  at  about  200  to  300  ms.  A  faint 
hemispherical  “shock”  is  observed  beyond  the  terminus  of  some  jets  traveling  at  the  same 
speed  as  the  earlier  rising  portion  of  the  jet.  The  typical  lateral  dimensions  of  the  jets  are  few 
kilometers  at  the  base  and  <  10  km  at  the  top.  Unlike  sprites,  blue  jets  may  be  characterized 
by  a  lack  of  strong  electromagnetic  emission  at  VLF  and  are  not  observed  to  be  associated 
with  simultaneous  cloud-to-ground  flashes  [Wescott  et  al.,  1995b;  1996], 
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1.3  SCIENTIFIC  CONTRIBUTIONS 

The  purpose  of  this  Ph.D.  dissertation  research  was  to  theoretically  investigate  new  forms  of 
electrodynamic  coupling  of  energy  originating  in  lightning  discharges  to  the  mesosphere  and 
the  lower  ionosphere.  The  scientific  contributions  resulting  from  this  effort  are  described 
in  Chapters  3, 4  and  5.  A  summary  of  the  specific  contributions  to  knowledge  is  presented 
below: 

•  Identification  of  a  new  type  of  electrodynamic  coupling  between  tropospheric  thun¬ 
derstorms  and  the  mesosphere  and  the  lower  ionosphere  based  on  heating,  ionization  and 
optical  emission  effects  produced  by  quasi-electrostatic  thundercloud  fields. 

•  Development  of  a  new  two  dimensional  and  selfconsistent  model  of  the  nonlinear 
upward  coupling  of  quasi-electrostatic  thundercloud  fields. 

•  Quantitative  evaluation  of  electron  heating,  ionization  of  the  ambient  neutral  atmo¬ 
sphere,  and  the  excitation  of  optical  emissions  due  to  intense  quasi-electrostatic  fields  fol¬ 
lowing  lightning  discharges.  Comparison  of  results  with  experimental  data  on  mesospheric 
optical  flashes  called  sprites. 

•  Investigation  of  the  formation  and  upward  propagation  of  streamer  type  ionization  chan¬ 
nels  produced  by  pre-discharge  quasi-electrostatic  fields  immediately  above  the  thunder¬ 
cloud.  Identification  of  this  process  as  the  underlying  physical  mechanism  for  the  observed 
upward  moving  highly  collimated  beams  of  blue  color  emanating  from  the  thundercloud 
tops  and  known  as  blue  jets. 

Most  of  the  results  presented  in  this  dissertation  were  previously  published  in  [Pasko  et 


al,  1995,  1996a,b], 
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Quasi-Electrostatic  Heating  Model 


In  this  chapter,  we  describe  our  model  formulation  of  the  interaction  of  quasi-electrostatic 
thundercloud  fields  with  the  mesosphere  and  lower  ionosphere. 


2.1  ELECTRIC  FIELD  AND  CHARGE  DENSITY 

The  quasi-electrostatic  field  is  established  in  the  mesosphere  and  lower  ionosphere  due  to 
the  accumulation  of  thundercloud  charge  and  its  evolution  in  time  as  a  portion  of  the  charge 
is  lowered  to  the  ground  by  a  positive  cloud-to-ground  (+CG)  or  negative  cloud-to-ground 
(-CG)  lightning  stroke,  or  neutralized  inside  the  cloud  due  to  an  intracloud  (IC)  discharge 
(Figure  2. 1 ).  In  addition,  (see  Section  3.8  for  more  details)  we  consider  the  specific  cases 
of  the  removal  or  deposition  of  monopole  charges  (see  cases  marked  Ml  and  M2  in  the  top 
panel  of  Figure  2. 1 ). 

The  quasi-electrostatic  fields  treated  in  our  model  are  the  slowly  varying  and  long  enduring 
electric  component  of  the  total  electromagnetic  field  which  is  generated  by  the  removal  of 
charge  from  the  cloud.  Specifically,  we  leave  out  the  short  duration  electromagnetic  pulses 
(EMP)  generated  primarily  by  the  return  stroke  currents  and  which  produce  bright  optical 
emissions  at  80-95  km  altitudes  emitted  in  a  thin  cylindrical  (doughnut-like)  shell  expanding 
to  radial  distances  up  to  150  km  [see  Inan  et  al.,  1996b,  and  references  cited  therein]  and 
referred  to  as  “elves’'  [Fukunishi  et  al.,  1996],  In  terms  of  calculations  of  the  quasi-static 
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+CG  -CG  1C  Ml  M2 


Fig.  2.1.  Thundercloud  charge  removal  models.  Positive  (+CG),  negative  (-CG),  intra¬ 
cloud  (1C),  monopole  removal  (M 1 )  and  monopole  deposition  (M2)  models  of  thundercloud 
“discharges". 


electric  fields  this  approach  appears  to  be  fully  sufficient,  as  was  shown  by  comparing  results 
obtained  from  the  same  quasi-static  equations  used  in  our  model  with  the  results  of  solutions 
of  a  full  set  of  Maxwell  equations  [Baginski  and  Hodel,  1 994], 

Any  effects  associated  with  magnetic  fields  are  necessarily  left  out  in  our  formulation. 
The  effects  of  the  earth’s  magnetic  field  can  be  neglected  in  view  of  the  extremely  high 
collision  rates  that  are  in  effect  even  under  mild  heating  conditions  (Section  3.4).  Any 
magnetic  fields  generated  due  to  the  slow  temporal  variations  of  the  quasi-static  electric 
fields  are  too  small  to  have  any  effect  on  the  electron  dynamics.  Indeed,  simple  estimations 
show  that  the  electric  field  with  amplitude  ~100  V/m  which  varies  on  a  characteristic  time 
scale  ~1  ms  and  has  spatial  scale  ~100  km  leads  to  amplitudes  of  induced  magnetic  fields 
which  are  at  least  two  orders  of  magnitude  lower  than  earth’s  magnetic  field  and  thus  can 
be  ignored. 

The  thundercloud  charges  +Q,  -Q  form  a  vertical  dipole  which  is  assumed  to  develop 
over  a  time  tj.  The  positive  or  negative  charges,  or  both  are  then  “discharged”  (or  removed) 
by  decreasing  the  magnitude  of  charge  with  a  time  constant  rs,  corresponding  respectively 
to  +CG,  -CG  or  IC  lightning  discharges,  as  illustrated  in  Figure  2.2.  Experimental  data  on 
lightning  discharges  indicates  that  the  charge  removal  time  ts  may  vary  between  ~  1  ms 
and  several  tens  of  ms  [Uman,  1987,  pp.  1 24,  196].  In  Chapter  3  we  concentrate  mainly  on 
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Fig.  2.2.  Thundercloud  charge  dynamics  corresponding  to  different  models.  The  as¬ 
sumed  behavior  of  the  source  charges  corresponding  to  +CG,  -CG,  and  IC  models. 
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CG  discharges  involving  relatively  rapid  removal  of  charge,  namely  cases  for  which  rs  ~1 
ms.  Slower  discharges  are  considered  in  Chapter  4. 

A  cylindrical  coordinate  system  ( r ,  4>,  z)  is  used  with  the  2  axis  representing  altitude,  as 
illustrated  in  Figure  2.3.  The  ground  and  ionospheric  boundaries  (taken  to  be  at  c=90  km) 
are  assumed  to  be  perfectly  conducting  and  the  entire  system  is  taken  to  be  cylindrically 
symmetric  about  the  z  axis.  Two  different  types  of  boundary  conditions  can  be  utilized  at 
the  cylindrical  boundary  (can  be  taken,  for  example,  at  r=60  km).  If  we  assume  a  perfect 
conductor  at  r=60  km,  we  have  dp/dz  =  0  (or  tp  =0),  where  p  is  the  electrostatic  potential. 
This  condition  was  used  by  Illingworth  [1972]  and  essentially  amounts  to  considering  the 
solution  of  the  electrostatic  fields  in  a  “tin-can”.  Alternatively,  we  can  assume  that  0p/dr=0, 
an  assumption  adopted  by  Tzur  and  Rohle  [1985].  Both  of  the  boundary  conditions  are 
rather  artificial  and  lead  to  unphysical  modification  of  the  electric  fields  near  the  boundary. 
However,  the  temporal  and  spatial  evolution  of  the  electrostatic  field  system  does  not  depend 
significantly  on  these  edge  fields  if  the  boundaries  are  chosen  to  be  far  enough  from  the 
charge  sources.  The  distance  r=60  km  used  in  our  calculations  represents  a  trade  off  between 
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the  resolution  which  is  possible  with  available  computer  resources  and  the  accuracy  of  the 
electric  field  calculation.  As  an  example,  the  presence  of  the  conducting  boundary  at  r=60 
km  leads  to  an  error  of  less  than  ~10  %  on  the  magnitude  of  the  electric  field  at  r=50  km 
and  2=10  km. 


Fig.  2.3.  Coordinate  system  used  in  calculations. 


The  electrostatic  field  E,  given  by  E  =  — V<p,  the  charge  density  p,  and  the  conduction 
current  .7  =  aE  are  calculated  using  the  following  system  of  equations: 


d(p  +  p.s) 

ot 


+  V  •  (J  +  J g)  =  0 


(2.1) 


v  •  E  =  ( p  +  ps)/Ec 


(2.2) 
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where  ps  and  Js  are  the  thundercloud  source  charge  and  current  which  satisfy  the  equation 
dps/dt  =  —V  •  J's,  where  V  •  J's  =  V  •  Js  +  psa/e0,  where  the  effective  source  current 
Jg  is  introduced  so  as  to  compensate  any  change  in  the  thundercloud  charge  ps  due  to  the 
conduction  current.  Equation  (2.1)  can  then  be  written  in  the  form: 

jfj  +  V(t  •  E  +  pa/e„  =  0  (2.3) 

If  in  (2.3)  one  simply  assumes  dps/dt  =  —  V  •  «7S,  (i.e.,  Js  -  J's) ,  it  is  necessary  that  pa/e0 
be  substituted  with  (p  +  ps)a/£0.  In  this  case,  if  one  ignores  the  effects  of  the  conductivity 
gradients,  the  quasi-stationary  solution  of  equation  (2.3)  is  simply  p  =  —  ps,  so  that  the 
source  charge  is  completely  compensated  by  charge  induced  in  the  conducting  medium 
over  a  time  scale  ~e„/o.  Strictly  speaking,  ps  is  a  function  of  time  [i.e.,  ps  =  p,,(()]  so 
that  the  quasi-stationary  solution  p  =  ps  has  physical  meaning  only  in  the  case  when  the 
characteristic  time  scale  of  ps(t)  is  much  longer  than  e0/ a.  As  can  be  seen  from  (2. 1 )  and 
(2.2),  the  introduction  of  J's  allows  us  to  specify  the  charge  dynamics  inside  the  cloud  as  an 
externally  determined  function  of  time. 

Physically,  the  current  J's  as  a  source  of  ps  is  associated  with  the  process  of  separation 
of  charges  inside  the  cloud  and  is  directed  opposite  to  the  resulting  electric  field  [Wilson. 
1956],  The  current  J's  may  be  related  to  the  small  and  light  positively  charged  ice  splinters, 
which  are  driven  upward  by  a  strong  stream  of  convecting  air,  and  heavy  negatively  charged 
hail  particles  which  are  driven  downward  by  gravity  [Malan,  1963,  p.  76;  Uman,  1987, 
p.  65].  We  note  that  the  model  of  thundercloud  charges  used  in  this  dissertation  is  rather 
simplified  in  comparison  with  the  actual  vertical  distribution  of  charge  density  in  mesoscale 
convective  systems  (MCS)  producing  sprites,  which  may  sometimes  have  up  to  six  layers 
of  charge  with  complicated  spatial  structure  [e.g.,  Marshall  and  Rust,  1993;  Rust  et  al., 
1996].  However,  as  we  show  later,  physical  results  of  our  model  depend  mostly  on  the 
magnitude  and  altitude  of  the  charge  removed  by  the  cloud-to-ground  lightning  discharge 
and  are  essentially  independent  of  the  details  of  the  initial  configuration  of  charges  with 
different  polarities  in  the  thundercloud. 
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We  assume  that  ps  =  p+  +  p_  and  that  the  separated  dipole  charges  (i.e.,  +Q  or  —Q) 
are  represented  as  distributed  charge  densities  p+(r,  t)  or  p_(r,  t)  having  a  Gaussian  spatial 
distribution  given  by  p±(r,t )  =  p0e~^z~z±^2+rna2  with  a= 3  km,  z-=  5  km,  2+=10km,  so 
that  the  total  source  charges  are  ±Q(t)  =  f  p±(r,  t)dV .  Our  formulation  can  also  be  used  to 
consider  of  horizontally  extended  charge  distributions  similar  to  those  existing  in  mesoscale 
convective  systems  [E.  Williams,  private  communication;  Boccippio  et  al.,  1995;  Marshall 
et  al,  1996].  In  such  a  case,  the  charge  density  corresponding  to  the  thundercloud  charges 
would  have  a  Gaussian  distribution  in  the  2  direction  while  being  uniform  in  the  r  direction, 
with  its  scale  determined  by  the  radius  of  the  “disk”  of  charge.  This  disk  charge  density  is 
normalized  accordingly  so  that  we  still  have  ±Q(t)  =  f  p±{r,  t)dV. 

Equation  (2.3)  is  integrated  in  time  using  an  explicit  two  step  method  of  solution  of 
ordinary  differential  equations  [Potter,  1973,  p.  34],  At  each  time  step  the  electric  field  E  to 
be  used  in  (2.3)  is  found  by  integration  of  equation  (2.2)  via  the  method  of  cyclic  reduction 
[Bunenum,  1 969]  using  a  FORTRAN  routine  published  by  Swarztrauber  and  Sweet  [1975]. 
The  conductivity  a  in  (2.3)  is  calculated  self-consistently  by  taking  into  account  the  effect  of 
the  electric  field  on  the  electron  component  through  changes  in  mobility  (due  to  heating)  and 
electron  density  (due  to  attachment  and/or  ionization).  A  fourth  order  Runge-Kutta  method 
and  different  conservative  and  nonconservative  numerical  schemes  are  also  employed  to 
test  the  accuracy  of  results  obtained. 


2.2  LIGHTNING  CURRENT  AM)  CHARGE 

The  important  input  parameters  that  need  to  be  specified  externally  for  each  model  calcu¬ 
lation  are:  (1)  the  altitude  profiles  of  ambient  ion  conductivity  and  electron  density;  (2) 
the  magnitude  ( Qa )  of  the  charge  transferred  (e.g.,  by  a  CG  lightning  discharge)  and  the 
time  constant  (rs)  which  determines  the  rate  at  which  the  charge  is  transferred;  and  (3)  the 
altitude  of  the  charge  removed.  None  of  these  parameters  are  directly  available  in  a  typical 
experimental  context,  except  for  the  fact  that  the  National  Lightning  Detection  Network 
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(NLDN)  records  lightning  intensity  (or  peak  current  Ip  [e.g.,  I  done  el  al.,  1993])  does  imply 
some  restrictions  on  the  magnitude  of  the  removed  charge  Q  and  the  time  scale  of  the  charge 
removal  rs. 


In  the  context  of  the  quasi-electrostatic  heating  model,  the  thundercloud  charging  and 
discharging  process  is  assumed  to  consist  of  two  stages.  During  the  predischarge  stage  with 
duration  r/,  charge  Q0  slowly  (hundreds  of  seconds)  accumulates  in  the  thundercloud.  The 
second  stage  is  the  lightning  discharge,  during  which  the  thundercloud  charge  is  quickly  and 
completely  removed  from  the  thundercloud  over  a  time  ts,  usually  of  the  order  of  several 
ms.  Mathematically,  the  continuous  thundercloud  charge  dynamics  can  be  represented  in 
the  form: 


Q(t)  —  Qo 


tanh(£/ry) 

tanh(l) 


0  <  t  <  Tf 


cm  =  Qo 


_  tanh(f/T„) 
tanh(l) 


Tf  <  t  <  Tf  +  Ts 


Q(t)  =  0  ,  Tf  +  TS  <  t 


Note  that  the  choice  of  the  functional  variation  as  tanh(-)  is  rather  arbitrary  but  is  not 
critical  for  the  physics  of  the  phenomena  modeled.  In  addition,  in  Chapter  4,  we  will 
consider  a  more  complicated  charge  removal  scenario  as  illustrated  in  Figure  4.1e. 

For  this  particular  assumed  functional  dependence  of  Q(t)  the  amplitude  of  the  peak 
current  at  t  =  Tf  can  be  calculated  as 

_  dQ  _  Qq 

p  dt  rtanh(l)  ^ 

Having  calculated  7p,  the  value  of  the  peak  electric  field  at  a  distance  of  100  km  from  the 
lightning  discharge  can  be  found  as  [Orville,  1991] 

I_pV 

2tDe0c? 


E  ioo  = 
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where  c  is  speed  of  light  in  free  space,  eQ  is  dielectric  permittivity  of  free  space,  .0=10* 
m,  and  v  =  1 .5  x  10x  m/s. 

For  example,  JE,]oo=’75  V/m  corresponds  in  this  model  to  a  peak  current  of  Ip=210  kA  and 
a  total  charge  removed  in  lightning  of  Q0= 200  C,  if  we  assume  rs=l  ms,  a  reasonable  value 
based  on  experimental  data  [ Uman ,  1987,  pp.  124,  199].  We  note  that  the  duration  of  the 
lightning  discharge  rs  as  well  as  the  functional  time  dependence  of  the  charge  removal  Q(t ) 
as  given  above  are  not  available  from  measurements  provided  by  NLDN.  Nevertheless,  these 
parameters  constitute  critically  important  inputs  for  modeling  the  corresponding  heating  and 
ionization  changes  in  the  lower  ionosphere. 

Another  important  input  to  the  model  is  the  ambient  ion  conductivity  profile,  which 
determines  the  relaxation  time  for  the  quasi-electrostatic  field  and  thus  the  amount  and 
duration  of  heating  for  given  values  of  Q0  and  rs.  Detailed  description  of  the  conductivity 
model  is  given  in  the  following  sections. 

2.3  AMBIENT  CONDITIONS  IN  THE  ATMOSPHERE  AND  LOWER 
IONOSPHERE 

For  the  various  model  calculations  in  this  thesis,  we  consider  three  different  models  of 
ambient  electron  number  density  (Figure  2.4a).  These  models  of  electron  density  were 
used  in  previous  studies  of  subionospheric  VLF  propagation  in  the  presence  of  localized 
ionospheric  disturbances  caused  by  lightning  induced  electron  precipitation  [e.g.,  Pasko  and 
I nan ,  1994],  They  were  employed  in  this  dissertation  to  represent  a  wide  range  of  possible 
variability  of  the  electron  density  in  the  nighttime  lower  ionosphere. 

For  the  ambient  ion  conductivity,  we  consider  three  different  altitude  profiles  as  illustrated 
in  Figure  2.4b.  For  illustrative  purposes  the  ambient  conductivity  profiles  in  Figure  2.4b 
also  include  the  electron  component  ae  which  is  given  by 


Oe  —  qeNeiie 


(2.4) 
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Fig.  2.4,  Ambient  electron  density  and  ion  conductivity.  Three  altitude  profiles  of  am- 
hienl  eleelron  number  density  (a),  and  ion  conductivity  (b)  used  in  the  various  calculations 
reported  in  this  thesis. 


where  qc  is  the  electronic  charge,  Ne  is  the  number  density  of  electrons,  and  //c  is  the  electron 
mobility  and  which  is  typically  dominant  at  altitudes  above  ~60  km.  The  total  conductivity 
profiles  shown  in  Figure  2.4b  were  calculated  with  the  electron  component  (dominant  at  >60 
km  altitude)  determined  by  the  electron  density  profile  1  shown  in  Figure  2.4a,  and  electron 
mobility  for  essentially  cold  electrons  of  //e  =  1 ,36Ar0/Arm2/V/s,  where  Na  =  2.688  x  If)23 
m“  3,  and  N  is  number  density  of  air  molecules  (see  Section  2.4. 1  for  further  details). 

A  altitudes  <60  km,  the  total  conductivity  is  dominated  by  the  ion  conductivity.  Profile 
A  (rr^  =  5  x  10“ 14  e~//6km  S/m)  is  taken  from  [Dejnakarintra  and  Park ,  1974],  while  Profile 
C  (ac  =  6  x  10-13  e2/1 lkm  S/m)  is  derived  from  experimental  measurements  of  Holzworth 
et  al.  [  1 985].  The  intermediate  profile  B  is  calculated  on  the  basis  of  the  known  mobilit)  of 
ions  m  as  a  function  of  altitude  (//,  ~  23Na/N  m2/V/s)  [Davies,  1983])  and  assuming  the 
number  density  of  ions  to  be  ~  103  cm-3  [Reid  et  al.,  1979],  We  note  that  A,  B  and  C  are 
used  to  mark  ambient  ion  conductivity  profiles.  In  further  calculations  we  will  use  different 
combinations  of  ions  conductivity  profiles  (A,  B,  or  C)  and  electron  conductivities  defined 
by  the  electron  density  profiles  (1,2,  or  3,  as  shown  in  Figure  2.4a). 

We  neglect  the  effects  of  the  magnetic  field  of  the  earth  on  the  conductivity  since,  for  the 
range  of  parameters  considered,  the  ionospheric  plasma  is  highly  collision  dominated  (see 
Section  3.4  for  detailed  justification  of  this  assumption). 
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2.4  SELFCONSISTENT  EVALUATION  OF  CONDUCTIVITY 

In  a  weakly  ionized  collisional  medium  such  as  the  lower  ionosphere  and  the  mesosphere, 
both  the  electron  mobility  and  the  electron  density  depend  nonlinearly  on  the  electric  field. 
The  heating  of  the  ambient  electrons  by  the  intense  fields  leads  to  substantial  modification 
of  the  atmospheric  conductivity  at  altitudes  >60  km.  The  modified  conductivity  in  turn 
determines  the  penetration  and  relaxation  of  the  electric  field,  so  that  it  is  necessary  to  self 
consistently  determine  the  evolution  of  the  conductivity  as  we  solve  for  the  evolution  of 
the  quasi-static  field  and  the  associated  consequences  (e.g.,  optical  emissions).  This  self 
consistent  solution  is  one  aspect  of  our  work  that  was  left  out  in  earlier  formulations  of 
penetration  of  thundercloud  fields  to  the  ionosphere  [e.g.,  Dejnakarintra  and  Park,  1974; 
Baginski  et  al.,  1988],  in  which  the  conductivity  was  simply  assumed  to  remain  equal  to 
the  assumed  ambient  profiles.  Note  that  the  heating  of  electrons  only  modifies  conductivity 
at  >~60  km  altitudes  where  the  electron  component  of  the  conductivity  is  dominant.  At 
altitudes  below  60  km,  where  ions  are  dominant,  the  conductivity  remains  equal  to  the 
ambient  levels,  since  the  ions  are  not  significantly  heated  by  the  electric  field. 

2.4.1  Electron  Mobility 

The  electron  mobility  fie  is  a  complicated  function  of  the  electric  field  as  determined  by  the 
various  loss  processes  each  having  different  cross  sections  [see  Taranenko  et  al.,  1993a,b 
and  references  therein]  In  this  dissertation,  we  adopt  a  functional  dependence  of  /xe  on  the 
magnitude  of  the  electric  field  E  which  is  derived  from  experimental  data  [ Davies ,  1983; 
Hegerberg  and  Reid,  1980].  This  function  is  shown  in  Figure  2.5  and  can  be  expressed  in 
the  following  analytical  form: 

2 

log (/xc7V)  =  ENo/N  >  1.62  x  103  V/m  (2.5) 

2=0 


fieN  =  1.36AT0,  EN0/N  <  1.62  x  103  V/m 
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where  x  =  \og{E/N),  aD  =  50.970,  a\  =  3.0260,  and  <22  =  8.4733  x  10~2.  The  analytical 
expression  (2.5)  can  be  used  forcomputationally  efficient  calculation  of  the  electron  mobility 
for  any  combination  of  electric  field  and  number  density  of  air  molecules  (i.e.,  altitude). 
This  analytical  method  provides  results  which  compare  well  with  those  obtained  from  fully 
kinetic  formulations,  as  discussed  in  further  detail  in  Appendix  A. 


Fig.  2.5.  Electron  mobility.  The  electron  mobility  /rf  is  shown  as  a  function  of  the  electric 
field  E  and  the  molecular  number  density  N  of  atmospheric  gas. 


Another  parameter  which  is  often  useful  in  specifying  the  degree  of  electron  heating  is 
the  effective  frequency  of  electron  collisions  with  ambient  neutral  gas.  This  parameter, 
denoted  t'eff.  is  given  as  a  simple  function  of  the  electron  mobility  /rr  as  /xeff  =  qe/(fieme). 

2.4.2  Ionization  Coefficient 

In  addition  to  the  nonlinear  dependence  of  /re  on  E,  the  local  value  of  the  electron  density 
Ne  depends  highly  nonlinearly  on  the  electric  field  due  to  impact  ionization  and  dissociative 
attachment  processes.  These  processes  begin  to  become  significant  as  the  electric  field 
exceeds  a  certain  threshold  and  eventually  lead  to  “dielectric  breakdown”  of  the  atmospheric 
gas. 
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For  quantitative  description  of  atmospheric  and  ionospheric  air  breakdown  we  use  the 
ionization  coefficient  Vi  given  in  [Papadopoulos  et  ai,  1993],  derived  from  recent  results 
of  computer  simulations  and  experimental  measurements,  which  for  our  quasi-static  case 
(cu  — »  0,  where  u  is  the  angular  frequency  with  which  the  electric  field  is  assumed  to  vary) 
is  illustrated  in  Figure  2.6  and  which  can  be  expressed  in  the  following  analytical  form 
[Papadopoulos  el  al.,  1993]: 


Vi  =  7.6  x  10_13Nx2/(x)e_4'7(l/x_1), 


(2.6) 


where 


f(x)  = 


1  +  6.3e-2'6/x 

L5 


x  -  E/Ek,  and  Ek  -  3.2  x  \06N/No  V/m  is  the  characteristic  air  breakdown  field  shown 
for  reference  in  Figure  2.6  by  a  dashed  line. 


Fig.  2.6.  Ionization  and  dissociative  attachment.  The  ionization  v,  and  dissociative 
attachment  ua  coefficients  are  shown  as  a  function  of  the  electric  field  E  and  the  molecular 
number  density  N  of  the  atmospheric  gas. 


It  should  be  noted  that  due  to  the  relatively  slow  variation  of  the  electric  field  the  con¬ 
ductivity  in  our  model  is  assumed  to  be  independent  of  frequency.  Similarly,  the  ionization 
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coefficient,  although  originally  derived  for  radio  frequency  (RF)  fields  of  certain  definable 
frequency  a ;  [ Papadopoulos  et  al,  1993],  is  used  in  this  dissertation  to  be  in  DC  limit 
(u;  — >  0).  Both  assumptions  are  well  justified  since  for  quasi-static  electric  fields  u>  <~  1 0^ 
rad/s  and  always  remains  much  less  than  z/eff  and  frequencies  in  RF  range  [ Papadopoulos 
et  al.,  1993], 

2,4,3  Attachment  Coefficient 

In  our  quasi-electrostatic  heating  formulation,  we  can  include  the  two-body  attachment  of 
electrons  to  O2  molecules  with  dissociation  (e_  +  O;  — +  O  +  O-).  The  effective  coefficient 
//„  of  this  reaction  is  shown  in  Figure  2.6  and  is  a  function  of  electric  field  which  is  derived 
from  experimental  data  [Davies,  1983].  This  quantity  can  be  expressed  in  the  following 
analytical  form: 


isa  =  —  2_^aj.T’  (2.7) 

0  »=0 

where  x  =  ENa/N,  a0  =  —2.41  x  10*,  «i  =  21  1.92,  and  ai  =  —3.545  x  10-?.  The 


attachment  coefficient  ua  given  in  (2.7)  is  an  effective  coefficient  which  is  a  function  of 
the  true  attachment  coefficient  (for  O-  formation),  the  charge  transfer  coefficient,  and  the 
detachment  coefficient  [Davies,  1983], 

At  very  low  values  of  //,  and  ua  (in  the  region  of  electric  fields  E  <  E /,.)  the  value  of  the 
attachment  coefficient  ua  as  given  by  (2.7)  becomes  smaller  than  the  ionization  coefficient 
u}  as  given  by  (2.6)  and  leads  to  unphysical  results.  For  such  low  values  of  uaN0/N ,  which 
correspond  to  EN0/N  <  1.628  x  106  V/m,  the  value  of  the  attachment  coefficient  can  be 
found  from  kinetic  calculations  [Taranenko  et  at.,  1993b]  as: 


(2.8) 


where  x  =  log(£'AT0/A'T),  and  the  series  of  coefficients  a0  =  —  1 073.8 , 0.1  =  465.99.  at  = 
—66.867,  and  03  =  3.1970. 
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In  previous  work  on  swarm  experiments  [Davies,  1983]  it  was  recognized  that  the  role 
of  dissociative  attachment  of  electrons  to  O2  (the  fastest  process  of  loss  of  electrons  under 
heated  conditions  which  may  lead  to  significant  modifications  of  the  lower  ionosphere 
[Taranenko  el  al.,  1993a;  Inan  et  al.,  1996b])  is  difficult  to  assess  since  the  measured  rates 
for  associative  detachment  and  for  collisional  detachment  are  known  only  within  one  to 
three  orders  of  magnitude,  respectively.  A  detailed  discussion  of  the  importance  of  the 
detachment  processes  in  the  analysis  of  swarm  experiments  is  given  in  [Wen  and  Wetzer, 
1988], 

2.4.4  Electron  Density 

The  ionization  vt  (2.6)  and  attachment  ua  (2.7),  (2.8)  coefficients  are  used  to  calculate  the 
temporal  dynamics  of  the  electron  density  using  following  equation: 

— 7T  =  fa  -  Va)Ne  -  rtJVe2  (2.9) 

at 

where  a  is  an  effective  recombination  coefficient  which  in  most  of  calculations  in  this 
dissertation  is  assumed  to  be  zero. 

As  will  be  discussed  further  in  Chapters  4  and  5  ,  both  sprites  and  jets  can  be  closely 
associated  with  streamer  type  plasma  processes,  which  can  be  briefly  described  as  a  special 
category  of  nonlinear  space  charge  waves  [e.g.,  Vitello  et  al.,  1994].  To  represent  streamers 
with  our  model  which  is  essentially  based  on  one  fluid  system  of  equations  (2.1)  and  ( 2.2)  we 
rely  upon  different  assumptions  derived  from  the  results  of  two-fluid  simulations  of  streamer 
processes  [e.g.,  Vitello  et  al.,  1994].  These  assumptions  most  concern  the  proper  treatment 
of  the  steep  boundary  at  the  streamer  head,  and  are  implemented  in  our  model  as  different 
types  of  limitations  on  the  growth  of  the  electron  density  due  to  breakdown  ionization. 
One  of  these  limitations  is  implemented  through  the  effective  recombination  coefficient  a 
in  (2.9)  which  in  this  context  has  straightforward  physical  meaning.  The  recombination 
effect  is  used  for  the  modeling  of  the  downward  propagation  of  channels  of  the  breakdown 
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ionization  to  altitudes  <50  km  as  discussed  in  Chapter  4.  Estimates  indicate  that  at  higher 
altitudes  three  body  attachment  processes  as  well  as  the  bulk  of  other  chemical  reactions 
in  the  lower  ionospheric  D-region  (including  different  recombination  processes)  have  time 
scales  >1  sec  [e.g Glukhov  et  al.,  1992 ;  Pasko  and  Inan,  1994]  and  can  be  neglected  on  the 
millisecond  time  scales  of  interest  in  this  dissertation.  Another  limitation  which  is  used  for 
modeling  of  blue  jets  simply  assumes  the  constancy  of  the  dielectric  relaxation  time  behind 
the  streamer  front.  This  assumption  will  be  additionally  discussed  in  Chapter  5. 

In  summary,  in  equation  (2. 1 )  the  conductivity  a  is  calculated  as  a  sum  of  the  ambient  ion 
component  <r,  which  remains  constant  in  time  and  the  electron  component  ae  (2.4)  which 
is  calculated  self-consistently  by  taking  into  account  the  effects  of  the  electric  field  through 
changes  in  mobility  fic  (2.5)  due  to  heating,  and  changes  in  the  number  density  of  electrons 
Nc  (2.9)  due  to  ionization  and  attachment  processes. 

2.5  OPTICAL  EMISSIONS 

Once  the  temporal  and  spatial  evolution  of  the  electric  field  is  determined  via  the  selfconsis- 
tent  solution  of  (2. 1 )  and  (2.2),  the  intensities  of  the  various  optical  emissions  are  determined 
on  the  basis  of  kinetic  solutions  for  the  electron  distribution  function  and  the  excitation  co¬ 
efficients  of  particular  bands.  We  consider  optical  emissions  from  the  1  st  {D^Y\g  —*  A  ’£■„) 
and  2nd  (C3FIU  — ♦  B2Y\q)  positive  bands  of  Nt,  the  Meinel  (A2FI  — *  X2Hq)  and  1st  nega¬ 
tive  bands  of  Nj  (B21.+  —>  „Y2Z+),  and  the  1st  negative  band  of  0}  (fr4X^  — *  a4T\u)  which 
have  short  lifetimes  and  are  expected  to  produce  the  most  intense  optical  output  during 
impulsive  heating  of  the  lower  ionospheric  electrons  [Taranenko  et  al.,  1993a,b]. 

Analysis  of  the  excitation  coefficients  14.  corresponding  to  these  bands  obtained  from 
kinetic  calculations  [Taranenko  et  al.,  1993a,b]  for  a  series  of  electric  field  values  and  air 
densities  shows  that  relatively  simple  and  accurate  analytical  approximations  similar  to 
those  developed  for  the  electron  mobility  (2.5),  ionization  (2.6)  and  attachment  (2.7),  (2.8) 
coefficients  can  be  derived  (see  Appendix  A.l  ).  The  results  are  shown  in  Figure  2.7  and 
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TABLE  2. 1 .  Approximation  Coefficients. 


Band 

do 

a  1 

a2 

G3 

1  St  pos.  N2 

-1301.0 

563.03 

-80.715 

3.8647 

2nd  pos.  N2 

-1877.1 

814.70 

-117.46 

5.6554 

1st  neg.  N* 
1st  neg.  O2 

-1760.0 

724.7 

-99.549 

4.5862 

-1802.4 

750.91 

-104.28 

4.8508 

Meinel  N* 

-2061.1 

870.25 

-122.43 

5.7668 

for  the  optical  bands  considered  can  be  expressed  in  a  general  form: 


?=0 


(2.10) 


where  x  =  \og(EN0/N),  and  the  series  of  coefficients  a0, 01,02,03  corresponding  to  the 
different  optical  bands  are  given  in  Table  2.1. 


Fig.  2.7.  Optical  excitation  coefficients.  The  excitation  coefficients  Vk  of  different  optical 
bands  arc  shown  as  a  function  of  electric  field  E  and  molecular  number  density  N  of 
atmospheric  gas. 
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The  intensity  of  each  optical  line  in  Rayleighs  is  given  by  the  expression  [ Chamberlain , 
1978,  p.  213]: 

4=I(T6  jAknkdl,  (2.11) 

where  n *  is  the  number  density  of  excited  particles  in  state  k,  A jt  is  the  radiation  transition 
rate,  and  the  integral  is  taken  along  L,  representing  the  horizontal  line  of  sight.  In  our 
calculations,  we  do  not  take  into  account  the  effects  of  radiative  transfer  between  the  altitude 
of  the  emission  and  the  ground  which  might  be  important  for  the  interpretation  of  the  total 
intensity  of  optical  emissions  and  their  fine  spectral  features  (i.e.,  the  strong  attenuation  of 
the  (3, 1 )  component  of  the  first  positive  band  of  N2  due  to  the  absorption  by  O2  at  762  nm 
over  slant  paths  [Mende  et  al.,  1995;  Hampton  et  al„  1996]).  The  quantity  n*.  is  governed 
by  the  relation  [Sipler  and Biondi,  1972]: 

-p~  = — ~  +  ^  timAm  +  Uf;Nc  (2.12) 

at  ti. 

m 

where  r^.  =  [zt*.  +  o:|  An2  +  c^Ao,]-1  is  the  total  lifetime  of  state  k,  Qi  and  a2  arc  the 
quenching  rates  due  to  collisions  with  N2  and  O2  molecules,  respectively,  A’n2  and  A’o, 
are  the  number  densities  of  N2  and  O2  molecules  respectively,  and  the  sum  over  the  terms 
nmAm  represents  increases  in  n*  resulting  from  cascading  from  higher  energy  states. 

The  intensity  of  the  optical  bands  in  equation  (2.1  I)  is  obtained  from  the  number  of 
excited  molecules  in  the  upper  state  n*.  which  is  given  by  equation  (2.12).  In  equation 
(2.12)  the  electron  number  density  Nc  and  the  excitation  coefficients  14.  are  the  variables 
which  evolve  in  time  and  give  rise  to  the  optical  emissions,  and  are  calculated  as  functions 
of  the  electric  field  using  equations  (2.9)  and  (2. 10),  as  previously  described. 

The  optical  emissions  from  the  1st  and  2nd  positive  bands  of  N2,  the  Meinel  and  1st 
negative  bands  of  Nt,  and  the  1st  negative  band  of  Oj  have  the  transition  rates  1 .7  x  105/s, 
2x  10"7s,  7x  l04/s,  1.4x  107s,  and  8.5 x  105/s  respectively  [ Valence-Jones ,  1974.  p.  1 19], 
Quenching  of  the  1st  positive  band  emissions  occurs  primarily  through  collisions  with  N2 
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molecules,  with  ctj  =  1 0— 1 1  cmVs  [Valence-Jones,  1974,  p.  119],  Quenching  of  the  2nd 
positive  band  occurs  primarily  through  collisions  with  O2  molecules  with  a2=3xlO-10 
cirr/s  [Valence-Jones,  1974,  p.  1 19].  For  the  quenching  of  the  Meinel  and  1st  negative 
bands  of  N£  we  have  aj  =  5x  10~10  cm  Vs  and  a]  =  4  x  10-,°  cmVs  respectively,  and  for 
the  1st  negative  band  of  Ot  we  have  a\  =  2  x  10_,°  ctrTVs  [Valence -Jones,  1974,  p.  1 19]. 


2.6  LIMITATIONS  OF  ANALYTICAL  MODELS 

The  simple  expressions  for  mobility  (2.5),  ionization  (2.6),  attachment  (2.7),  (2.8),  and 
optical  excitation  (2. 1 0)  coefficients  are  supported  by  experimental  measurements  as  well  as 
by  results  of  kinetic  simulations  (see  Appendix  A.  1  for  details).  These  analytical  expressions 
are  useful  from  a  computational  point  of  view  since  they  provide  fast  and  reasonably  accurate 
solutions  for  a  wide  range  of  input  parameters.  However,  these  models  are  only  valid  only 
for  EN0/N  <  1.6  x  107  V/m  since  for  fields  above  this  value  the  mean  electron  energy 
approaches  the  ionization  energy  and  the  ionization  process  becomes  nonstationary  (see 
Papadopoulos  et  al.  [1993]  and  references  cited  therein).  The  criteria  EN0/N  <  1.6  x  Iff 
V/m  is  equivalent  to  the  condition  E/E^  <  5  given  by  Papadopoulos  et  al.  [1993].  At 
70  km  altitude,  this  criteria  requires  that  E  <1200  V/m,  a  condition  comfortably  satisfied 
for  all  the  cases  considered  in  this  work.  In  addition,  the  above  expressions  cannot  be  used 
on  time  scales  less  than  ~10  ps  (this  estimate  is  valid  for  altitudes  ~  80-90  km)  which 
is  the  time  required  for  the  establishment  of  the  stationary  electron  distribution  function 
[Taranenko  et  al.,  1 993a]. 


3 

Heating,  Ionization  and  Optical  Flashes 
in  the  Mesosphere 


In  this  chapter,  we  present  results  obtained  using  the  quasi-electrostatic  heating  model 
described  in  the  previous  chapter  for  a  variety  of  ambient  atmospheric/lower  ionospheric 
conditions. 


3.1  PHYSICAL  MECHANISM 

The  basic  mechanism  of  the  penetration  of  the  thundercloud  electric  fields  to  the  higher 
altitude  regions  is  illustrated  in  Figure  3.1.  As  the  thundercloud  charges  slowly  build  up 
before  a  lightning  discharge,  high  altitude  regions  are  shielded  from  the  quasi-electrostatic 
fields  of  the  thundercloud  charges  by  space  charge  induced  in  the  conducting  atmosphere 
at  lower  altitudes.  The  appearance  of  this  shielding  charge  is  a  consequence  of  the  finite 
vertical  conductivity  gradient  of  the  atmosphere  above  the  thundercloud.  When  one  of  the 
thundercloud  charges  (e.g.,  the  positive  one  as  shown  in  Figure  3. 1 )  is  quickly  removed  by 
a  lightning  discharge,  the  remaining  charges  of  opposite  sign  in  and  above  the  thundercloud 
produce  a  large  quasi-electrostatic  field  that  appears  at  all  altitudes  above  the  thundercloud, 
and  endures  for  a  time  equal  to  approximately  (see  discussion  in  Section  3.10)  the  local 
relaxation  time  (rr  =  £0/v)  at  each  altitude.  This  large  electric  field  can  be  thought  of 
as  being  the  difference  between  the  electrostatic  field  dictated  (via  Coulomb’s  law)  by  the 
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dipole  configuration  of  thundercloud  charges  and  the  polarization  charge  which  are  in  effect 
before  the  discharge,  and  that  required  by  the  combination  of  the  single  thundercloud  charge 
remaining  after  the  discharge  and  the  polarization  charge.  These  temporarily  existing  electric 
fields  lead  to  the  heating  of  ambient  electrons  and  the  generation  of  ionization  changes  and 
optical  emissions. 

In  cases  of  more  complex  charge  distributions  in  the  thundercloud  which  sometimes 
involve  up  to  six  charge  layers  in  the  vertical  direction  [e.g.,  Marshall  and  Rust,  1993], 
each  of  the  charge  centers  can  be  viewed  as  generating  its  own  polarization  charge  in 
and  above  thundercloud,  and  the  resultant  configuration  of  the  electric  field  and  charge 
density  can  be  obtained  by  using  principle  of  superposition.  This  consideration  is  helpful 
in  visualization  of  the  fact  that  the  electric  field  appearing  at  mesospheric  altitudes  after  the 
charge  removal  by  cloud-to-ground  lightning  discharge  is  defined  mostly  by  the  absolute 
value  and  altitude  of  the  removed  charge  and  is  essentially  independent  of  the  complexity 
of  the  charge  configuration  in  the  cloud.  The  charge  removal  can  also  be  viewed  as  the 
“placement”  of  an  identical  charge  of  opposite  sign.  The  initial  field  above  the  cloud  is 
simply  the  free  space  field  due  to  the  “newly  placed"  charge  and  its  image  in  the  ground 
which  is  assumed  to  be  perfectly  conducting.  This  effect  will  be  further  illustrated  in  the 
context  of  the  results  of  our  calculations  (Section  3.8). 


3.2  MODEL  PARAMETERS 

In  most  of  the  calculations  in  this  chapter  we  assume  the  lightning  discharge  duration  to 
be  ~1  ms.  We  evaluate  the  electron  density  changes  due  to  the  dissociative  attachment  of 
electrons  to  O2  and  the  impact  ionization  of  O2  and  N2,  as  well  as  the  dynamics  and  the  two- 
dimensional  distribution  of  the  optical  emission  intensity  of  the  1  st  and  2nd  positive  bands  of 
N2,  Meinel  and  1st  negative  bands  of  N|,  and  the  1st  negative  band  of  O2,  caused  by  quasi¬ 
electrostatic  heating  and  ionization  in  the  lower  ionosphere.  We  examine  the  dependence 
of  the  heating  and  ionization  upon  the  altitude  profile  of  ion  conductivity  and  the  number 
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Sprites 


Before +CG  -  1-10  ms  after  +CG  -  1  sec  after +CG 


Time 


Fig.  3.1.  Physical  mechanism  of  mesospheric  optical  flashes.  Illustration  of  the  mecha¬ 
nism  of  penetration  of  large  quasi-electrostatic  thundercloud  fields  to  mesospheric  altitudes, 
leading  to  electron  heating,  breakdown  ionization  and  excitation  of  optical  emissions. 


density  of  electrons.  Some  effects  associated  with  the  duration  of  the  lightning  discharge  are 
also  discussed,  although  further  discussion  of  this  subject  is  provided  in  Chapter  4.  We  also 
discuss  the  dependence  of  the  quasi-electrostatic  heating  on  the  lightning  discharge  type 
(i.c.,  intracloud  versus  positive  or  negative  cloud-to-ground)  as  well  as  on  the  thundercloud 
charge  geometry.  Finally,  we  analyze  the  relaxation  properties  of  the  electric  field  in  the 
conducting  medium. 

In  all  of  the  calculations  leading  to  the  results  presented  below  in  Sections  3.3,  3.4  and 
3.5  ,  we  use  the  ambient  electron  density  profile  1  (Figure  2.4a),  the  ambient  ion  conductivity 
profile  A  (Figure  2.4b),  and  a  lightning  discharge  duration  of  r,=  l  ms.  In  these  sections 
we  consider  the  +CG  type  of  lightning  discharge  (see  Figures  2.1,  2.2).  In  Sections 
3.3,  3.4  and  3.5  ,  the  distribution  of  thundercloud  charges  is  assumed  to  be  spherically 
symmetric  and  Gaussian  with  a  spatial  scale  a  =3  km.  The  effects  produced  for  different 
values  of  above  mentioned  parameters  are  investigated  in  Sections  3.6  through  3.9.  In  all 
cases,  we  assume  that  the  thundercloud  charge  is  removed  by  a  lightning  discharge  from 
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10  km  altitude.  Results  obtained  for  given  charge  values  at  10  km  altitude  are  valid  for 
different  combinations  of  charge  value  and  altitude  as  long  as  the  total  thundercloud  dipole 
moment  is  kept  constant.  From  this  point  of  view,  the  variation  of  the  altitude  of  the  charge 
is  equivalent  (in  terms  of  the  resultant  physical  effects)  to  the  variation  of  the  value  of  charge 
at  a  given  altitude,  and  vise  versa.  For  instance,  a  1 00  C  charge  removed  from  1 0  km  altitude 
produces  approximately  the  same  effects  in  the  lower  ionosphere  as  a  50  C  charge  removed 
from  20  km  altitude. 


3.3  ELECTRIC  FIELD  AND  CHARGE  DENSITY 

Figures  3.2  and  3.3  show  the  results  of  model  calculations  for  the  case  of  dipole  thun¬ 
dercloud  charges  of  ±200  C.  These  results  are  similar  to  the  monopole  case  considered  in 
[Pasko  et  al.,  1995]  and  marked  Ml  in  Figure  2.1.  The  slow  {jf= 0.5  s)  build-up  of  the 
charges  leads  to  the  formation  of  small  but  finite  quiescent  quasi-electrostatic  fields,  which 
are  similar  in  nature  to  the  fair  weather  electric  field  produced  by  the  charges  maintained 
on  the  earth  [e.g.,  Urnan,  1974;  Ogawa,  1985].  The  removal  of  a  positive  charge  in  rs=l 
ms  leads  to  the  appearance  of  electric  fields  of  up  to  ~  100  V/m  at  ionospheric  altitudes. 
Figure  3.3  shows  the  spatial  distribution  of  the  charge  density  f>  and  the  electric  field  E  just 
before  ( t  =0.5  s)  and  right  after  ( t  =0.501  s)  the  cloud-to-ground  lightning  discharge  and 
at  /  =1  s.  The  physical  causes  of  the  “explosive”  appearance  of  the  electric  field  at  higher 
altitudes  at  f  =  0.501  s  (as  already  briefly  discussed  in  Section  3.1)  can  be  understood  as 
follows: 

As  the  thundercloud  charges  accumulate  before  the  lightning  discharge,  free  charges  in 
the  mesosphere  and  lower  ionosphere  redistribute  (Figure  3.3,  top  left  panel)  to  create  an 
electric  field  E  approximately  equal  and  opposite  to  the  “applied”  electric  field  E  due  to  the 
large  thundercloud  charges.  Once  the  thundercloud  charge  of  any  polarity  is  removed,  only 
the  electric  field  E  due  to  the  charge  remaining  in  thundercloud  and  the  space  charge  above 
it  remains,  and  this  large  field  endures  for  a  time  determined  approximately  (see  discussion 
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Fig.  3.2.  Temporal  variation  of  the  quasi-static  electric  field.  The  cstablishmcnl  and 
relaxation  of  the  z  component  of  the  electric  field  at  various  altitudes  at  r-0  km  for  the  +CG 
and  Ml  models  of  charge  removal  (see  Figures  2.1  and  2.2). 


in  Section  3.10)  by  the  local  relaxation  time  (rr  =  s0/a)  at  each  altitude.  As  we  discussed 
in  Section  3.1,  since  charge  removal  can  alternately  be  viewed  as  the  "placement"  of  an 
identical  charge  of  opposite  sign,  the  initial  field  at  ionospheric  altitudes  is  approximately 
the  free  space  field  due  to  the  “newly  placed"  charge  and  its  image  in  the  ground.  Thus, 
the  important  physical  consequences  depend  primarily  on  the  magnitude  and  altitude  of  the 
removed  charge  and  are  essentially  independent  of  the  configuration  of  all  other  charges  in 
the  system.  In  this  regard  we  note  that  +CG  amd  M  1  models  of  charge  removal  (Figure  2. 1 ) 
lead  to  essentially  the  same  electric  fields  at  ionospheric  altitudes  as  is  discussed  futher  in 
Section  3.8. 


Altitude  (km)  Altitude  (km) 
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Fig.  3.3.  Electric  field  and  charge  density.  A  cross  sectional  view  of  the  distribution  of 
the  absolute  values  of  the  charge  density  p  (left  panels)  and  the  total  electric  field  E  (right 
panels)  at  selected  instants  of  time.  The  top  panels  show  altitude  scans  at  r=0. 
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3.4  IONIZATION  AND  ATTACHMENT 

Figure  3.4  shows  the  altitude  profile  of  electron  density  Nc  and  the  conductivity  a  for 
different  values  of  removed  thundercloud  charge  Q  (Q= 50,  100,  150,  200  and  250  C) 
at  the  termination  of  a  lightning  discharge  (7=1  ms).  Corresponding  profiles  of  effective 
electron  collision  frequency  izetf,  which  provides  a  measure  of  electron  heating,  are  shown 
in  Figure  3.5. 

During  the  charge  accumulation  stage  prior  to  the  lightning  discharge,  the  electric  field 
remains  below  the  critical  breakdown  field  E 7  (see  Section  2.4.2)  everywhere  above  the 
thundercloud  (sec  top  right  panel  in  Figure  3.3).  However,  the  large  electric  field  which 
appears  after  the  removal  of  the  charge  exceeds  the  breakdown  threshold  E above  ~7() 
km  altitude,  leading  to  intense  ionization  changes  in  the  lower  ionosphere  (Figure  3.4  )  and 
optical  emissions  (Figure  3.9). 

For  higher  values  of  the  removed  thundercloud  charge  Q  we  see  from  Figure  3.4  a 
transition  from  a  regime  in  which  attachment  dominates  leading  to  electron  density  decreases 
of  several  orders  of  magnitude  below  ambient  values  to  an  ionization  regime  in  which 
dramatic  increases  in  electron  density  occur.  For  low  Q  values  such  as  50  C  (Figure  3.4) 
the  attachment  and  ionization  thresholds  are  not  reached,  but  nevertheless  a  large  change  in 
effective  electron  collision  frequency  occurs  (Figure  3.5),  indicating  significant  changes  in 
electron  conductivity  (and  thus  the  total  conductivity  at  altitudes  >60  km)  purely  due  to  the 
heating  of  the  electrons. 

The  changes  in  the  electron  density  Nc  result  in  corresponding  changes  in  conductivity  a 
(right  hand  panels  in  Figure  3.4),  although  the  resultant  conductivity  changes  are  determined 
by  the  interplay  between  changes  in  electron  density  Ne  (due  to  ionization  and  attachment) 
and  electron  mobility  fir  (due  to  heating)  (2.4).  For  example,  the  factor  of  ~10  lower 
conductivity  at  ~80  km  altitude  for  (^=200  C  versus  100  C  is  dominantly  due  to  electron 
heating,  since  the  electron  density  changes  are  relatively  small. 

The  altitude  scan  of  the  effective  electron  collision  frequency  i/cff,  modified  by  heating. 
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Fig.  3.4.  Model  results  for  electron  density  and  conductivity.  The  altitude  distribution 
of  the  electron  density  Ne  and  conductivity  cr  for  different  values  of  thundercloud  charge  Q 
at  the  moment  of  time  f=lms.  The  top  panels  show  altitude  scans  at  r=0. 
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Fig.  3.5.  Model  results  for  effective  electron  collision  frequency  The  effective  electron 
collision  frequency  vt ff  is  shown  as  function  of  altitude  for  the  same  thundercloud  charge 
values  as  in  Figure  3.4. 

is  shown  in  Figure  3.5  for  the  same  thundercloud  charge  values  as  in  Figure  3.4.  For 
purposes  of  comparison,  the  electron  gyrofrequency  is  shown  in  Figure  3.5  by  a  vertical 
dashed  line.  We  note  that  ff  >  u>h  below  ~87  km  altitude  for  all  charge  values  considered 
here.  Our  assumption  concerning  the  neglect  of  the  effects  associated  with  the  external 
geomagnetic  field  (Section  2.1)  is  justified  since  the  most  important  region  of  the  optical 
emissions  associated  with  sprites  is  well  below  this  altitude  (see  Section  3.5). 

The  results  in  Figure  3.4  indicate  that  quasi-electrostatic  fields  can  significantly  change 
the  mesospheric/lower  ionospheric  electron  density  Ne  and  conductivity  a  at  ~50  to  90  km 
altitudes  and  in  regions  with  ~50-60  km  transverse  extent.  One  of  the  manifestations  of  such 
conductivity  changes  in  experimental  data  may  be  some  early/fast  VLF  signal  perturbations 
observed  in  association  with  lightning  discharges  and  sprites  [Inan  et  ai,  1988,  1993.  1995, 
1996a], 


Often  the  time  period  between  consecutive  discharges  from  thunderstorm  centers  [see 
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Uman,  1987,  p.  19]  is  much  less  than  the  typical  D-region  relaxation  time  for  electron 
density  enhancements  (10  to  100  s  [e.g.,  Pasko  and  Inan,  1994]).  In  such  cases  quasi¬ 
electrostatic  fields  associated  with  successive  lightning  discharges  as  well  as  electromagnetic 
pulses  (EMP)  generated  by  lightning  strokes  can  lead  to  accumulation  of  the  ionization  or 
depletion  of  the  electron  density  in  the  lower  ionosphere  leading  to  the  build  up  of  ionization 
“bubbles”  or  attachment  “holes”  above  thunderstorms.  Such  an  effect  was  predicted  by  Inan 
et  al.  [1991]  and  also  discussed  by  Taranenko  [1993]  in  relation  to  heating  and  ionization 
produced  by  electromagnetic  impulses  generated  by  lightning. 


3.5  OPTICAL  EMISSIONS 

The  dynamics  of  the  spatial  distribution  of  optical  emissions  is  shown  in  Figure  3.6  for 
Q  =250  C.  The  local  optical  emissivities  of  the  first  (red)  positive  band  of  N2  excited  via 
impact  excitation  at  different  points  in  a  cylindrically  symmetric  volume  were  integrated 
along  the  horizontal  line  of  sight  to  determine  the  emission  intensities  in  Rayleighs  observed 
when  viewed  horizontally  from  a  distant  point. 

For  direct  comparison  with  video  observations  of  sprites  [ Sentman  et  al.,  1 995]  the  optical 
emissions  are  averaged  over  the  total  duration  of  sprite  luminosity  (~  1 .6  ms),  which  in  this 
case  is  substantially  less  than  the  duration  of  one  video  field  (~  16  ms).  Results  of  the 
averaging  are  shown  in  the  lower  two  panels  and  are  in  good  agreement  with  observations 
of  the  upper  part  (“head”  and  “hair”  [Sentman  et  al.,  1995])  of  the  sprites.  Following 
the  occurrence  of  the  cloud-to-ground  discharge,  the  optical  emission  initially  starts  at  high 
altitudes  and  then  moves  down,  its  motion  being  controlled  by  the  development  of  ionization 
and  the  relatively  more  rapid  relaxation  of  the  electric  field  at  higher  altitudes.  We  note  that 
since  the  total  duration  of  the  optical  emission  is  approximately  one-tenth  the  duration  of 
one  video  field,  the  maximum  intensity  of  ~105  kR  corresponds  to  ~104  kR  when  averaged 
over  a  1 6  ms  field,  a  value  larger  than  the  maximum  reported  intensity  (~500  kR)  of  sprites 
[Sentman  et  al.,  1995]  by  a  factor  ~  1 0- 1 00  when  viewed  through  the  center  of  the  heated 
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Fig.  3.6.  Space-time  dynamics  of  optical  emissions.  The  distribution  of  intensities  (in 
units  of  Rayleighs)  of  emissions  of  the  1st  positive  band  of  N2  as  a  function  of  time  for  the 
case  of  250  C  charge  removed  from  altitude  10  km  in  1  ms.  The  bottom  panel  shows  the 
corresponding  quantities  averaged  over  a  time  of  1 .6  ms. 


region  (at  ~70  km).  This  apparent  discrepancy  is  likely  due  to  the  response  of  the  optical 
instruments  used  to  observe  sprites. 

It  is  important  to  note  that  the  optical  excitation  rates  for  selected  bands  (i.e.,  1  st  positive 
band  of  N2)  are  calculated  using  integral  cross  sections  which  do  not  provide  information 
about  the  particular  wavelength  of  emitted  photons  [Taranenkn  et  al.,  1993].  Typical  video 
cameras  have  well  defined  red,  green,  and  blue  filters  [e.g.,  Wescott  et  al.,  1995a]  (an 
example  is  shown  in  upper  left  corner  of  Figure  3.7)  which  extract  only  a  narrow  part  of 
the  spectrum  in  the  visible  range  of  wavelengths.  To  determine  the  number  of  photons 
corresponding  to  the  1  st  positive  band  of  N2  which  produces,  for  example,  the  red  response 
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of  the  video  camera,  one  needs  to  integrate  the  normalized  distribution  of  emissions  in 
different  wavelengths  over  the  band  corresponding  to  the  red  filter.  This  consideration  is 
important  since  most  of  the  energy  of  the  N2  1st  positive  band  is  distributed  in  the  infrared 
region  of  the  spectrum  with  only  a  small  fraction  in  the  visible  region  covered  by  the  red 
video  filter. 
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Fig.  3.7.  Optica)  emission  intensity  as  a  function  of  wavelength.  Normalized  (to  unity) 
distribution  of  optical  emission  intensity  as  a  function  of  wavelength  (assumed  to  be  the 
same  as  in  aurora)  corresponding  to  different  optical  bands.  Top  left  corner  show's  typical 
red,  green,  and  blue  video  filters. 


The  distribution  of  the  intensities  as  a  function  of  wavelength  within  all  of  the  different 
optical  bands  corresponding  to  sprites  is  not  yet  known.  However,  on  the  basis  of  the  first 
spectral  observations  of  the  1st  positive  band  of  N2  associated  with  sprites  [Mende  et  ai, 
1995;  Hampton  et  al.,  1996]  (see  also  Section  1.2.8)  the  distribution  of  intensities  of  the 
different  wavelengths  within  this  band  appear  similar  to  that  observed  in  auroral  emissions 
[  Valance  Jones,  1974,  pp.  129- 132, 139]  (Figure  3.7).  Although  the  wavelength  distribution 
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for  the  various  components  of  other  emission  bands  is  not  yet  known,  we  assume  in  our 
further  calculations  below  that  the  distribution  within  these  bands  to  be  also  similar  to  that 
in  aurorae.  Calculations  based  on  this  assumption  indicate  the  N2  1st  and  2nd  positive 
bands  produce  the  largest  (in  comparison  with  other  optical  bands)  red  and  blue  video 
response,  respectively.  Only  1%  of  the  total  number  of  emitted  photons  corresponding  to 
the  N2  1  st  positive  band  appear  in  the  pass  band  of  the  red  filter  of  a  typical  video  camera. 
The  corresponding  scaling  factor  for  the  predominantly  blue  N2  2nd  positive  band  is  2.9%. 
Based  on  these  estimates  the  optical  intensities  shown  in  Figure  3.6  are  in  good  agreement 
with  video  observations  of  sprites  [Sentman  et  al.,  1 995].  Scaling  factors  for  selected  optical 
bands  and  different  colors  calculated  on  the  basis  of  video  filter  responses  given  in  [Wescott 
et  al.,  1995a]  are  summarized  in  Table  3.1.  These  scaling  factors  are  provided  here  for 
illustrative  purposes  to  stress  the  importance  of  this  type  of  scaling  for  the  interpretation  of 
the  experimental  video  data,  and  may  be  different  for  video  cameras  with  different  technical 
characteristics.  Furthermore,  the  scaling  factors  may  strongly  depend  on  the  assumption 
on  the  auroral-like  distribution  of  optical  emissions  as  a  function  of  wavelength  within  the 
bands  considered. 


TABLE  3. 1 .  Video  response  scaling  factors. 
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1.04xl()-- 
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Figure  3.8  shows  the  results  of  time  averaging  of  the  optical  emission  intensity  of  the 
N2  1st  positive  band  for  different  values  of  removed  thundercloud  charge.  Values  of  the 
maximum  intensity  of  the  first  positive  band  corresponding  to  Q=250,  200,  150,  100,  50  C 
are  105  kR,  104  kR,  103  kR,  100  kR,  2  kR,  respectively,  demonstrating  a  highly  nonlinear 
dependence  of  the  peak  intensity  on  the  removed  charge  value  Q,  as  expected  on  the  basis 
of  the  threshold-like  behavior  of  the  optical  excitation  cross  sections.  This  result  is  also 
consistent  with  the  experimental  fact  that  sprites  are  associated  with  only  the  most  intense 
(usually  positive  [e.g.,  Boccippio  et  al.,  1995])  lightning  discharges. 
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Fig.  3.8.  The  distribution  of  time-averaged  optical  emission  intensity.  The  distribution 
of  the  time-averaged  intensity  (in  units  of  Rayleighs)  of  emissions  of  the  1st  positive  band 
of  N2  is  shown  on  linear  scale  and  for  different  values  of  removed  thundercloud  charge.  It 
was  assumed  in  all  cases  that  the  charge  w'as  removed  by  a  lightning  discharge  from  10  km 
altitude  in  1  ms. 


The  highly  nonlinear  dependence  of  the  optical  emission  intensity  on  Q  as  demonstrated 
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in  Figure  3.8  is  also  supported  by  the  altitude  profiles  of  resultant  optical  emission  intensity 
at  the  termination  of  the  lightning  discharge  ( t  =  1  ms)  as  shown  in  Figure  3.9a  for  different 
values  of  thundercloud  charges. 


Fig.  3.9.  Instantaneous  optical  emission  intensity,  (a)  Allilude  scans  at  r=0  of  optical 
emission  intensity  corresponding  to  1st  positive  band  of  No  for  different  values  of  removed 
thundercloud  charge;  (b)  Altitude  scans  at  r=0  of  optical  emission  intensity  corresponding 
to  different  optical  bands  for  Q=200  C. 


The  altitude  distribution  of  the  optical  intensities  corresponding  to  the  first  and  second 
positive  bands  of  N2,  Meinel  and  first  negative  bands  of  Nj,  and  first  negative  band  of  Ot 
are  shown  in  Figure  3.9b.  The  red  band  (first  positive  band  of  N2)  is  clearly  the  dominant 
optical  emission,  with  its  intensity  being  ~10  times  larger  than  that  of  the  blue  (N2  second 
positive  band),  consistent  with  the  dominant  red  color  of  the  sprites  [Pasko  et  al.,  199?]. 
The  Meinel  band  intensity  is  approximately  three  orders  of  magnitude  below  than  that  of 
the  N;  first  positive  band.  The  first  negative  bands  of  Nj  and  0|  are  approximately  four 
orders  of  magnitude  lower  than  the  intensity  of  the  dominant  N2  1st  positive  band.  These 
definitive  predictions  concerning  relative  excitation  levels  of  different  bands  provide  a  basis 
for  direct  comparisons  with  spectroscopic  measurements. 
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3.6  ELECTRON  AND  ION  CONDUCTIVITY 

Due  to  the  highly  nonlinear  nature  of  the  interaction  of  quasi-electrostatic  thundercloud 
fields  with  the  mesosphere  and  the  lower  ionosphere  it  is  expected  that  both  the  electron 
density  changes  and  the  optical  emission  intensities  are  sensitively  dependent  on  the  ambient 
electron  and  ion  conductivity  profiles. 

Figure  3. 1 0  compares  the  electron  density  changes  and  optical  emission  intensities  corre¬ 
sponding  to  the  three  models  of  ambient  electron  number  density  shown  in  Figure  2.4a.  All 
of  the  results  shown  up  to  this  point  have  been  for  the  most  tenuous  ambient  electron  den¬ 
sity  (profile  1 ).  Two  major  effects  associated  with  larger  values  of  ambient  electron  density 
(profiles  2,3).  Firstly,  the  relative  ionization  changes  are  lower,  with  the  maximum  changes 
occurring  at  lower  altitudes.  Secondly,  the  corresponding  optical  emission  intensities  are 
dramatically  reduced  at  the  higher  altitudes.  Both  effects  are  due  to  the  shorter  relaxation 
time  of  the  electric  field  at  high  altitudes  for  higher  ambient  electron  densities.  The  peak 
intensity  of  optical  emissions  at  ~  70  km  is  not  affected  by  the  ambient  electron  number 
density  since  this  quantity  is  defined  mostly  by  the  amplitude  of  the  electric  field  (in  other 
words  by  the  value  and  amplitude  of  the  removed  thundercloud  charge). 

The  dependence  of  ionization  and  optical  effects  on  the  ambient  ion  conductivity  is 
illustrated  in  Figure  3.1 1.  For  the  same  value  of  the  removed  thundercloud  charge  (100  C), 
but  for  different  ion  conductivity  profiles,  we  find  either  depletion  in  electron  density  due  to 
attachment  or  increased  density  due  to  ionization.  In  terms  of  optical  emissions,  the  results 
obtained  with  profiles  B,  C  and  A  differ  by  more  than  two  orders  of  magnitude  at  ~  80  km. 
We  note  that  the  ion  conductivity  profiles  A  and  C  are  similar  to  ion  conductivity  profiles 
experimentally  measured  and/or  used  in  previous  studies  [e.g.,  Dejnakarintra  and  Park , 
1974;  Holzworth  et  al.,  1985].  Our  results  for  these  two  profiles  can  thus  be  considered 
respectively  as  lower  and  upper  bound  estimates  of  expected  levels  of  optical  emissions  and 
ionization  changes. 


Altitude  (km)  Altitude  (km) 
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First  positive  band  of  N2 


I  (R) 


Fig.  3.10.  Dependence  on  the  ambient  electron  density,  (a)  Results  of  calculations  for 
electron  density  changes  corresponding  to  different  models  of  ambient  electron  density 
(Figure  2.4a).  The  ambient  profiles  of  electron  density  are  shown  for  reference  by  dashed 
line;  (b)  Dependence  of  the  optical  emission  intensity  on  the  model  of  ambient  electron 
number  density  (Figure  2.4a).  It  was  assumed  in  all  cases  that  a  charge  of  200  C  was 
removed  by  the  lightning  discharge  from  10  km  altitude  in  1  ms. 


Fig.  3.11.  Dependence  on  the  different  ambient  ion  conductivity.  Dependence  of  the 
optical  emission  intensity  (right  panel)  on  ambient  ion  conductivity  (Figure  2.4b).  Left  panel 
shows  results  of  calculations  for  electron  number  density  changes.  The  assumed  ambient 
profile  of  electron  density  is  shown  for  reference  by  dashed  line.  It  was  assumed  in  all  cases 
that  a  charge  of  1 00  C  was  removed  by  the  lightning  discharge  from  10  km  altitude  in  1  ms. 


3.7  LIGHTNING  DISCHARGE  DURATION 


The  production  of  intense  optical  emissions  by  quasi-electrostatic  fields  does  not  necessarily 
require  large  charge  transfer  rates  from  cloud  to  ground.  In  this  context,  it  is  important  to 
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First  positive  band  of  N2 


Fig.  3.12.  Model  results  for  different  discharge  durations  and  ambient  conductivity 
profiles.  Dependence  of  the  optical  emission  intensity  (right  panel)  on  ambient  conductivity 
profile  (left  panel)  and  discharge  duration  ( 1  ms  and  5  ms).  It  was  assumed  that  a  charge  of 
200  C  was  removed  by  the  lightning  discharge  from  10  km  altitude. 


note  that  in  all  of  the  results  shown  up  to  now  the  duration  of  the  lightning  discharge  was 
assumed  to  be  1  ms.  The  rate  of  charge  transfer  required  for  achieving  observable  levels 
of  optical  emissions  depends,  in  addition  to  the  magnitude  of  the  total  removed  charge, 
upon  the  conductivity  profile  in  the  60-90  km  altitude  range,  and  the  corresponding  relation 
between  the  duration  of  the  discharge  and  the  specific  relaxation  time  of  the  electric  field  in 
the  conducting  medium.  If  we  assume  that  the  value  of  the  total  removed  charge  remains 
constant  the  observable  levels  of  optical  emissions  considered  in  previous  sections  can  still 
be  achieved  with  longer  charge  transfer  times  as  long  as  they  do  not  significantly  exceed  the 
relaxation  time  at  50-90  km  altitudes. 

To  illustrate  this  point,  we  show  in  Figure  3.12  results  for  two  different  ambient  ion 
conductivity  profiles  which  differ  by  a  factor  of  five  at  <60  km  altitudes  where  the  ion 
conductivity  is  dominant.  Profile  A/5  shown  in  Figure  3.12  represents  conductivity  values 
reduced  by  a  factor  of  5  with  respect  to  profile  A.  We  assume  the  electron  component  of 
the  ambient  conductivity  to  be  the  same  for  both  profiles,  and  determined  by  the  electron 
density  profile  1  shown  in  Figure  2.4a. 

Optical  emission  intensities  corresponding  to  two  different  durations  of  the  lightning 
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discharge  ( 1  ms  and  5  ms)  for  profile  A,  and  for  5  ms  duration  for  profile  A/5  are  shown  in 
the  right  hand  panel  of  Figure  3.12.  Results  shown  correspond  to  times  at  the  termination 
of  the  lightning  discharges  (i.e.,  1  ms  or  5  ms).  The  intensity  of  the  optical  emission  is 
significantly  reduced  for  profile  A  and  the  5  ms  discharge  duration  in  comparison  with  the 
1  ms  case.  However,  the  intensity  distribution  is  similar  for  profile  A  with  1  ms  discharge 
duration  and  profile  A/5  with  the  5  ms  duration.  The  fast  drop  of  the  optical  emission 
intensity  above  85  km  altitude  for  the  latter  case  is  due  to  the  increase  of  conductivity  of 
the  medium  resulting  from  ionization.  Since  the  range  of  known  ion  conductivity  profiles 
span  more  than  an  order  of  magnitude  at  mesospheric  altitudes  (Figure  2.4b),  one  can 
conclude  that  similar  optical  emission  intensities  would  occur  with  charge  removal  times 
substantially  longer  than  the  1  ms  used  in  most  of  the  calculations  reported  in  previous 
sections.  For  lower  ion  conductivities,  the  predicted  optical  output  would  be  different 
primarily  in  terms  of  the  duration  of  the  emission  at  ionospheric  altitudes.  The  peak  optical 
emission  intensities  do  not  depend  significantly  on  the  ambient  conductivity  profile  since  this 
quantity  is  determined  primarily  by  the  magnitude  of  the  electric  field  (and  thus  the  value  and 
the  altitude  of  the  removed  thundercloud  charge).  The  latter  point  was  already  mentioned 
as  part  of  the  discussion  of  results  corresponding  to  different  ambient  electron  density 
profiles  (Figure  3. 1 0),  and  is  also  illustrated  by  the  results  shown  in  Figure  3. 1 1  indicating 
a  saturation  effect  corresponding  to  profiles  B  and  C.  Further  reduction  of  conductivity  at 
mesospheric  altitudes  (below  those  represented  by  profile  C)  would  not  produce  a  significant 
increase  in  the  optical  emission  intensity  or  the  ionization. 

Although  up  to  now  we  have  primarily  considered  the  removal  of  <300  C  of  charge  in 
~  1  ms,  removal  of  larger  amounts  of  charge  (e.g.,  ~400  C)  over  longer  times  (e.g.,  ~  1 0  ms) 
may  lead  to  the  formation  of  streamer-type  channels  of  breakdown  ionization  at  altitudes  ~ 
50-70  km.  This  regime  of  excitation  and  various  consequences  are  discussed  in  Chapter  4. 
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3.8  LIGHTNING  DISCHARGE  TYPE 


Figure  3. 1 3  shows  postdischarge  distributions  of  electric  field  and  optical  emission  intensity 
as  a  function  of  altitude  for  three  different  lightning  discharge  models,  M 1 ,  M2  and  +CG  as 
shown  in  Figure  2. 1 .  We  note  that  the  M2  case  corresponds  to  the  deposition  of  a  monopole 
charge,  in  contrast  to  removal  of  charge  as  represented  by  the  Ml  and  +CG  cases. 


optical  emission  intensity  (right  panel)  on  the  model  oflightning  discharge  (Ml,  M2,  +CG) 
(Figure  2.1).  Left  panel  shows  the  corresponding  distribution  of  postdischarge  electric  field 
at  r  =  0.  All  calculations  were  performed  for  the  removed  charge  value  200  C. 


The  results  obtained  for  three  different  lightning  discharge  models  are  identical  above 
~  60  km  altitude,  as  expected  since  the  electric  field  at  high  altitudes  is  defined  mainly 
by  the  value  of  the  removed/deposited  charge  and  its  altitude,  independent  of  the  charge 
configuration  inside  the  cloud  before  or  after  the  lightning  discharge.  The  range  of  altitudes 
over  which  the  electric  field  is  independent  of  the  lightning  discharge  model  extends  lower 
with  increasing  charge  accumulation  time.  In  this  Chapter  we  use  a  charge  accumulation 
time  of  r,=0.5  sec  [Pasko  et  ai,  1 995]  since  optical  emissions  appear  at  altitudes  well  above 
60  km.  Calculations  with  rs  longer  then  0.5  sec  lead  to  identical  results,  and  rs=  0.5  sec 
value  was  chosen  from  the  point  of  view  of  computational  efficiency. 
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Figure  3.14  shows  results  for  +CG,  -CG  and  IC  discharges.  The  +CG  flashes  clearly 
produce  larger  electric  fields  and  optical  emission  intensities,  with  the  optical  intensity  being 
three  orders  of  magnitude  above  that  for  -CG  and  IC.  This  result  is  a  simple  consequence  of 
the  fact  that  +CG  discharges  remove  charge  from  the  higher  altitude,  producing  larger  values 
of  postdischarge  electric  fields.  Note  that  for  the  IC  case  the  field  of  the  lower  removed 
charge  significantly  reduces  the  field  of  the  upper  charge  and  the  resulting  optical  output.  It 
should  be  noted  that  we  model  a  +CG  discharge  as  the  removal  of  the  thundercloud  positive 
charge  which  is  assumed  to  be  located  at  a  higher  altitudes  than  the  accompanying  negative 
charge.  In  a  real  situation,  positive  charge  may  well  be  accumulated  at  and  removed  from 
the  lower  parts  of  the  thundercloud  [e.g.,  Boccippio  et  al.,  1995;  Marshall  et  al.,  1996]. 


Fig.  3.14.  Dependence  on  the  type  of  lightning  discharge.  Results  are  shown  in  the  same 
format  as  in  Figure  3.13  but  for  +CG,  -CG,  and  IC  (Figure  2.1)  models  of  the  lightning 
discharge. 
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3.9  THUNDERCLOUD  CHARGE  GEOMETRY 

The  dependence  of  optical  emission  characteristics  on  the  geometry  of  the  thundercloud 
charge  is  of  particular  interest  since  sprites  usually  appear  above  mesoscale  convective 
systems  with  spatially  distributed  charges  [E.  Williams,  private  communication;  Boccippio 
et  al. ,  1995;  Marshall  et  al.,  1996].  Figure  3.15  shows  the  postdischarge  distribution  of 
electric  fields  and  optical  emissions  corresponding  to  +CG  discharge  for  three  different 
geometries  of  the  thundercloud,  namely,  ( 1 )  Gaussian  spherically  symmetric  charge  density 
distributions  with  scale  a=3  km  (used  in  most  of  the  calculations  in  this  thesis);  (2)  Also 
Gaussian  spherically  symmetric  distribution  but  with  a=6  km;  (3)  Gaussian  distribution  in 
the  z  direction  with  a=3  km  and  flat  (uniform)  in  the  horizontal  direction  with  a  “disk"  of 
radius  30  km. 


Fig.  3.15.  Dependence  on  the  thundercloud  charge  geometries.  Results  of  calculations 
of  postdischarge  electric  fields  (left  panel)  and  optical  emission  intensities  of  the  first  positive 
band  of  N2  (right  panel)  for  thundercloud  charge  distributions  of:  (i)  Gaussian  spherically 
symmetric  charge  density  distributions  with  scale  a=3  km  (used  in  most  of  the  calculations 
in  this  thesis);  (ii)  Also  Gaussian  spherically  symmetric  distribution  but  with  a=6  km;  (iii) 
Gaussian  distribution  in  the  2  direction  with  a=3  km  and  flat  (uniform)  in  the  horizontal 
direction  with  a  “disc”  of  radius  30  km 


The  results  in  Figure  3.15  show  that  horizontally  spread  charges  generally  lead  to  a 
lower  value  (for  the  same  value  Q  of  the  removed  thundercloud  charge)  of  the  electric  fields 
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everywhere  above  the  thundercloud.  The  size  of  the  spherical  thundercloud  charges  (3  and 
6  km  scales)  does  not  appear  to  significantly  affect  the  electric  fields  and  optical  outputs 
at  mesospheric  and  ionospheric  altitudes.  In  this  connection  it  is  interesting  to  note  that 
if  there  were  no  upper  conducting  boundary  (ionosphere),  increasing  the  effective  radius 
of  the  disk  of  charge  would  tend  to  approach  the  one  dimensional  geometry  of  the  electric 
field  due  to  a  sheet  of  charge.  Since  the  thundercloud  charge  always  has  an  image  charge 
below  the  ground,  the  field  would  be  enhanced  between  the  disk  and  ground  but  would 
rapidly  decrease  above  it  with  increasing  disk  radius.  Since  the  ionospheric  boundary  also 
generates  an  image  charge,  the  amplitude  of  the  field  is  relatively  enhanced.  Results  of 
calculations  shown  in  Figure  3.15  suggest,  however,  that  the  spherical  geometry  gives  the 
largest  possible  amplitudes  of  electric  fields  above  the  thundercloud.  On  the  other  hand,  the 
disk  projects  a  significant  electrc  field  over  a  wider  (in  the  radial  direction)  region  of  space. 

The  electric  field  above  the  startiform  clouds  associated  with  mesoscale  convective  sys¬ 
tems  have  been  measured  by  Marshall  et  al.  [1996]  using  balloon-borne  instruments.  The 
measured  amplitudes  of  electric  fields  were  explained  by  considering  the  disk  charge  with 
charge  density  ~3  nC  m~3  centered  at  4  km  altitude,  with  thickness  400  m  (effective  surface 
charge  density  ~1200  nC  m“2)  and  diameters  ranging  from  20  to  200  km.  The  authors 
accounted  for  image  charge  in  the  ground,  but  ignored  the  effects  of  ionospheric  images. 
The  effects  of  the  upper  ionospheric  boundary  must  be  taken  into  account  for  a  correct 
determination  of  the  electric  field  values  at  ionospheric  altitudes  and  the  total  charge  ac¬ 
cumulated  in  the  thundercloud,  especially  for  horizontally  spread  charge  distributions.  We 
note  that  even  relatively  small  values  of  charge  (200  C)  distributed  over  a  disk  with  diameter 
60  km  (surface  charge  density  only  ~  100  nC/m2)  would  lead  to  ~1  kV/m  electric  field 
change  above  the  thundercloud  (Figure  3.15),  in  excellent  agreement  with  observations  of 
Marshall  et  al.  [  1 996].  Thus  relatively  large  values  ~  1 200  nC/m2  suggested  by  Marshall  et 
al.  [  1 996]  and  derived  without  accounting  for  the  effects  of  the  upper  ionospheric  boundary 
may  not  in  fact  be  required  for  production  of  observed  amplitudes  of  the  electric  field. 
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3.10  RELAXATION  TIMES  OF  QUASI-ELECTROSTATIC  FIELDS 

The  dynamics  of  the  electric  field  in  the  conducting  medium  above  thunderclouds  directly 
control  the  ionization  changes  and  the  duration  of  optical  emissions  associated  with  sprites. 
An  important  aspect  of  our  model  with  respect  to  previous  work  on  the  relaxation  of  quasi¬ 
electrostatic  fields  [ Dejnakarintra  and  Park,  1974;  Greifinger  and  Greifinger,  1976;  Hale 
and  Baginski,  1987;  Bagitiski  et  al„  1988;  Hale,  1994]  is  the  self-consistent  modeling  of 
the  conductivity  a,  as  discussed  in  Section  2.4.  To  demonstrate  the  importance  of  the 
self-consistent  treatment  of  conductivity  in  comparison  with  previous  work,  we  compare 
in  Figure  3.16  the  temporal  variation  of  the  electric  field  \EZ  \  at  70  and  80  km  for  three 
different  treatments  of  conductivity  a  for  the  case  of  removal  of  a  monopole  charge  (case 
M 1  shown  in  Figure  2. 1 ). 


Fig.  3.16.  Relaxation  of  the  quasi-electrostatic  field.  The  vertical  component  of  the 
electric  field  at  two  selected  altitudes  at  r  =  0  as  a  function  of  time  since  the  beginning  of 
the  lightning  stroke  for  a  200  C  source  charge  for  different  solutions  of  conductivity  cr.  (a) 
self-consistent  formulation;  (b)  conductivity  kept  constant  in  time;  (c)  relaxation  with  the 
local  ambient  relaxation  time  (e0/o)\  (d)  same  as  (b)  but  for  charge  deposited  during  1  ms 
(M2,  Figure  2a)  (rather  than  removal)  at  the  beginning  of  the  simulation. 


At  2  =80  km,  the  electric  field  \EZ\  for  case  (a)  is  ~4  times  that  for  (b),  and  initially 
exhibits  ~3  times  longer  time  duration,  followed  by  a  faster  decrease,  resulting  from  the 
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initial  decrease  of  conductivity  a  due  to  heating  followed  by  its  increase  due  to  ionization 
at  later  times.  For  case  (d)  at  both  altitudes  (70  and  80  km)  the  relaxation  rate  of  the  electric 
field  slows  down  even  in  comparison  with  cases  (a)  and  (b);  this  effect  is  due  to  the  formation 
of  the  quiescent  electric  field  around  the  source  charge.  This  steady  field  can  be  found  as 
the  stationary  solution  of  the  problem  of  a  point  charge  in  a  medium  with  exponentially 
increasing  (in  altitude)  conductivity  [e.g.,  Holzer  and  Saxon,  1952],  In  all  cases,  we  see 
significant  slowing  down  of  the  relaxation  rate  of  the  electric  field  in  comparison  with  the 
local  ambient  e0ja  (marked  (c)  in  Figure  3.16)  rate.  The  most  recent  discussion  of  this 
effect,  which  has  been  observed  experimentally  and  extensively  studied  theoretically,  is 
given  by  Hale  [1994]  (see  also  references  cited  therein). 

Since  the  slowing  down  of  the  electric  field  relaxation  rate  is  observed  on  a  millisecond 
time  scale  [e.g.,  Dejnakarintra  and  Park,  1974;  Baginski  et  al.,  1988;  Hale,  1994]  our 
quasi-electrostatic  formulation  which  is  valid  for  the  near  zone  fields  is  fully  sufficient  for 
the  investigation  of  the  temporal  evolution  of  the  quasi-static  fields  and  its  consequences  in 
the  form  of  heating,  ionization  and  optical  emissions.  That  this  is  so  can  also  be  confirmed 
by  recent  results  of  Baginski  and  Hodel  [1994]  who  provide  a  detailed  comparison  of  the 
solution  of  the  system  of  equations  (2. 1 )  and  (2.2)  with  the  full  system  of  Maxwell  equations. 

The  slowing  down  effect  can  be  physically  understood  as  follows.  The  atmospheric 
conductivity  exponentially  increases  with  altitude.  There  exists  an  altitude  below  which 
we  can  consider  charges  (both  thundercloud  and  induced  polarization  charges  in  the  con¬ 
ducting  atmosphere)  to  be  “stationary”  on  a  certain  time  scale.  The  slowing  down  of  the 
electric  field  relaxation  can  be  interpreted  as  the  establishment  of  a  quasi-stationary  solu¬ 
tion  above  this  “stationary”  charge.  As  time  progresses,  the  imaginary  boundary  below 
which  we  can  consider  charges  and  associated  fields  to  be  stationary  moves  dow  nward  [e.g.. 
G r elfin ger  and  Greifinger,  1976];  however,  the  rate  of  the  field  and  charge  relaxation  above 
this  boundary  always  remains  slower  than  the  local  relaxation  rate  £0/a.  At  each  altitude 
above  this  boundary  there  is  a  nonvanishing  value  of  induced  charge  and  current  density 
which  is  a  consequence  of  the  inhomogeneity  of  the  atmospheric  conductivity  [ Greifinger 
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and  Greifinger,  1976].  The  relaxation  rate  would  be  even  further  slowed  down  if  there  is  a 
source  charge  which  is  maintained  by  an  external  effect  (e.g.,  regeneration  of  charge  within 
the  thundercloud  due  to  updraft).  Such  a  circumstance  is  represented  by  the  deposition  of 
charge  over  1  ms,  as  in  the  case  marked  (d)  in  Figure  3.16. 

The  imaginary  downward  moving  boundary  mentioned  above  is  essentially  the  “moving 
capacitor  plate”  model  introduced  first  by  Greifinger  and  Greifinger  [197 6]  and  successfully 
used  by  Hale  and  Baginski  [1987]  to  estimate  the  current  which  flows  to  the  ionosphere 
following  a  lightning  stroke.  Greifinger  and  Greifinger  [1976]  define  this  boundary  as 
separating  two  regions  dominated  by  conduction  (above)  and  displacement  (below)  currents. 

Any  persisting  charges  at  low  altitudes  constitute  a  source  of  electric  field  at  high  altitudes, 
so  that  the  electric  field  at  high  altitudes  completely  relaxes  to  zero  only  when  the  last  re¬ 
maining  charges  dissipate  at  low  altitudes,  a  process  which  occurs  with  long  relaxation  times 
due  to  low  atmospheric  conductivity.  Until  such  time  that  the  polarization  charges  at  the 
lowest  altitudes  dissipate  away,  the  electric  fields  at  high  altitudes  continuously  goes  through 
a  series  of  quasi-stationary  states  corresponding  to  persisting  at  low  altitudes  charges.  This 
effect  leads  to  the  apparent  slowing  down  of  the  relaxation  of  the  electric  fields. 

An  important  role  is  played  by  surface  charges  induced  at  the  highly  conducting  ground 
and  lower  ionospheric  boundaries  which  are  assumed  to  be  maintained  at  a  constant  po¬ 
tential.  These  charges  change  dynamically  in  time  consistent  with  the  electric  fields  and 
correspondingly  affect  the  dynamics  of  the  fields.  The  slowing  down  effect  is  not  a  property 
of  two-dimensional  system  and  is  also  observed  in  one  dimensional  geometry. 

To  illustrate  the  above  points,  we  show  in  Figure  3.17  results  obtained  for  cases  similar 
to  those  shown  in  Figure  3.16b,d,  but  for  the  deposition  and  removal  of  a  “disk”  of  charge 
of  10  C  with  radius  150  km  and  vertical  extent  3  km.  The  radial  size  of  the  system  was 
assumed  to  be  also  150  km  by  adopting  the  boundary  condition  d<p/dr  =  0  at  r  =  150 
km.  The  electric  field  in  this  case  is  essentially  one-dimensional  and  is  directed  vertically 
everywhere  in  the  system.  Similar  to  those  of  Figure  3.16  ,  results  of  Figure  3.17  show 
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slowing  down  of  the  electric  field  relaxation.  Note  that  the  field  shown  is  that  at  50  km 
altitude  so  that  the  relaxation  of  the  field  in  this  case  is  entirely  determined  by  the  gradients 
of  conductivity  and  persisting  charges  (no  heating  effects  included). 


Fig.  3.17.  One  dimensional  relaxation  of  the  electric  field.  The  electric  field  at  altitude 
50  km  as  a  function  of  time:  (i)  relaxation  w  ith  the  local  ambient  relaxation  lime  (e0/a)\ 
(ii)  deposition  in  10  ms  of  a  disk  of  charge  of  10C  with  diameter  300  km  and  vertical  extent 
3  km  at  altitude  20  km:  (iii)  the  same  as  (ii)  but  the  charge  w'as  accumulated  during  1  sec 
and  then  removed  in  10  ms. 


In  summary,  we  can  identify  three  separate  physical  effects  which  contribute  to  the  slow¬ 
ing  down  of  the  relaxation  rate  of  the  quasi-static  electric  field:  (i)  the  decrease  of  the  ambient 
electron  component  of  the  conductivity  due  to  selfconsistent  heating  by  quasi-electrostatic 
fields;  (ii)  the  generation  of  additional  surface  charges  at  the  highly  conducting  ground 
and  the  lower  ionospheric  boundaries;  (iii)  the  establishment  of  quasi-stationary  distribu¬ 
tions  of  the  electric  fields  and  charge  densities  around  persisting  thundercloud  charges  and 
polarization  charges  in  the  conducting  atmosphere  above  the  thundercloud. 


4 

Sprites  as  Luminous  Columns  of 
Ionization 


In  this  chapter  we  investigate  the  formation  of  columnar  channels  of  breakdown  ionization 
and  carrot-like  vertical  luminous  structures  depending  on  the  dynamics  of  the  thundercloud 
charge  removal  by  continuing  currents  associated  with  positive  cloud-to-ground  lightning 
flashes.  We  show  that  the  quasi-electrostatic  thundercloud  fields  can  potentially  lead  to 
intense  heating  and  ionization  in  structures  with  typical  transverse  dimension  ~5-10  km 
and  spanning  an  altitude  range  from  ~80km  to  well  below  ~50  km.  These  carrot-like  forms 
closely  resemble  those  observed  in  sprites.  We  also  show  that  the  appearance  of  optical 
emissions  can  be  significantly  delayed  in  time  (~  1  -20  ms)  with  respect  to  the  causative 
lightning  discharge.  We  compare  results  of  model  calculations  with  recent  video,  optical, 
and  photometric  observations  of  sprites. 

4.1  PHYSICAL  MECHANISM 

In  the  previous  chapter  it  was  illustrated  that  large  electric  field  transients  capable  of  causing 
breakdown  ionization  at  mesospheric  altitudes  can  he  produced  in  the  region  between  the 
thundercloud  tops  and  the  lower  ionosphere  following  the  removal  of  thundercloud  charge  by 
a  cloud-to-ground  discharge  (see  Figure  3. 1  and  accompanying  discussion  in  Section  3. 1 ). 

During  the  process  of  charge  removal  from  the  thundercloud  (typically  1-20  ms),  the 
intensities  of  the  quasi-electrostatic  fields  slowly  increase,  and  may  exceed  the  threshold 
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of  breakdown  ionization  in  a  localized  region.  The  altitude  at  which  this  would  occur, 
and  the  radial  extent  of  the  ionization  region,  depends  on  the  value  of  the  threshold  field, 
which  in  turn  depends  on  the  atmospheric  neutral  density.  Once  the  ionization  threshold 
is  exceeded,  the  conductivity  in  this  localized  region  rapidly  increases  by  typically  several 
orders  of  magnitude  as  was  shown  in  the  previous  chapter.  The  enhanced  conductivity  in 
such  a  localized  region  in  turn  leads  to  faster  rates  of  electric  field  relaxation  within  it  and 
to  significantly  enhanced  electric  field  in  its  vicinity,  much  like  the  case  of  highly  curved 
surfaces  (sharp  edges)  of  conducting  bodies  placed  in  an  electric  field.  If  the  relaxation  rate 
of  the  surrounding  medium  is  not  high  enough  to  compensate  for  the  enhanced  electric  field, 
the  field  can  further  exceed  the  local  threshold  for  breakdown  ionization  and  spatially  extend 
the  region  of  high  conductivity.  If  the  curvature  of  the  resultant  conductivity  inhomogeneity 
is  still  large  enough,  continuous  development  and  propagation  of  the  ionization  channel  may 
thus  occur.  Formation  of  these  ionization  channels  may  be  closely  associated  with  streamer 
type  processes  characterized  as  filamentary  plasmas,  the  dynamics  of  which  are  controlled 
by  highly  localized  nonlinear  space  charge  waves  [e.g.,  Uman ,  1984,  p.  202;  Vitello  et  ai, 
1994],  In  general,  we  might  expect  the  curvature  of  the  inhomogeneity  to  sharpen  (i.e.,  the 
effective  radius  of  the  ionized  column  to  decrease)  as  the  streamer  propagates  downward. 
However,  as  we  show  later,  depending  on  the  temporal  profile  of  charge  removal  from  the 
thundercloud  ionization  columns  with  a  variety  of  shapes  can  be  created,  including  radial 
convergence  or  divergence  of  the  luminous  regions  as  the  streamer  moves  downward. 

On  the  basis  of  these  arguments,  the  regions  of  air  breakdown  may  have  the  shape  of  an 
ionization  column(s)  which  electrically  connects  the  highly  conducting  lower  ionosphere  to 
lower  atmospheric  regions.  An  ionization  channel  may  develop  from  a  “seed”  enhancement 
of  local  conductivity  of  any  origin  (see  next  sections)  and  its  propagation  is  physically  con¬ 
trolled  by  the  relaxation  of  the  electric  field  inside  its  conducting  body  and  the  enhancement 
of  the  field  in  the  vicinity  of  its  tip  as  described  above. 
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4.2  MODEL  PARAMETERS 

For  the  calculations  in  this  chapter  we  assume  the  ambient  electron  density  profile  1  shown 
in  Figure  2.4a  and  ambient  ion  conductivity  profile  A  illustrated  in  Figure  2.4b.  We  use 
the  quasi-electrostatic  model  described  in  Chapters  2  and  3  for  five  different  models  of 
charge  removal  shown  in  Figure  4. 1  with  time  f= 0  corresponding  to  the  start  of  the  charge 
removal  process.  These  lightning  discharge  models  represent  the  removal  of  larger  amounts 
of  total  charge  while  being  consistent  with  reasonable  lightning  current  values  of  <300  kA. 
As  a  result,  generally  larger  peak  magnitudes  of  electric  fields  are  produced,  leading  to  a 
different  class  of  effects  than  those  discussed  in  Chapter  3. 

We  note  that  our  model  is  effectively  a  one  fluid  model,  and  in  contrast  with  the  two-fluid 
approaches  [e.g.,  Vitello  et  al.,  1 994]  it  is  based  on  the  continuity  equation  for  charge  density 
(2.3)  and  may  not  accurately  predict  the  small  scale  details  of  the  electron  density  near  the 
steep  boundaries  associated  with  the  front  of  streamers.  Calculations  indicate  that  downward 
propagating  ionization  channels  tend  to  focus  to  very  sharp  edges,  an  effect  which  often 
leads  to  numerical  instability  at  low  altitudes  (~50  km).  For  modeling  of  the  downward 
propagation  we  somewhat  arbitrarily  introduce  in  our  model  an  “effective”  recombination 
constant  a=l  cm  Vs  at  50  km  which  varies  exponentially  in  altitude  with  a  scale  height  of 
6  km,  accounting  for  the  increasing  electron  losses  at  low  altitudes,  and  also  serves  as  a 
smoothing  factor  preventing  unphysically  large  oscillations  in  the  electron  number  density. 
For  modeling  of  upward  propagating  streamers  which  can  be  launched  from  seeds  of  plasma 
at  low  altitudes  (see  next  section)  we  assume  that  the  dielectric  relaxation  time  (£0/<t)  behind 
the  front  of  the  channel  remains  constant  in  accordance  with  the  results  of  extensive  streamer 
simulations  [e.g.,  Vitello  et  al.,  1994].  This  assumption  also  is  used  in  the  next  chapter  for 
modeling  of  blue  jets.  The  thundercloud  charges  +Q ,  -Q  ( +Q  placed  at  10  km)  form  a 
vertical  dipole  which  develops  over  several  tens  of  seconds.  The  +Q  charge  is  then  rapidly 
discharged  to  the  ground.  The  response  of  the  system  at  mesospheric  altitudes  depends  on 
the  altitude  and  magnitude  of  +Q  and  does  not  depend  on  the  initial  charge  configuration 
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inside  the  cloud  (i.e.,  magnitude  and  altitude  of  - Q )  as  was  discussed  in  Sections  3.1  and 
3.3. 

The  top  panels  of  Figure  4.1  show  five  different  models  of  charge  removal  used  for 
calculations  in  this  chapter.  The  time 


Fig.  4.1.  Models  of  charge  removal.  Five  models  of  ihundercloud  charge  removal  used 
in  calculations:  (a)  removal  of  200  C  of  charge  in  1  ms;  (b)  300  C  in  5  ms;  (c)  400  C  in  10 
ms;  (d)  475  C  in  20  ms;  (c)  fast  initial  removal  of  200  C  in  1  ms  and  then  slow  removal  of 
275  C  in  19  ms. 


4.3  ELECTRON  DENSITY 

Results  of  calculations  of  the  electron  number  density  as  a  function  of  time  corresponding  to 
case  e  in  Figure  4. 1  are  shown  in  Figure  4.2a.  The  case  e  is  chosen  for  detailed  discussion 
since  it  most  closely  represents  the  typical  time  dynamics  of  the  charge  removal  from  a 
thundercloud  consisting  of  fast  initial  removal  of  a  portion  of  the  charge  by  a  lightning 
stroke  followed  by  a  more  prolonged  period  of  continuing  current  [e.g.,  Liman ,  1987,  p. 
200].  At  time  t= 2  ms,  a  region  of  significant  ionization  with  characteristic  radius  >  10  km 
is  produced  at  70-80  km  altitude  by  the  fast  (1  ms)  removal  of  the  first  200  C,  in  a  manner 
similar  to  the  cases  studied  earlier  in  Section  3.4.  Although  the  electric  field  is  enhanced 
just  outside  of  this  localized  region  of  enhanced  conductivity,  this  increase  is  compensated 
by  the  relatively  high  ambient  conductivity  at  these  altitudes,  inhibiting  the  formation  of  an 


SPRITES  AS  LUMINOUS  COLUMNS  OF  IONIZATION 


75 


ionization  channel.  If  the  removal  of  charge  in  this  case  were  limited  to  200  C,  no  ionization 
channel  would  be  formed  and  the  totality  of  the  resultant  effects  would  be  those  presented 
in  Chapter  3.  In  most  cases,  however,  the  fast  impulsive  removal  of  charge  is  followed  by  a 
period  of  continuing  current  which  usually  lasts  tens  of  ms  and  can  account  for  a  significant 
fraction  of  the  total  charge  lowered  by  a  cloud  to  ground  discharge  [Uman,  1987,  pp.  196, 
200],  Such  low  rates  of  charge  removal  can  in  some  cases  further  increase  the  curvature  of  the 
highly  conducting  region,  leading  to  the  formation  of  streamer-type  channels  of  breakdown 
ionization  propagating  downward  as  shown  in  Figure  4.2a.  Focusing  (Figure  4.2a,  f  =10 
ms)  and  subsequent  expansion  (Figure  4.2a,  (=20  ms)  of  ionization  toward  lower  altitudes 
is  caused  by  the  fact  that  the  additional  slow  charge  removal  enhances  the  electric  field 
and  w  idens  the  region  of  breakdown  ionization.  This  behavior  occurs  at  altitudes  (<60km) 
where  the  ambient  conductivity  is  low  enough  so  that  it  does  not  prevent  the  enhancement 
of  the  electric  field  near  the  tip  of  the  channel  in  contrast  with  that  at  high  altitudes  (~  70 
km). 

4.4  ELECTRIC  FIELD 

The  spatial  distribution  of  the  enhancement  of  the  electric  field  above  its  value  at  (=  1 5  ms  is 
shown  in  Figure  4.2b  at  different  time  intervals.  The  purpose  of  this  kind  of  normalization 
is  to  illustrate  the  enhancement  of  the  electric  field  in  regions  near  the  tip  of  the  streamer. 
The  characteristic  breakdown  field  Ef.  [ Papadopoulos  et  al.,  1993]  is  also  shown  in  the  top 
panel  of  Figure  4.2b.  The  electric  field  remains  approximately  equal  to  £V  at  the  tip  of 
the  streamer.  The  physical  reason  for  this  behavior  is  the  fact  that  when  the  electric  field 
exceeds  the  characteristic  breakdown  field  E ^.,  ionization  increases  lead  to  enhanced  local 
conductivity  which  in  turn  leads  to  the  rapid  relaxation  of  the  electric  field.  As  a  result, 
the  electric  field  tends  to  remain  near  E Thus,  the  ionization  channel  proceeds  with  the 
electric  field  at  its  tip  staying  just  barely  above  the  threshold  field. 
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Fig.  4.2.  Formation  and  downward  propagation  of  an  ionization  channel.  A  crossec- 
tional  view  of  the  distribution  of  (a)  the  electron  number  density  Ne\  (b)  absolute  values  of 
the  electric  field  E\  and  (c)  optical  emission  intensity  I  of  the  1st  positive  band  of  molecular 
nitrogen  N2  at  selected  instants  of  time.  Top  panels  show  altitude  scans  at  r  =  0. 
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4.5  OPTICAL  EMISSIONS 

The  emission  intensity  corresponding  to  the  first  positive  band  of  N2  is  shown  in  Figure  4.2c, 
and  is  significant  only  in  regions  where  the  electric  field  is  above  the  threshold  for  optical 
excitation  of  this  emission  band  of  N2.  The  significantly  reduced  electric  field  inside  the 
ionized  regions  (Figure  4.2a)  results  in  insignificant  optical  output.  The  optical  emission 
region  is  rather  wide  (radial  extent  is  ~20  km)  at  t= 2  ms  at  ~70  km  altitude  but  subsequently 
becomes  narrower  (~5  km)  at  ~60  km  altitude  at  i=  1 0  ms,  and  expands  downward  at  later 
times.  Downward  propagation  slows  down  below  ~50  km  altitudes  (not  shown)  due  to  the 
exponential  increase  in  electron  loss  at  these  altitudes.  The  average  speed  of  downward 
propagation  of  the  column  in  this  case  is  ~10^  km/s. 

The  optical  emissions  shown  in  Figure  4.2c  endure  for  times  less  than  or  comparable  to 
the  duration  of  one  video  frame  (1/30  sec).  For  comparison  with  video  observations,  it  is 
thus  useful  to  average  over  ~30  ms.  Figure  4.3  shows  the  distribution  of  the  intensity  of 
the  first  positive  band  of  N2  averaged  over  32  ms  for  cases  a,b,c.,d,  and  e.  Case  a  is  similar 
to  the  cases  of  fast  charge  removal  studied  in  the  previous  chapter  of  this  dissertation.  The 
increase  of  the  electric  field  outside  the  conducting  region  in  this  case  is  compensated  by  the 
relatively  high  ambient  conductivity,  inhibiting  the  formation  of  an  ionization  channel.  For 
case  c  the  electric  field  reaches  the  ionization  threshold  Ek  only  at  the  very  end  of  the  charge 
removal  at  t  ~10  ms,  so  that  we  see  the  self  driven  development  of  an  ionization  channel 
propagating  downward  which  tends  to  focus  toward  lower  altitudes.  For  case  d ,  the  total 
amount  of  removed  charge  (475  C)  is  greater  than  for  case  c  (400  C)  so  that  the  electric  field 
exceeds  E 'f.  over  a  larger  region  and  the  corresponding  structure  of  the  optical  emissions 
appears  to  be  wider  and  does  not  show  significant  focusing  toward  lower  altitudes.  Case  e 
shows  structure  consisting  of  two  parts  separated  by  a  gap  at  ~  67  km  altitude.  This  feature 
is  consistent  with  the  typical  structure  of  sprites  as  a  high  altitude  “plume”  and  adjacent 
lower  altitude  “dendritic”  (forked)  forms  [e.g.,  Winckler  et  ai,  1996]. 


In  case  b  we  placed  a  “ball”  of  highly  conducting  plasma  with  radius  ~5  km  and  electron 


Altitude  (km) 
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Fig.  4.3.  Time  averaged  optical  emission  intensity.  The  distribution  of  time  averaged 
optical  emission  intensity  I  (Rayleighs)  of  first  positive  band  of  N2  is  shown  on  linear  scale 
and  for  models  of  charge  removal  as  illustrated  in  Figure  4.1.  Note  1/8  scale  modification 
for  case  a. 


density  Nc  ~  1 04cm-3  at  45  km  altitude.  This  somewhat  artificial  initial  condition  is  used  to 
demonstrate  how  preexisting  inhomogeneities  of  the  ambient  conductivity  may  dramatically 
modify  the  properties  of  the  breakdown.  The  ionization  channel  is  initiated  in  case  b  at  ~50 
km  altitude  due  to  the  enhanced  conductivity  and  the  associated  electric  field  and  propagates 
upward  showing  a  slight  expansion  toward  higher  altitudes.  The  propagation  speed  is  very 
high  (~  104  km/s)  and  is  mostly  defined  by  the  almost  simultaneous  enhancement  of  the 
electric  field  above  the  characteristic  breakdown  field  in  a  region  between  the  charge 
induced  in  the  tip  of  the  column  and  the  upper  conducting  boundary  assumed  to  be  at  80 
km. 
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4.6  CURRENT  AND  CHARGE  TRANSFER  TO  THE  IONOSPHERE 

The  current  which  flows  within  the  ionization  channel  is  dominantly  due  to  the  electron 
component  of  conductivity.  Form  Figure  4.2  the  radius  of  the  column  can  be  estimated  to 
be  ~5  km,  and  the  electron  density  is  ~102cm-3.  The  electron  drift  velocity  in  an  electric 
field  E  comparable  to  the  critical  breakdown  field  E ^  (Section  2.4.2)  can  be  estimated  as 
Vd  =  where  for  E  ~  E *  ~  4  x  10~2Aro/Ar  m2/V/s  (see  equation  (2.5)).  The 

resultant  value  of  the  electron  drift  velocity  is  i'd  ~  105  m/s,  so  that  the  current  flowing 
within  the  ionization  channel  is  ~100  A.  The  total  charge  transfer  to  the  ionosphere  by  this 
current  on  a  time  scale  of  ~30  ms  (roughly  the  duration  of  optical  emissions)  is  ~3  C.  The 
optical  emission,  however,  is  not  a  direct  indication  of  current  flow  in  the  ionization  channel 
and  most  of  the  charge  to  the  ionosphere  is  transferred  on  much  longer  time  scales  (hundreds 
of  ms)  well  after  the  initial  transient  appearance  of  the  optical  emission.  Since  the  electric 
field  which  appears  at  mesospheric  altitudes  follow  ing  a  cloud-to-ground  lightning  is  mostly 
defined  by  the  magnitude  and  altitude  of  the  removed  (monopole)  charge  and  its  image  in 
the  ground,  it  is  expected  that  the  total  transfer  of  negative  charge  (for  +CG  discharge  case) 
to  the  ionosphere  is  comparable  with  the  amount  of  positive  charge  originally  removed  from 
the  cloud.  Physically,  conduction  currents  flowing  in  the  atmosphere  and  mesosphere  after 
a  +CG  discharge  tend  to  bring  positive  charge  down  from  the  ionosphere  to  compensate  the 
postdischarge  electric  field.  The  amount  of  charge  required  to  achieve  this  compensation 
is  comparable  to  the  charge  removed  from  the  thundercloud.  It  is  illustrative  to  note  in 
this  context  that  a  static  charge  placed  in  an  inhomogeneous  conducting  medium  generates 
around  itself  an  integral  amount  of  screening  charge  (with  opposite  sign)  which  is  equal  in 
magnitude  to  the  originally  placed  charge.  Ionization  channels  associated  with  sprites  act 
as  highly  conducting  wires  connecting  the  ionosphere  and  lower  atmospheric  regions  which 
contribute  to  faster  compensation  of  the  postdicharge  electric  fields.  However,  it  should 
be  stressed  that  the  electric  field. can  be  relaxed  simply  due  to  the  ambient  atmospheric 
conductivity,  and  that  the  associated  current  flow  from/to  the  ionosphere  in  case  of  +CG/- 
CG  must  normally  accompany  of  every  lightning  discharge.  This  current  is  enhanced  by 
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atmospheric  ionization,  and  is  accompanied  by  intense  optical  emissions  known  as  sprites, 
only  for  the  most  intense  discharges  which  transport  hundreds  of  coulombs  of  charge  from 
cloud-to-ground. 

The  conservation  of  current  within  the  streamer  as  it  propagates  in  the  atmosphere  requires 
focusing  of  the  streamer  propagating  downward  and  expansion  of  the  streamer  propagating 
upward  as  is  discussed  in  more  detail  in  the  next  chapter  (Section  5.3).  The  most  probable 
physical  process  which  may  stop  the  focusing  of  a  downward  propagating  ionization  channel 
is  diffusion.  For  example,  the  ambipolar  diffusion  of  electrons  may  stop  the  focusing  process 
as  soon  as  the  column  has  a  cross  section  <  1  O^rrr  at  ~50  km  and  may  lead  to  the  splitting  of 
the  sprite  into  many  small  scale  filaments.  The  splitting  process  cannot  be  properly  modeled 
with  our  a7imuthally  symmetric  model;  nevertheless  the  model  results  are  in  reasonable 
agreement  with  the  dynamics  of  the  decay  of  luminosity  which  first  starts  at  high  altitudes 
(Figure  4.2c,  top  panel)  as  observed  by  Rairden  and  Mende  [1995]. 

4.7  COMPARISON  WITH  OBSERVATIONS  OF  SPRITES 


4.7.1  Comparison  With  Video  Observations 

Figure  4.4b  shows  the  expected  red  and  blue  video  response  at  t=20  ms  for  case  e,  calcu¬ 
lated  using  the  procedure  described  in  Section  3.5.  The  blue  color  becomes  increasingly 
significant  at  lower  altitudes,  being  only  a  factor  of  2  below  the  red  at  ~50  km.  The  range 
of  the  optical  emission  intensity  ( 1 02- 1 0^  kR)  shown  in  Figure  4.4b  is  consistent  with  video 
observations  [Sentman  et  al.,  1995]. 

The  quenching  altitude  of  the  1st  positive  band  of  N2  is  ~53  km  and  that  of  the  2nd 
positive  band  is  ~30  km  [Valance  and  Jones,  1974,  p.  1  19].  Additional  calculations  (see, 
for  example,  results  for  blue  jets  for  altitude  range  20-50  km,  in  the  next  chapter)  indicate 
that  blue  can  become  the  dominant  color  below  ~  40  km  due  to  the  quenching  effects 
combined  with  particular  properties  of  video  filters  which  generally  overemphasize  blue 
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Fig.  4.4.  Optical  emission  intensity  for  selected  bands  and  corresponding  video  re¬ 
sponse.  Altitude  scans  at  (r=0)  of  optical  emission  intensity  /  corresponding  to  different 
optical  bands;  (b)  Expected  red  and  blue  video  response  corresponding  to  optical  emissions 
shown  in  (a). 


color  (Section  3.5),  in  agreement  with  video  observations  of  bluish  lower  altitude  tendrils 
in  most  of  sprites  [Serttrnan  et  al.,  1995], 

The  variety  of  sprite  shapes  shown  in  Figure  4.3  are  within  the  range  of  the  measured 
transverse  and  horizontal  sizes  of  sprites.  However,  the  five  cases  shown  in  Figure  4.3  do 
not  represent  the  complete  diversity  of  shapes  which  have  been  experimentally  documented 
or  which  can  occur  under  different  conditions.  Model  calculations  (not  shown)  indicate  that 
the  transverse  size  of  sprites  can  be  very  large  (e.g.,  ~  50x50  km)  depending  primarily  on 
the  total  amount  of  charge  removed  from  the  thundercloud  (i.e.,  effective  volume  of  space 
in  which  the  post  discharge  electric  field  exceeds  the  characteristic  breakdown  field  E ;.),  the 
dynamics  of  charge  removal,  and  the  ambient  conductivity  profile.  Localized  conductivity 
enhancements  (e.g.,  case  b)  which  occur  due  to  external  sources  (i.e.,  electron  precipitation 
from  the  magnetosphere,  meteors,  runaway  electrons,  density  perturbations  left  from  a 
previous  sprite,  etc.)  may  significantly  enhance  the  electric  field  and  control  the  transverse 
scale  of  ionization  channels  and  may  lead  to  the  creation  of  clusters  of  sprites  by  a  single 
cloud-to-ground  discharge.  Under  certain  conditions,  such  preexisting  enhancements  of 
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conductivity  may  stimulate  both  upgoing  (toward  the  ionosphere)  and  downgoing  channels 
of  ionization. 

Video  observations  [Rairden  and  Mende,  1995]  suggest  the  duration  of  sprites  may  be  as 
large  as  ~40  ms.  On  the  other  hand,  photometric  measurements  of  [Fukunishi  et  al,  1 995] 
suggest  that  the  full-width-half-  maximum  duration  is  ~2  ms.  The  longer  persistence  of  the 
optical  emissions  observed  in  video  may  be  due  to  the  spectral  response  of  a  video  camera  in 
the  near  infrared  region  of  the  spectrum  where  the  emission  from  the  sprite  might  last  longer 
in  time  than  in  the  visible  range  [Winckler  et  al,  1996],  consistent  with  the  fact  that  the 
threshold  of  excitation  for  the  red  emission  (N2  1  st  positive  band)  is  lower.  Another  reason 
for  the  persistence  of  optical  emissions  might  be  connected  with  the  flow  of  the  current  in 
the  body  of  the  ionization  channel.  The  conduction  current  always  exist  in  the  body  of  the 
ionization  channel  as  it  propagates  due  to  the  finite  electric  field  E  <  E *  (e.g..  Figure  4.2b, 
t  =2,10,15,20  ms),  and  decays  as  soon  as  the  channel  stops  (e.g.,  Figure  4.2b,  t—  32  ms). 
The  observed  persistence  of  the  luminosity  of  sprites  for  several  tens  of  ms  [Rairden  and 
Mende ,  1995]  might  be  a  strong  indication  of  the  persistence  of  the  electric  field  and  the 
associated  current  flow  in  the  body  of  the  sprite  (Section  4.6). 

4.7.2  Comparison  With  Photometric  Observations 

Figure  4.5  shows  results  of  the  spatially  integrated  optical  emission  intensity  as  a  function 
of  time  t  for  the  different  cases  of  charge  removal  a,b,c,d,  and  e  shown  in  Figure  4.1 .  The 
onset  of  the  optical  emission  does  not  necessarily  coincide  with  the  onset  of  the  causative 
cloud-to-ground  discharge  (/=0)  and  can  in  fact  be  significantly  delayed  in  time  (up  to  ~20 
ms  for  case  d).  This  result  is  in  agreement  with  photometric  measurements  by  Winckler  et  al 
[  1 996]  who  reported  delays  of  ~3  ms,  and  by  Fukunishi  et  al  [1996]  who  reported  delays 
ranging  from  several  ms  up  to  several  tens  of  ms  (see  Section  1.2.7).  Physically,  this  effect 
can  be  understood  in  terms  of  the  threshold  dependence  on  the  electric  field  of  the  intensity  of 
optical  emission  and  ionization  processes.  Optical  emissions  appear  nearly  instantaneously 
in  time  when  the  corresponding  threshold  is  reached.  Results  in  Figure  4.5  for  cases  a  and 
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e  show  that  short  pulses  with  ~3  ms  duration  which  immediately  follow  the  causative  cloud 
to  ground  discharge  are  created  by  the  initial  rapid  removal  of  200  C  of  charge  in  1  ms. 
These  few  ms  pulses  are  different  than  the  shorter  duration  (few  hundred  //s)  ionospheric 
flashes  called  “elves”  which  are  likely  to  be  produced  by  heating  of  the  ambient  electrons 
by  electromagnetic  pulses  (EMPs)  radiated  by  the  discharge  and  are  observed  within  a  few 
hundred  microseconds  of  the  causative  cloud-to-ground  discharge  [Fukunishi  et  al.,  1996; 
hum  et  al.,  1996b].  As  was  discussed  in  Section  2.1,  the  EMP-induced  interactions  are  not 
included  in  our  quasi-static  formulation. 


Fig.  4.5.  Spatially  integrated  optical  emission  intensity  as  a  function  of  time.  Spatial l> 
integrated  optical  emission  intensity  of  the  first  positive  band  of  N2  I  as  a  function  of  time 
t  corresponding  to  five  models  of  thundercloud  charge  removal  described  in  Figure  4.1 


4.7.3  Comparison  With  Spectroscopic  Observations 

The  altitude  distribution  of  emission  intensities  (at  r=0)  corresponding  to  five  different 
optical  bands  are  shown  in  Figure  4.4a  for  case  e  and  at  (=20  ms.  The  emission  levels  in 
these  bands  depend  on  the  value  of  the  electric  field  (i.e.,  the  level  of  electron  heating)  and 
the  electron  density  (ambient  plus  that  produced  due  to  ionization)  and  might  be  different 
for  different  altitudes  and  regimes  of  charge  removal.  However,  our  calculations  show  that 
under  quasi-electrostatic  heating  conditions  the  Meinel  band  of  Nt  and  1st  negative  band 
of  Ctf  are  typically  excited  at  intensities  well  below  those  of  the  ambient  nighttime  airglow 
(~104  R  and  ~102  R,  respectively  [Chamberlain,  1978,  pp.  214-218]).  The  first  and 
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second  positive  bands  of  N2  and  the  1st  negative  band  of  N£  have  ambient  levels  of  ~104 
R,  ~5  R,  and  ~5  R,  respectively  [Chamberlain,  1978,  pp.  214-218],  and  under  different 
quasi-electrostatic  heating  conditions  can  be  excited  at  levels  well  above  the  ambient  levels. 

Recent  spectral  measurements  of  sprites  confirmed  the  dominant  presence  of  the  first 
positive  band  of  N2  [ Mende  et  al„  1995;  Hampton  et  al.}  1996].  Most  of  the  energy  of 
the  second  positive  band  of  N2  is  concentrated  in  the  wavelength  region  of  ~300-400  nm 
[Vallance  Jones,  1 974,  Table  4. 1 1 ,  p.  131]  (see  also  Figure  3.7)  which  is  below  the  spectral 
range  resolved  in  Mende  et  al.  [1995]  (450-800  nm)  and  Hampton  et  al.  [1996]  (540- 
840  nm)  measurements,  consistent  with  the  fact  that  this  band  was  not  detected.  As  was 
discussed  in  Section  1.2.8,  with  the  available  signal-to-noise  ratio  of  the  measurements  of 
Mende  et  al.  [1995]  and  Hampton  et  al.  [1996],  in  order  to  be  detectable  other  spectral 
bands  would  have  had  to  be  at  levels  no  less  than  a  factor  of  ten  of  that  of  the  N2  first  positive 
band.  In  accordance  with  our  model  results,  the  intensities  of  the  Meinel  and  1st  negative 
bands  of  N£  and  the  1st  negative  band  of  Oj  are  generally  produced  at  levels  several  orders 
of  magnitude  lower,  also  fully  consistent  with  the  fact  that  these  bands  were  not  detected 
by  Mende  et  al.  [  1 995]  and  Hampton  et  al.  [  1 996]. 


4.8  SPRITES  AND  RUNAWAY  ELECTRONS 

The  possible  role  of  runaway  electrons  in  producing  X-rays  that  are  sometimes  observed 
prior  to  lightning  discharges  and  in  the  initiation  of  the  lightning  discharge  itself  was  consid¬ 
ered  by  McCarthy  and  Parks  [  1 992],  Gurevich  et  al.  [  1 992, 1 994]  and  Roussel-Dupre  et  al. 
[1994]  (see  also  references  cited  therein),  who  showed  that  secondary  cosmic  ray  electrons 
of  energy  ~1  MeV  can  become  runaways  whenever  the  local  electric  field  E  exceeds  the 
threshold  field  Et  =  ( 4TTNmZqja)/(mc 2),  where  Nm  is  the  number  density  of  airmolecules, 
Z  ~  1 4.5,  a  ~1 1,  m  is  the  electron  rest  mass,  and  c  is  the  velocity  of  light  [Gurevich  et  al., 
1994]. 


Bell  et  al.  [  1 995]  considered  the  effects  of  an  upward  traveling  beam  of  ~  1  MeV  runaway 
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electrons  accelerated  by  large  quasi-electrostatic  fields  above  thunderclouds  as  a  possible 
mechanism  for  production  of  sprites.  The  main  energy  loss  for  the  runaway  electrons  occurs 
from  inelastic  collisions  with  O2  and  NS  molecules,  resulting  in  optical  excitation  of  the 
neutrals,  and  subsequent  optical  emissions.  Results  of  a  one-dimensional  model  suggest 
that  the  runaway  electrons  can  produce  optical  emissions  similar  in  intensity  and  spectra  to 
those  observed  in  sprites,  but  only  for  large  quasi-electrostatic  fields  produced  by  positive 
cloud-to-ground  discharges  lowering  250  C  or  more  to  ground  from  an  altitude  of  at  least 
10  km.  In  general,  the  threshold  runaway  field  Et  is  ~10  times  lower  than  the  threshold 
field  of  the  conventional  breakdown  Ef.  (see  Section  2.4.2)  and  on  this  basis  alone  we 
might  expect  the  runaway  mechanism  to  proceed  at  lower  quasi-electrostatic  field  levels 
[Taranenko  and  Roussel-Dupre,  1996].  However,  in  the  lower  atmosphere,  the  runaway 
avalanche  breakdown  process  starts  from  very  low  ambient  production  rates,  requiring 
many  e-folding  distances  before  observable  effects  can  be  produced  [Belief  al.,  1995].  As  a 
result,  production  of  significant  optical  emission  levels  via  the  runaway  mechanism  requires 
charge  transfers  similar  to  quasi-electrostatic  heating  mechanism  [Bell  et  al,  1995]. 

Figure  4.6a  from  [Bell  et  al.,  1 995]  shows  the  optical  intensity  per  unit  volume  of  the  1  st 
positive  band  of  NS  produced  by  the  runaway  breakdown  process.  The  five  curves  shown 
correspond  to  specific  values  of  the  total  charge  Q  removed  from  a  thundercloud  at  10  km 
altitude  in  1  ms  and  indicate  that  the  optical  output  sensitively  depends  on  the  removed 
charge  value.  The  same  result  holds  for  all  the  other  excited  bands,  not  shown  here  for  the 
sake  of  brevity. 

Figure  4.6b  shows  the  optical  intensity  per  unit  volume  of  the  1st  and  2nd  positive  bands 
of  Nt,  the  Meinel  and  1st  negative  bands  of  N,,  and  the  1st  negative  band  of  O,  produced 
by  runaway  electrons  accelerated  by  quasi-electrostatic  thundercloud  fields  for  Q  =275  C 
[Bell  et  al.,  1995].  The  two  positive  bands  of  N;  dominate  the  spectrum  up  to  ~80  km. 
above  which  the  Meinel  band  of  Nt  also  becomes  important.  Video  observations  indicate 
that  the  brightest  region  of  a  sprite  lies  between  66  and  74  km  and  that  some  sprites  exhibit 
significant  intensity  at  altitudes  as  low  as  40-50  km  [e.g.,  Sentman  et  al.,  1 995].  As  shown 
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1st  positive  band  of  N2 


Fig.  4.6.  Optical  emissions  produced  by  runaway  electrons,  (a)  The  volume  emission 
rale  (in  photons  cm-3s-1)  profile  of  lhe  lsl  posilive  band  of  N2  for  selecled  Ihundercloud 
charge  values,  and  (b)  lhe  volume  emission  rale  of  selecled  bands  for  N2,  N£  and  0£  for 
a  thundercloud  charge  of  275  C,  obtained  from  runaway  electron  model;  (c)  The  optical 
emission  intensity  (in  Ray  leighs)  profiles  of  lhe  firsl  positive  band  of  N2  obtained  from 
quasi-electrostatic  heating  model  for  differenl  thundercloud  charge  values.  For  a  typical 
~10  km  (~106  cm)  diameter  of  a  sprite  the  horizontal  scale  in  (a)  and  (b)  can  be  interpreted 
as  lhe  inlensily  in  Rayleigh  (see  (2.1 1)).  Figure  is  reproduced  from  [Bell  et  al,  1995). 


in  Figure  4.6,  the  optical  emission  intensities  produced  by  the  runaway  breakdown  process 
tend  to  peak  in  the  55-70  km  altitude  range,  and  substantial  intensity  can  be  found  at  altitudes 
as  low  as  45  km. 

Figure  4.6c  shows  the  intensity  of  the  first  positive  band  of  N2  produced  by  the  quasi - 
electrostatic  heating  mechanism  for  different  values  of  the  thundercloud  charge  Q  removed 
from  10  km  altitude  in  1  ms.  These  simulation  conditions  are  identical  to  those  used  to  obtain 
results  of  the  runaway  model  shown  in  Figure  4.6a.  Two  important  features  are  apparent: 
Firstly,  the  runaway  process  produces  optical  emissions  which  span  a  wider  range  of  altitudes 
than  those  produced  by  the  quasi-electrostatic  heating;  Secondly,  the  emission  intensities 
generated  in  the  runaway  process  sensitively  depends  on  the  removed  thundercloud  charge 
value  and  become  negligible  for  charges  Q  <250  C. 

Bell  et  al.  [1995]  recognized  two  other  notable  differences  between  the  runaway  break¬ 
down  and  quasi-electrostatic  heating  model:  i)  the  ratio  R\2  of  the  intensity  of  the  1st 
positive  band  of  N2  to  the  2nd  positive  band  of  N2  is  ~2  at  all  altitudes  above  60  km  for 
the  runaway  model,  while  for  the  quasi-electrostatic  heating  model  R\2  —10  above  80  km 


SPRITES  AS  LUMINOUS  COLUMNS  OF  IONIZATION 


87 


altitude  (see  Figures  3.9b  and  4.4a);  and  ii)  in  the  runaway  model  the  intensity  of  the  Meinel 
band  of  Nt  and  the  2nd  positive  band  of  N2  are  roughly  equal  at  80  km  altitude,  while  for 
the  quasi-electrostatic  heating  mechanism  the  intensity  of  the  Meinel  band  of  Nj  is  smaller 
than  that  of  the  2nd  positive  band  of  N2  by  a  factor  of  ~  1 (see  Figure  3.9b).  As  a  result  of 
these  differences,  Bell  el  al.  [1995]  suggested  that  spectral  analysis  of  sprites  may  delineate 
the  conditions  under  which  the  two  different  generation  mechanisms  are  operative.  For 
example,  altitude  profiles  of  spectral  ratios  could  determine  the  regions  dominated  by  each 
mechanism. 

As  a  discussion  point  and  a  subject  for  future  studies  it  is  important  to  realize  that  the 
temporal  behavior  of  the  runaway  electrons  may  significantly  complicate  the  interpretation 
of  spectral  measurements  of  sprites  (see  Section  1.2.8).  For  example,  impulsive  (e.g.,  <  1 
ms)  production  of  a  large  number  of  energetic  secondary'  electrons  would  lead  to  the  rapid 
cooling  of  these  electrons  (e.g..  for  times  >  1  ms).  Since  the  electric  field  values  required  for 
the  runaway  process  are  comparable  with  those  for  conventional  quasi-electrostatic  heating 
these  cooled  electrons  might  remain  at  a  heated  level  defined  by  the  quasi-electrostatically 
heated  electron  distribution,  instead  of  the  energetic  distribution  initially  created  by  the 
runaway  beam.  The  recent  report  of  Winckler  el  al.  [1996]  that  sprites  begin  with  an 
intense,  w  ideband  flash  covering  both  the  blue  and  red  portions  of  the  spectrum,  with  a  fast 
decaying  blue  portion  (see  Section  1.2.7  for  more  details),  may  indicate  the  presence  of 
a  cooling  process  in  which  initially  energetic  secondary  electrons  produced  by  a  runaway 
beam  are  cooled  to  the  energy  level  of  conventional  heating  as  described  above. 

Although  in  the  stationary  case  the  secondary  electron  distribution  produced  by  runaway 
electrons  and  the  resultant  complexity  of  optical  emissions  are  similar  to  that  created  by 
energetic  precipitating  particles  in  aurora  (see  Appendix  B  )  we  note  that  in  view  of  our 
previous  arguments  the  impulsive  nature  of  the  runaway  process  may  lead  to  a  highly  variable 
(on  a  time  scale  of  ~  1  ms)  spectral  content.  This  subject  requiress  further  investigation  (see 
Section  6.2). 


5 

Blue  Jets  Produced  by  Pre-Discharge 
Thundercloud  Fields 


In  this  chapter  we  show  that,  under  certain  conditions,  pre-discharge  quasi-electrostatic 
thundercloud  fields  immediately  above  the  thundercloud  can  lead  to  the  formation  and 
upward  propagation  of  streamer  type  ionization  channels  with  features  in  good  agreement 
w  ith  recent  video  observations  of  blue  jets. 


5.1  PHYSICAL  MECHANISM 

We  propose  that  blue  jets  are  also  associated  with  streamer  type  processes  occurring  on  an 
atmospheric  spatial  scale.  The  difference  between  sprites  and  jets  in  the  framework  of  our 
quasi-electrostatic  heating  model  can  be  understood  as  follows: 

Sprites  are  relatively  fast  processes  (time  scale  ~1-10  ms)  occurring  at  mesospheric 
altitudes  due  to  the  impulsive  appearance  of  large  quasi-electrostatic  fields  follow  ing  cloud- 
to-ground  lightning  discharges.  The  time  scale  of  sprites  and  their  energetics  are  thus  closely 
associated  with  properties  of  the  causative  cloud-to-ground  discharge.  Ionization  regions 
associated  with  sprites  enhance  atmospheric  conductivity  and  serve  as  additional  factors 
leading  to  more  rapid  relaxation  of  the  post-discharge  electric  fields.  Jets  are  much  slower 
processes  (time  scale  ~300  ms)  which  are  not  associated  with  cloud-to-ground  lightning 
discharges  and  can  be  initiated  near  the  thundercloud  tops  in  cases  of  slow  accumulation 
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of  large  amounts  of  charge  at  altitudes  ~20  km.  The  upward  propagation  of  jets  can  only 
be  sustained  if  thundercloud  acts  as  a  “battery”  which  is  energetic  enough  to  continuously 
compensate  lost  of  thundercloud  charge  due  to  ionization  and  conduction  current  associated 
with  jet.  Slower  time  scales  of  jets  thus  are  mostly  controlled  by  the  dynamics  of  charge 
separation  processes  inside  the  cloud. 


5.1.1  Blue  Jets  as  Streamers 

As  was  already  discussed  in  the  previous  Chapter,  streamers  can  be  defined  in  general  terms 
as  transient  filamentary  plasmas,  the  dynamics  of  which  are  controlled  by  highly  localized 
nonlinear  space  charge  waves  [see  Vitello  et  al.,  1994;  Grange  et  al.,  1995;  Lagarkov  and 
Rutkevich ,  1994,  and  references  therein].  The  morphology  of  streamers  is  very  complex 
even  for  the  simplest  geometries  [e.g.,  Vitello  et  al.,  1994].  However,  a  series  of  rigorous 
numerical  simulations  [Dhali  and  Williams,  1987;  Vitello  et  al.,  1994]  has  revealed  several 
robust  features  of  streamers  which  are  discussed  below  in  terms  of  their  application  to 
blue  jets.  The  streamer  velocity  vs  usually  lies  in  the  range  ~102-103  krn/s  depending  on 
the  streamer  polarity,  the  applied  electric  field  ( E ),  the  initial  radius,  the  density  of  the 
seed  ionization  which  initiates  the  streamer  process,  the  background  electron  density  [e.g., 
Dhali  and  Williams,  1987],  and  the  geometry  of  the  boundaries  (e.g.,  point  to  plane  corona 
discharge)  [Grange  et  al.,  1995],  The  polarity  of  a  streamer  is  defined  by  the  sign  of  the 
space  charge  induced  in  its  “head”.  A  negative  streamer  propagates  due  to  ejection  of 
electrons  from  its  head.  Positive  streamers  require  a  population  of  ambient  electrons  (e.g., 
produced  by  cosmic  rays)  for  their  development,  have  relatively  low  velocities  of  ~  100 
km/s,  comparable  to  electron  drift  velocities  (vd)  in  the  body  of  the  streamer,  and  occur 
for  low  applied  electric  fields  [Dhali  and  Williams,  1987]  and  in  point-to-plane  corona 
discharges  [e.g.,  Grange  et  al.,  1995].  The  minimum  propagation  velocity  of  the  positive 
streamers  is  equal  to  100  km/sec  in  air  [e.g.,  Petrov  et  al.,  1994].  The  streamer  velocity 
generally  depends  on  the  steepness  of  the  applied  voltage  pulse  in  time,  remains  ~  100  km/s 
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for  relatively  low  voltage  rise  rates  (<  50  kV/jis)  but  becomes  larger  for  rise  rates  >50 
kV//xs  [e.g.,  Petrov  et  al.,  1994], 


Assuming  that  the  radius  rs  of  the  streamer  column  is  approximately  equal  to  that  of  its 
head,  the  streamer  velocity  vs,  its  radius  rs  and  the  dielectric  relaxation  time  rs  within  the 
streamer  ( rs  =  £0/vs,  where  os  is  the  streamer  conductivity)  are  related  simply  by  [Vitello 
et  al.,  1994]: 


b_ 

rs 


Vs  Ts 


Et, 

Eh 


(5.1) 


where  E/t,  Eh  are  the  electric  fields  at  the  streamer  head  and  in  the  streamer  body,  respec¬ 
tively,  and  b  =  vsrs  is  an  effective  “skin  depth”  for  the  electric  field  inside  the  streamer.  In 
most  cases  0. 1<  Eh/ Eh  <  1  so  that  6  is  comparable  to  the  streamer  radius  rs  [ Lagarkov 
and  Rutkevich,  1 994,  p.  1 6].  Since  diffusion  effects  do  not  appear  to  play  a  significant  role  in 
the  formation  of  the  streamer  [Dhali  and  Williams,  1987;  Vitello  et  al.,  1994],  the  effective 
skin  depth  b  is  the  only  parameter  which  characterizes  the  spatial  scale  of  the  streamer.  On 
this  basis,  we  suggest  that  streamer  type  solutions  are  possible  for  a  variety  of  transverse 
dimensions,  including  the  kilometer  spatial  scales  typical  for  blue  jets,  provided  that  the 
appropriate  initial  conditions  can  occur  near  the  top  of  the  thundercloud.  In  this  context, 
blue  jets  are  similar  to  previously  studied  positive  streamers  which  initiate  point-to-plane 
corona  discharges  in  laboratory  experiments  [e.g..  Grange  et  al.,  1995],  except  for  the 
differences  in  the  electric  field  production,  geometry  and  maintenance  of  the  field  by  the 
thundercloud  charge  distribution,  and  the  altitude  variation  of  atmospheric  neutral  density 
and  conductivity. 


5.1.2  Formation  of  Blue  Jets 

The  positive  (top)  and  the  negative  (bottom)  thundercloud  charges  (Figure  5.1)  slowly 
(seconds)  accumulate  due  to  the  current  J's  (see  Section  2.1)  associated  with  the  separation 
of  charges  inside  the  cloud,  and  directed  opposite  to  the  resulting  electric  field.  The  current 
may  be  related  to  the  small  and  light  positively  charged  ice  splinters  driven  by  updrafts,  and 
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heavy  negatively  charged  hail  particles  driven  downward  by  gravity  [e.g.,  liman,  1987,  p. 
65].  Unusually  intense  hail  activity  was  indeed  observed  in  association  with  jets  [  Wescott 
et  al.,  1 9956;  1 996],  and  is  a  strong  indication  of  intense  electrical  activity  inside  the  cloud. 


Electrosphere  r  -v 
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Fig.  5.1.  Mechanism  of  Blue  Jets.  Currents,  charges  and  electric  fields  associated  with 
the  blue  jet.  The  Figure  is  partly  adapted  from  [ Uman ,  1987,  p.  30]. 


Under  certain  relatively  rare  conditions,  a  charge  of  (^=300-400  C  can  accumulate  in 
thunderclouds  at  ~20  km  altitude.  This  charge  may  have  an  extended  dish-like  distribution 
and  can  lead  to  breakdown  ionization  in  regions  with  lateral  extent  102-103  m.  The 
ionization  may  also  be  connected  with  the  accumulation  of  long  lived  (~l-2  sec  at  ~20  km 
[Lowke,  1992])  metastable  O2  molecules  due  to  upward  directed  discharges  preceding  the 
development  of  jets  [Wescott  et  al.,  1995a],  corona  discharges  initiated  by  hydrometeors 
[Griffiths  and  Phelps ,  1976]  or  highly  charged  hailstones  triggering  instability  of  charged 
droplets  [Grigor’ev  and  Shiryaeva,  1989].  Once  formed,  the  ionized  region  can  act  as 
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the  seed  for  the  development  of  an  upward  propagating  positive  streamer.  We  propose 
that  jets  are  the  mesospheric  analog  of  the  positive  type  streamers  produced  in  laboratory 
point-to-plane  corona  discharge  experiments. 

The  electric  field  is  non  zero  inside  the  streamer  as  it  propagates  and  an  upward  directed 
current  (Is)  always  flows  in  the  streamer  body  (Figure  5.1).  The  streamer  propagates 
freely  when  the  current  Is  is  maintained  by  a  source  in  or  near  the  cloud.  Otherwise,  the 
negative  charge  flowing  within  the  streamer  body  towards  the  positive  thundercloud  charge 
would  reduce  the  source  charge  and  the  electric  field  and  would  eventually  suppress  the 
propagation  of  the  streamer.  Thus,  the  streamer  can  propagate  as  long  as  the  source  current 
J'  can  deliver  sufficient  positive  charge  to  the  top  of  the  thundercloud.  An  equal  amount 
of  negative  charge  is  accumulated  at  the  thundercloud  base  so  that  the  overall  charge  in  the 
cloud-streamer  system  is  conserved.  The  thundercloud  thus  plays  the  role  of  a  “battery” 
and  the  net  effect  of  the  streamer  attached  to  this  “battery”  is  to  deliver  positive  charge  to 
the  conducting  electrosphere,  contributing  to  the  maintenance  of  the  fair  weather  potential. 
The  streamer  propagation  can  also  cease  when  it  reaches  the  critical  altitude  at  which  the 
dielectric  relaxation  time  rs  within  the  streamer  body  is  comparable  to  the  relaxation  time  r 
of  the  ambient  conducting  atmosphere  (r  =  e0/o).  At  this  altitude,  all  of  the  current  flowing 
through  the  channel  created  by  the  streamer  can  be  supplied  by  the  flow  of  atmospheric  ions, 
and  the  electric  field  at  the  tip  of  the  streamer  is  not  enhanced  since  the  ambient  conductivity 
a  is  comparable  to  the  conductivity  as  inside  the  streamer  head. 


5.2  MODEL  PARAMETERS 

For  calculations  in  this  Chapter  we  use  the  ambient  electron  density  profile  1  shown  in 
Figure  2.4a  and  the  ion  conductivity  profile  A  illustrated  in  Figure  2.4b. 

In  our  model,  the  effective  source  current  J'  (see  Figure  5. 1  and  discussion  in  Section  2. 1 ) 
is  introduced  in  such  a  way  as  to  compensate  any  changes  in  the  thundercloud  charge  ps 
due  to  the  conduction  losses  (which  in  the  case  of  blue  jets  also  include  the  current  flow  and 
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associated  charge  accumulation  due  to  the  streamer  attached  to  the  top  of  the  cloud).  As 
discussed  in  Section  2.1  the  introduction  of  J's  allows  us  to  separately  specify  the  charge 
dynamics  inside  of  the  cloud  as  a  given  function  of  time. 

We  assume  that  a  positive  thundercloud  charge  is  located  at  20  km  altitude,  has  a  spher¬ 
ically  symmetric  Gaussian  distribution  with  a  spatial  scale  of  a=3  km,  builds  up  over  a 
time  scale  ry,  and  is  dissipated  over  a  time  of  r^.  This  positive  charge  is  presumed  to  be 
accompanied  by  negative  charge  at  lower  altitudes,  as  shown  in  Figure  5.1.  In  determining 
the  dynamics  of  the  resulting  blue  jet,  we  neglect  the  effects  of  the  lower  thundercloud 
charge  due  to  its  proximity  to  the  ground.  This  negative  charge  may  in  fact  also  be  removed 
by  a  series  of  negative  lightning  discharges  which  sometimes  are  observed  several  seconds 
before  the  appearance  of  blue  jets  [Wescott  et  al.,  1 9956;  1 996]. 

Equations  (2.2)  and  (2.3)  which  are  used  for  the  calculations  of  the  electric  field  and 
charge  density  are  contained  in  the  more  general  two-fluid  streamer  equations  [e.g.,  Vitello 
et  al.,  1 994],  In  comparison,  equation  (2.3)  does  not  contain  a  term  describing  exponential 
growth  of  the  electron  density  due  to  ionization  processes  (see,  for  example,  equation  (2.9) 
in  Section  2.4.4).  The  increase  in  the  electron  density  Ne  is  self-consistently  included  in 
our  model  through  changes  in  conductivity  a.  However,  since  for  typical  parameters  used 
we  have  vj  ~  vs,  our  formulation  cannot  be  used  to  accurately  determine  the  small  scale 
details  of  the  electron  number  density  Ne  near  the  steep  boundary  at  the  streamer  head. 
Results  of  extensive  modeling  [e.g.,  Vitello  et  al.,  1994]  indicate  that  for  a  well  developed 
streamer,  the  relaxation  time  rs  behind  the  streamer  front  remains  approximately  constant. 
We  use  this  result  as  an  additional  external  condition  in  our  model,  placing  a  bound  on  the 
growth  of  the  electron  density  Nc  in  the  streamer  front. 

A  characteristic  breakdown  electric  field  E%  =  3.7  x  106iV/iVo  V/m  (N0  =  2.688 x  1025 
m-3)  is  defined  for  calculations  in  this  chapter  as  the  field  at  which  the  effective  ionization 
coefficient  Vi  is  equal  to  the  attachment  coefficient  ua  (see  Figure2.6),  ~  1 5%  larger  than  the 
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characteristic  air  breakdown  value  E =  3.2  x  \0(>N/No  V/m  defined  by  Papadopoulos  et 
al.  [  1 993]  and  used  in  sprite  ionization  column  calculations  reported  in  the  previous  chapter. 

For  electric  fields  E  >  E£,  the  electron  density  Ne  grows  exponentially  in  time  with 
a  rate  (vj-va)>0  determined  as  a  function  of  the  electric  field  by  the  equations  (2.6)  and 
(2.7).  For  electric  fields  E  <  E%  (i.e.,  (ly-i^cO  as  is  the  case  just  behind  the  streamer 
front)  electrons  may  be  quickly  removed  by  dissociative  attachment  to  CF.  However,  such 
may  not  always  be  the  case  due  to  the  extremely  effective  detachment  process  of  electrons 
in  the  presence  of  metastable  O2  [Lowke,  1992],  leading  to  a  reduction  of  the  breakdown 
field  Ek  by  factor  of  about  6,  in  some  cases.  As  already  discussed  in  Section  2.4.3,  the 
role  of  dissociative  attachment  of  electrons  to  O2  is  difficult  to  assess  since  the  measured 
rates  for  associative  detachment  and  for  collisional  detachment  are  known  only  within  one 
to  three  orders  of  magnitude,  respectively  [Davies,  1983].  In  view  of  the  complexity  of  the 
role  of  detachment,  we  conservatively  assume  for  calculations  in  this  chapter  that  (ui-ua)=0 
for  electric  fields  E  <  E '£.  In  addition  to  the  above  arguments,  the  more  conservative 
threshold  level  of  E%  (i.e.,  3.7x  106AT/./Vo  instead  of  3.2x  \0bN/No)  is  used  here  for  blue 
jets  calculations  in  view  of  the  sensitive  dependence  of  the  optical  excitation  coefficients 
on  the  electric  field  at  low  values  of  E.  The  chosen  value  of  E^  places  the  optical  emission 
production  at  electric  field  levels  where  the  ratio  of  the  excitation  coefficients  of  the  1  st  and 
2nd  positive  bands  of  N2  is  ~10  rather  than  several  orders  of  magnitude  (Figure  2.7). 

Model  calculations  indicate  that  the  streamer  morphology  above  a  thundercloud  can 
be  complex,  exhibiting  large  variability  in  lateral  extent  rs,  stopping  altitude  and  vertical 
velocity  of  propagation  vs,  depending  mostly  on  the  initial  conditions  and  the  source  function 
(e.g.,  the  current  J's).  Such  variability  is  in  general  agreement  with  the  degrees  of  freedom 
allowed  for  streamer  solutions  in  equation  (5.1)  and  with  recent  observations  of  a  new 
phenomenon  known  as  “blue  starters”  [Wescott  et  al.,  1995ft;  1996],  which  exhibit  shorter 
lifetimes,  wider  range  of  velocity  variations  (~27- 1 57  km/s)  and  smaller  transverse  and 
vertical  extents  in  comparison  with  jets  [Wescott  et  al.,  1995a;  1996]. 
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We  assume  tj= 1  sec  and  r^= 0.5  sec,  and  consider  rs  as  an  input  parameter.  The  initial 
ionization  region  and  the  initial  radius  of  the  jet  rs  are  then  defined  mostly  by  the  effective 
volume  of  space  over  which  the  electric  field  E  exceeds  E%,  depending  on  the  magnitude 
and  distribution  of  the  thundercloud  charge  Q.  The  initial  ionization  that  forms  naturally 
above  the  thundercloud  plays  the  same  role  as  the  seed  ionization  in  streamer  simulations 
[e.g.,  Vitello  et  al.,  1994], 

Our  model  results  indicate  that  the  dependence  of  the  streamer  characteristics  on  rs  and 
Q  are  in  agreement  with  (5.1 ).  The  model  solutions  are  very  sensitive  to  initial  conditions 
and  demonstrate  large  variability  in  possible  outputs.  This  sensitive  dependence  makes  it 
difficult  to  identify  a  physical  mechanism  which  automatically  and  selfconsistently  selects 
initial  conditions  characteristic  for  blue  jets.  It  is  important  to  stress,  however,  that  when 
we  select  the  “correct”  initial  conditions,  model  calculations  reproduce  the  selfconsistent 
development  of  the  streamer  and  demonstrates  agreement  with  almost  every  observed  feature 
of  blue  jets.  Such  a  case  showing  agreement  w  ith  the  observed  characteristics  of  jets  [Wescott 
et  al,  1995a]  is  described  in  the  following  section. 

5.3  ELECTRON  NUMBER  DENSITY,  ELECTRIC  FIELD  AND  CURRENT 

Figure  5.2  shows  the  electric  field  and  the  electron  number  density  Ne  corresponding 
to  a  jet  produced  by  a  thundercloud  charge  Q  =375  C  assuming  rs=0.05  sec.  Altitude 
scans  of  the  electric  field  (at  r=0)  are  shown  in  the  left  hand  panel  at  selected  instants 
of  time,  and  in  comparison  with  E The  enhancement  of  the  electric  field  at  the  tip  of 
the  streamer  is  apparent  in  the  middle  panel  of  the  Figure  5.2a.  The  characteristic  spatial 
scale  <5  is  comparable  with  rs.  The  electric  field  barely  exceeds  £*,  since  the  ionization 
coefficient  becomes  very  large  rs_1)  for  slight  increases  in  the  electric  field  above 

Ek-  Accordingly,  it  can  be  assumed  that  Eb  ~  Eh  in  (5.1).  The  electron  number  density  Ne 
shown  in  the  right  hand  panel  illustrates  the  formation  of  a  beam-like  distribution  expanding 
upward  with  a  radius  ~2  km  at  the  cloud  top  and  ~7  km  at  ~45  km  altitude. 
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Fig.  5.2.  Electric  field  and  electron  density  associated  with  a  blue  jet.  A  crossectional 
view  of  the  distribution  of  the  absolute  values  of  the  electric  field  E  and  the  electron  density 
Ne  at  selected  instants  of  time.  Right  hand  panel  shows  altitude  scans  at  r=0. 


The  results  of  the  numerical  model  can  be  qualitatively  interpreted  based  on  simple 
analytical  estimates.  Behind  the  front  of  the  jet  the  electric  held  E  ~  £*.  (Figure  5.2a), 
so  that  a  current  ~  qcNeV(j7rrl  flows  in  the  streamer  body  where  qe  is  the  electronic 
charge,  and  v d  =  / ieEk .  This  current  Is  is  estimated  to  be  ~200  A,  and  can  deliver  ~100C 
to  the  thundercloud  during  the  lifetime  of  the  jet.  The  electron  mobility  /j,e  for  E  ~  E is 
[je  ~  4  x  IO_2Ar„/A’’  (see  equation  (2.5)),  so  that  the  electron  drift  velocity  v<j  ~100  km/s, 
independent  of  altitude.  The  current  Is  is  conserved  along  the  jet.  Since  vj  is  constant,  the 
conservation  of  Is  implies  that  the  quantity  Ner 2  should  also  be  conserved.  In  our  model,  we 
assume  ts  =  £a/as  to  be  constant,  where  the  conductivity  in  the  streamer  body  as-qeNef.ic. 
which  is  equivalent  to  requiring  Nc  ~  N.  The  conservation  of  streamer  current  ls  thus 
implies  that  r,~Ar-7,  consistent  with  the  general  tendency  of  the  jet  to  expand  toward 
higher  altitudes,  in  good  agreement  with  our  numerical  results,  especially  below  40  km. 
Above  40  km,  the  expansion  slows  down  because  of  the  increasing  role  of  ionospheric 
ions  in  carrying  the  streamer  current,  as  well  as  the  reduction  in  the  source  charge  Q  (so 
that  the  current  Js  is  reduced  in  later  stages  of  streamer  propagation).  A  comparison  of 
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the  conditions  rs=const  and  rs~N~i  with  equation  (5.1)  shows  that  in  good 

agreement  with  our  numerical  results.  At  the  initial  stage  of  the  jet  formation  at  24  km 
altitude  we  have  km/s;  the  streamer  speed  subsequently  reaches  ~100  km/s  at  ~35 

km  altitude  and  remains  in  the  range  100-130  km/s  until  the  jet  stops.  The  numerical  values 
of  the  electron  density  and  streamer  radii  at  24  km  and  35  km  are  103  and  1 80  cm-3  and  2 
and  5  km,  respectively.  These  values  are  reasonably  consistent  with  the  conservation  of  the 
current  Is  and  the  condition  (5.1). 

5.4  OPTICAL  EMISSIONS 

The  optical  emission  intensity  of  the  second  positive  band  of  N2  is  shown  in  Figure  5.3. 
The  jet  propagates  upward  with  an  average  velocity  of  88  km/s  and  stops  at  ~50  km 
altitude  mainly  due  to  our  chosen  temporal  variation  of  the  source  thundercloud  charge 
Q(t)  with  decay  time  T<f= 0.5  sec.  If  the  charge  Q  is  instead  kept  constant  after  1  sec,  the 
jet  propagates  to  a  slightly  higher  altitude  (~53  km)  and  stops  when  rs  ~  e/ct,  which 
for  the  ion  conductivity  profile  used  occurs  at  49  km  altitude.  The  optical  emission  has 
the  shape  of  an  upward  expanding  beam  with  transverse  dimensions  being  essentially  the 
same  as  the  ionization  column  shown  in  Figure  5.3.  The  corresponding  cone  angle  is  ~20°. 
Figure  5.3  shows  that  the  optical  emission  intensity  fades  away  at  <  =1.18  sec  simultaneously 
at  all  points  along  the  jet. 

Figure  5.4a  shows  the  altitude  distribution  of  the  optical  emission  intensity  at  the  axis  of 
symmetry  (r  =  0)  for  selected  bands  at  t  =1.04  sec.  The  first  and  second  positive  bands 
of  N2  dominate  and  have  similar  intensities.  Figure  5.4b  shows  the  red  and  blue  video 
responses  at  /=1 .04  sec. 

Figure  5.5  shows  the  spatially  integrated  value  of  the  optical  emission  intensity  as  a 
function  of  time.  The  intensity  peaks  ~200  ms  after  the  start  of  the  jet  and  disappears  over 
~1 50-200  ms,  defined  primarily  by  the  chosen  value  of  r^.  Additional  calculations  (not 
shown)  indicate  that  the  jet  can  last  as  long  as  the  thundercloud  charge  can  support  the  flow 
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Fig.  5.3.  Optical  emission  intensity  associated  with  a  blue  jet  A  crosscctional  view 
of  the  distribution  of  optical  emission  intensity  of  the  first  positive  band  of  Nt  at  selected 
instants  of  time.  Right  hand  panel  shows  altitude  scans  at  r= 0. 
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Fig.  5.4.  Optical  emission  intensity  for  different  optical  bands,  (a)  Altitude  scans  at 
(r= 0)  of  optical  emission  intensity  I  corresponding  to  different  optical  bands;  (b)  Expected 
red  and  blue  video  response  corresponding  to  optical  emissions  shown  in  (a). 


of  current  from  the  cloud  to  the  ionosphere.  In  this  case,  the  jet  is  stationary  and  has  the 
shape  of  a  significantly  expanded  upward  luminous  column  which  electrically  connects  the 
thundercloud  and  the  higher  conducting  regions  of  the  atmosphere. 
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Fig.  5.5.  Spatially  integrated  optical  emission  intensity  as  a  function  of  time. 


5.5  COMPARISON  WITH  VIDEO  OBSERVATIONS  OF  BLUE  JETS 


5.5.1  Geometrical  Shape  and  Velocity  of  Blue  Jets 

The  dynamics  of  the  source  charge  Q  may  play  a  dominant  role  in  defining  the  radial  extent 
rs  and  velocity  vs  of  jets.  Equation  (5.1)  and  the  continuity  of  the  streamer  current  Is 
indicate  that  the  velocity  vs  can  stay  constant  with  altitude  if  rs~iV-1  and  NC~N2.  The 
initial  radii  of  jets  in  such  a  case  should  he  very  small  (100-500  m)  to  produce  the  observed 
collimated  (with  cone  angle  <  30°)  beam-like  shape  of  the  lower  parts  of  jets  with  altitude 
extent  of  several  atmospheric  scales.  Such  a  result  may  be  consistent  with  the  fan-like  upper 
parts  of  blue  jets  [Wescott  et  al,,  1995a]  (see  also  lower  image  in  Figure  1.12);  however,  the 
optical  emission  intensity  (7  ~  Ne)  in  this  case  dramatically  decreases  with  altitude  with  a 
scale  height  of  rs*/  3  km. 

5.5.2  Brightness  and  Color  of  Blue  Jets 

Figure  5.4b  indicates  that  the  optical  emission  intensity  changes  by  a  factor  of  ~  1 02  between 
the  lower  and  upper  ends  of  the  jet,  in  good  agreement  with  observations  [  Wescott  et  al., 
1995a],  However,  the  predicted  absolute  intensities  ~5x104  kR  (Figure  5.4b)  are  ~l-2 
orders  of  magnitude  above  those  observed  (~10^  kR  [ E .  Wescott,  private  communication]). 
In  this  regard,  we  note  that  the  excitation  rates  of  the  optical  emissions  depend  sensitively  on 
the  electric  field  (e.g.,  £p  and  may  change  several  orders  of  magnitude  in  response  to  small 
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variations  in  the  electric  field.  In  addition,  the  uncertain  role  of  detachment  processes  (as 
discussed  in  previous  sections)  may  have  lead  to  an  overestimation  of  the  optical  emission 
intensities. 

Concerning  the  relation  between  red  and  blue  colors  we  note  that  the  distribution  as 
a  function  of  wavelength  of  the  first  and  second  positive  bands  of  N2  may  significantly 
differ  at  low  altitudes  from  that  observed  in  aurora  at  ionospheric  altitudes  and  used  in 
our  calculations  of  the  red/blue  video  response  for  measurements  of  sprites  in  Figure  5.4b. 
For  example,  the  Av=3,4  sequences  of  the  first  positive  band  of  N2,  producing  red  visible 
emissions,  may  be  significantly  degraded  in  comparison  with  that  shown  in  Figure  5.4b 
due  to  a  reduction  in  the  populations  of  high  vibrational  levels  with  increasing  pressure 
[e.g.,  Morrill  et  ol.,  1 988].  Since  the  first  and  second  positive  bands  are  the  only  significant 
emissions  (Figure  5.4a)  the  second  positive  band  is  the  most  probable  source  of  the  blue 
color  of  jets. 

5.6  BLUE  JETS  AND  RUNAWAY  ELECTRONS 

In  the  context  of  blue  jets,  runaway  electrons  may  possibly  be  involved  in  producing  the 
initial  ionization  in  the  streamer  head.  As  was  previously  discussed,  the  runaway  breakdown 
mechanism  appears  attractive  since  it  has  a  threshold  field  Et  which  is  ~ 1 0  times  lower 
than  the  conventional  breakdown  field  E^.  However,  in  the  lower  atmosphere,  the  runaway 
electron  process  starts  from  very  low  ambient  production  rates  of  energetic  electrons  pro¬ 
duced  by  cosmic  rays  ( S0 ),  requiring  many  e-folding  distances  before  observable  effects 
can  be  produced  [Bell  et  til.,  1 995].  As  a  result,  the  electric  field  should  in  fact  significantly 
exceed  (>5  times)  the  threshold  field  Et  in  order  to  reach  a  regime  of  significant  ionization 
and  optical  output. 

To  elaborate  further,  consider  the  case  E  =  lEt  and  two  different  altitudes  of  30  and  40 
km.  Runaway  electron  e-folding  distances  for  these  altitudes  are  If  ~50  m  and  ~200  m, 
respectively  [Bell  et  ol.,  1995].  To  support  a  streamer  with  t5=0.05  s  and  r.s  =2  km,  the 


102 


BLUE  JETS  PRODUCED  BY  PRE-DISCHARGE  THUNDERCLOUD  HELDS 


runaway  process  should  produce  sufficient  secondary  electrons  (e.g.,  Ns  ~  101  cm-3;  see 
Figure  5.2)  on  a  time  scale  much  shorter  than  rs,  say  At  ~1  ms.  Using  Ns  ~  cN rN cr0At, 
where  c  is  speed  of  light,  <r0  =  2.3  x  10-18  cm2  [Bell  et  al.,  1995]  the  required  number  of 
runaway  electrons  can  be  estimated  as  Nr  ~  3.9  x  10-s  and  1 .8  x  10-4  cm-3  at  altitudes 
30  and  40  km,  respectively.  The  number  of  required  e-foldings  nj  can  be  estimated  from 
Nji~S0lf/c  e1lf,  and  is  equal  to  20.1  and  21.6  for  30  and  40  km,  respectively,  where 
5„=4.4  x  10-7  cm-3  s-1,  at  30  km  and  So=10-7  cm-3  s-1  at  40  km  [Bell  et  al.,  1995]. 
Thus,  runaway  electrons  should  move  at  least  ~1  km  at  30  km  and  ~4  km  at  40  km  altitude 
to  produce  sufficient  density  Ns  of  secondaries.  These  values  are  lower  than  observed  scales 
of  jets,  so  that  the  runaway  mechanism  may  well  contribute  to  their  production.  In  other 
words,  secondary  electrons  Ns  may  play  a  major  role  in  establishing  the  current  and  space 
charge  balance,  as  well  as  the  general  shape  of  the  streamer.  Since  the  thermalization  time  of 
secondaries  is  expected  to  be  very  fast  at  ~30  km  altitudes,  the  streamer  development  in  case 
of  ionization  by  runaway  electrons  may  be  similar  to  that  caused  by  the  quasi-electrostatic 
heating,  including  the  propagation  velocity  us~  100  km/s.  However,  the  optical  spectra  of 
the  resultant  emissions  may  differ  since  the  secondary  electrons  produced  by  the  runaway 
process  would  in  general  be  more  energetic  [Bell  et  al.,  1995]. 


6 

Summary  and  Suggestions  for  Future 
Research 


6.1  SUMMARY  OF  RESULTS 

A  new  form  of  electrodynamic  coupling  between  thunderstorms  and  the  mesosphere  and 
lower  ionosphere  has  been  identified,  involving  heating,  ionization  and  optical  emission 
effects  produced  by  quasi-electrostatic  thundercloud  fields,  and  accounting  for  the  basic 
features  of  mesospheric  optical  emissions  called  sprites  and  jets. 

Sprites  are  proposed  to  be  a  consequence  of  large  electric  field  transients  capable  of 
causing  electron  heating,  breakdown  ionization  and  excitation  of  optical  emissions  at  meso¬ 
spheric  altitudes  following  the  removal  of  a  thundercloud  charge  by  a  cloud-to-ground 
discharge.  These  transients  can  be  produced  by  the  combined  action  of  induced  spatially 
distributed  space  charge  which  exists  in  the  region  between  the  thundercloud  tops  and  the 
lower  ionosphere  and  the  charge  left  in  the  cloud  after  the  lightning  discharge.  Results 
obtained  with  a  two  dimensional  and  selfconsistent  model  of  the  nonlinear  upward  coupling 
of  quasi-electrostatic  thundercloud  fields  are  consistent  with  the  observed  spectra,  temporal 
duration,  and  spatial  extent  of  sprites. 

The  heating  of  the  mesosphere/lower  ionosphere  by  quasi-electrostatic  fields  can  lead  to 
significant  modification  of  the  conductivity  of  these  regions.  Model  calculations  predict 
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significant  (several  orders  of  magnitude)  modification  of  the  lower  ionospheric  conductivity 
in  a  region  with  50-60  km  transverse  extent.  This  modification  can  be  in  the  form  of  deple¬ 
tions  of  electron  density  due  to  dissociative  attachment  to  O2  molecules  or  enhancements  of 
electron  density  due  to  breakdown  ionization  of  O2  and  N2.  Ionization  changes  produced 
by  consecutive  lightning  discharges  can  accumulate,  leading  to  the  build  up  of  ionization 
“bubbles”  or  attachment  “holes”  above  thunderstorms. 

Results  of  the  model  calculations  corresponding  to  fast  (~  1  ms)  removal  of  100-300  C 
of  thundercloud  charge  from  10  km  altitudes  are  in  good  agreement  with  observations  of 
the  upper  part  (“head”  and  “hair”  [ Sentman  et  al.,  1 995])  of  the  sprites.  The  typical  region 
of  brightest  optical  emission  has  horizontal  and  vertical  dimensions  ~10  km,  centered  at 
altitudes  ~70  km  and  is  interpreted  as  the  head  of  sprite.  Our  model  results  also  indicate 
the  formation  of  low  intensity  glow  (“hair”)  above  this  region  due  to  the  excitation  of 
optical  emissions  at  altitudes  ~  85  km  during  the  ~500  /is  at  the  beginning  of  the  lightning 
discharge.  Comparisons  of  the  optical  emission  intensity  of  the  1st  and  2nd  positive  bands 
of  N2,  Meinel  and  1st  negative  bands  of  Nt,  and  1st  negative  band  of  O}  show  that  the 
1st  positive  band  of  N2  is  the  dominant  emission  under  conditions  of  heating  by  quasi¬ 
electrostatic  thundercloud  fields. 

The  dependencies  of  quasi-electrostatic  heating  upon  the  altitude  profile  of  ion  conduc¬ 
tivity,  number  density  of  electrons  and  duration  of  the  lightning  discharge  indicate  that  the 
level  of  optical  emissions  is  defined  mainly  by  the  relation  between  the  lightning  discharge 
duration  and  the  conductivity  of  the  upper  atmosphere  (i.e.,  relaxation  time  of  the  electric 
field  in  the  conducting  medium).  For  low  ambient  conductivities  the  lightning  discharge 
duration  can  be  significantly  longer  than  1  ms  and  still  produce  similar  levels  of  optical 
emissions.  The  maximum  possible  intensity  of  optical  emissions  is  determined  mainly  by 
the  value  of  the  removed  thundercloud  charge  and  its  altitude. 

Comparison  of  the  model  results  for  different  types  of  lightning  discharges  (+CG,  -CG, 
IC)  indicate  that  positive  cloud-to-ground  discharges  lead  to  the  largest  electric  fields  and 
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optical  emissions  at  ionospheric  altitudes,  but  only  because  they  are  typically  associated 
with  the  removal  of  larger  amounts  of  charge  from  higher  altitudes.  Studies  of  different 
geometries  of  thundercloud  charges  suggest  that  spherically  distributed  charges  produce 
more  intense  electric  fields  and  optical  emissions  than  horizontally  spread  charges. 

The  quasi-electrostatic  model  of  lower  ionospheric  heating  and  ionization  indicates  that 
most  of  the  observed  features  of  sprites  with  vertically  oriented  columnar  structure  can  be 
explained  in  terms  of  the  formation  and  self-driven  propagation  of  streamer  type  channels  of 
breakdown  ionization.  Optical  emissions  of  the  dominant  first  positive  and  second  positive 
bands  of  molecular  nitrogen  may  form  carrot-like  vertical  structures  with  typical  transverse 
dimensions  of  ~5-10  km  which  can  span  the  altitude  range  from  ~80  km  to  well  below  ~50 
km,  with  detailed  structure  depending  primarily  on  the  time  history  of  charge  removal  from 
the  thundercloud  by  the  cloud-to-ground  lightning  discharge.  The  appearance  of  optical 
emissions  associated  with  sprites  can  be  delayed  in  time  (~  1  -20  ms)  with  respect  to  the 
causative  cloud-to-ground  discharge. 

Our  results  indicate  that  blue  jets  are  streamer  type  processes  occurring  on  atmospheric 
spatial  scales,  with  most  features  of  jets  as  observed  in  video  [Wescott  et  al„  1995a]  being 
reproduced  in  the  context  of  the  quasi-electrostatic  model,  including  the  velocity  of  up¬ 
ward  propagation  of  ~100  km/s,  the  general  shape  of  jets  as  upward  expanding  beams  of 
luminosity  with  cone  angles  <  30°,  and  the  stopping  altitudes  (~50  km)  for  jets.  Model 
results  are  in  agreement  with  video  observations  showing  that  the  optical  brightness  drops 
~2  orders  of  magnitude  between  the  lower  part  and  the  tip  of  the  jet,  and  the  fact  that  the 
luminosity  fades  away  simultaneously  everywhere  along  the  cone  of  the  jet.  The  blue  color 
of  jets  as  observed  in  video  is  interpreted  as  the  emission  from  the  second  positive  band  of 
molecular  nitrogen  excited  upon  electron  impact.  Jets  are  produced  during  the  accumulation 
phase  of  the  thundercloud  charge,  so  that  it  is  not  necessary  for  large  amounts  of  charge  to 
be  removed.  Thus,  jets  are  not  necessarily  associated  with  lightning  discharges  and  may 
appear  only  in  relatively  rare  cases  of  large  (~300-400C)  thundercloud  charge  accumula¬ 
tion  at  high  (~20  km)  altitudes.  It  is  conceivable  that  a  horizontally  extended  thundercloud 
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charge  distribution,  consisting  of  many  charge  centers  located  at  lower  altitudes  (<  20km), 
may  produce  the  electric  field  equivalent  to  that  of  simple  monopole  of  ~300-400C  at  ~20 
km  altitude. 


6.2  SUGGESTIONS  FOR  FUTURE  RESEARCH 

The  existing  quasi-electrostatic  heating  model  should  be  extended  to  treat  two  fluid  equa¬ 
tions  for  the  electron  and  ion  components  to  allow  a  more  general  treatment  of  streamer 
processes  associated  with  sprites  and  jets.  Important  questions  not  yet  addressed  concern  the 
chemical  and  thermal  effects  produced  by  quasi-electrostatic  fields  and  runaway  electrons 
associated  with  sprites  and  jets  in  the  mesosphere  and  the  selfconsistent  modification  of 
the  quasi-electrostatic  fields  by  runaway  electrons.  The  existing  model  should  be  modified 
to  include  the  two  dimensional  dynamics  of  runaway  electrons  and  their  secondary  elec¬ 
tron  products  (see  Sections  4.8,  5.6).  Relative  excitation  intensities,  spatial  distribution, 
and  duration  of  various  emissions  (ranging  from  infrared  to  ultraviolet)  stimulated  simul¬ 
taneously  by  quasi-electrostatic  heating  and  runaway  processes  need  to  be  evaluated,  and 
criteria  should  be  formulated  for  comparisons  with  future  spectroscopic  measurements  of 
sprites  (see  discussion  in  section  4.8).  It  is  important  to  relate  these  characteristics  to  the 
lightning  characteristics  and  ambient  conditions  in  the  atmosphere,  mesosphere  and  lower 
ionosphere. 

Sprites  and  jets  occur  under  conditions  of  exceptionally  intense  electrical  and  lightning 
activity  which  is  characteristic  of  large  mesoscale  convective  systems  [e.g.,  Lyons,  1996]. 
Recent  theoretical  and  experimental  findings  [Inan  et  al.,  1996d]  suggest  that  even  rela¬ 
tively  small  quasi-electrostatic  fields  above  thunderclouds  are  capable  of  maintaining  the 
ionospheric  electrons  at  a  persistently  heated  level  well  above  their  ambient  thermal  en¬ 
ergy.  Changes  in  the  thundercloud  charge  distribution  (e.g.,  in  lightning  discharges)  lead 
to  heating/cooling  above/below  this  quiescent  level,  and  can  be  experimentally  detected 
as  changes  in  the  amplitudes  and  phases  of  subionospherically  propagating  VLF  signals 
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[hum  et  ai,  1996a,d].  These  changes  occur  simultaneously  with  lightning  discharges,  and 
are  referred  to  as  early/fast  VLF  events.  These  rapid  ionospheric  conductivity  changes 
may  thus  be  evidence  for  thermal  and  aeronomic  coupling  between  the  troposphere  and 
the  mesosphere/lower  ionosphere  regions,  with  potentially  more  important  (longer  term) 
implications  than  the  spectacular  optical  displays  (i.e.,  sprites  and/or  elves)  observed  during 
the  brief  transients  immediately  after  lightning  discharges.  In  view  of  these  findings,  an  im¬ 
portant  next  step  in  future  research  should  be  to  determine  the  spatial  structure  and  temporal 
dynamics  of  long  term  energy  transfer  from  thunderstorm  generated  electric  fields  to  the 
mesosphere/lower  ionosphere  and  to  determine  the  consequences  of  this  deposited  energy 
on  the  long  term  chemical  balance  of  the  D-region  of  the  ionosphere.  Important  topics 
to  investigate  here  would  be:  (i)  the  spatial  distributions  of  the  temperature  and  number 
density  of  the  different  kinds  of  charged  particles  and  the  number  density  of  rotationally, 
vibrationally  and  optically  excited  molecules,  and  (ii)  the  importance  of  these  effects  in  the 
global  thermal  and  chemical  balance  of  the  upper  atmosphere. 

It  is  also  important  to  investigate  whether  the  variety  of  observed  shapes  of  sprites  can  all  be 
explained  in  terms  of  the  quasi-electrostatic  and  runaway  mechanisms.  Lower  ionospheric 
conductivity  inhomogeneities  (e.g.,  created  by  particle  precipitation)  may  play  a  key  role 
in  formation  of  sprite  structures  (see  section  4.1);  possible  sources  of  the  electron,  ion  or 
neutral  density  inhomogeneities  in  the  mesosphere/lower  ionosphere  above  thunderstorms 
need  to  be  investigated.  The  existing  electrodynamic  models  should  be  modified  to  treat 
typical  variety  of  spatial  configurations  of  such  inhomogeneities.  The  most  recent  theoretical 
developments  [Pasko  et  ai,  1 996d]  suggest  that  heating  of  the  neutral  atmosphere  above 
thunderstorms  (for  which  there  is  some  experimental  evidence  [e.g.,  Taylor,  1979])  may 
play  an  important  role  in  the  modulation  of  the  neutral  atmospheric  density  and  may  lead 
to  the  creation  of  the  observed  shapes  of  sprites.  The  possible  means  of  heating  the  neutral 
atmosphere  above  thunderstorms,  for  example  by  internal  gravity  waves  [Taylor,  1979]  and 
through  collisional  transfer  of  energy  from  heated  ionospheric  electrons  to  neutrals  [Pasko 
et  al.,  1996d]  should  also  be  investigated. 


A 

Comparison  with  Kinetic  Solutions  and 
Swarm  Experiments 


In  this  appendix,  we  evaluate  the  validity  of  parametric  approximations  for  mobility,  ion¬ 
ization,  attachment  and  optical  excitation  rates  (Sections  2.4. 1,  2.4.2,  2.4.3,  2.5)  with  the 
results  of  kinetic  solutions  [Taranenko  et  al.,  1993a,b]  and  with  the  results  of  the  electron 
swarm  experiments  where  available  [e.g.,  Davies,  1983].  For  this  purpose,  we  consider  the 
altitude  profile  (at  r  =  0)  of  the  electric  field  which  approximately  corresponds  to  a  charge 
of  300  C  placed  at  10  km  altitude  km  (M2  model  from  Figure  2.1).  We  then  divide  the 
altitude  range  from  10  to  90  km  in  equidistant  intervals  with  1  km  steps.  At  each  altitude, 
the  value  of  the  electric  field  is  divided  into  50  linear  intervals  so  that  the  first  value  is 
0  and  last  value  is  the  maximum  value  of  the  electric  field  at  this  altitude.  The  resulting 
discretization  of  the  electric  fields  and  altitudes  is  illustrated  in  Figure  A.l.  The  stationary 
electron  distribution  function  f(v)  corresponding  to  the  each  of  the  points  shown  in  Fig¬ 
ure  A.  1  (i.e.,  for  different  electric  fields  and  densities  of  atmospheric  gas)  is  obtained  from 
the  solution  of  the  Boltzman  equation,  taking  into  account  inelastic  collisions  consisting  of 
rotational,  vibrational,  optical,  dissociative,  dissociative  with  attachment  and  ionizational 
losses  [Taranenko  etal.,  1993a].  The  mobility,  ionization,  attachment  and  optical  excitation 
rates  are  calculated  using  the  obtained  distribution  functions  and  known  cross  sections  for 
these  processes  [Taranenko  et  al.,  1 993a,b] .  Results  obtained  with  this  approach  are  di¬ 
rectly  comparable  to  swarm  laboratory  experiments  [e.g.,  Davies,  1983]  in  which  when  the 
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electron  mobility  is  measured  for  different  magnitudes  of  applied  electric  fields  and  ambient 
neutral  gas  densities. 

A.l  ELECTRON  MOBILITY 

The  electron  mobility  can  be  calculated  as 


(A.l) 


where  the  distribution  function  is  assumed  to  be  normalized  so  that  Ne  =  f  f(v)d?v.  Re¬ 
sults  are  shown  in  Figure  A. lb  by  dots.  We  also  show  the  corresponding  values  obtained 
from  experimental  measurements  by  Davies  [1983]  (open  circles)  and  those  obtained  using 
the  analytical  approximation  (2.5)  (dashed  line).  For  lower  values  of  the  electric  fields 
( EN0/N  <  104  V/m)  the  analytic  approximation  (2.5)  is  supported  by  measurements  of 
Hegerberg  and  Reid  [  1980]  (not  shown).  The  normalized  mobility  values  of  (j,eN/N0  =  1.36 
m2/V/s  for  EN0/N  <  1 .62  x  1 0^  V/m  corresponds  to  the  thermal  state  of  the  electrons  (with 
no  heating).  The  value  of  the  mobility  in  this  region  of  electric  fields  is  found  to  be  in  good 
agreement  with  that  derived  from  hydrodynamic  equations  using  effective  collision  frequen¬ 
cies  and  average  fractions  of  energy  transferred  in  one  collision  from  electrons  to  neutrals 
at  an  electron  temperature  Te  =  200  K  [Gurevich,  1978,  pp.  32, 99,  103],  The  results  of  the 
three  different  approaches  (kinetic,  experimental,  analytical)  shown  in  Figure  A.  1  b  agree 
with  each  other  reasonably  well. 

A.2  IONIZATION  AND  ATTACHMENT  COEFFICIENTS 

The  ionization  Uj  and  attachment  va  coefficients  can  be  calculated  as: 


(A.2) 


(A. 3) 


N0  /N  (s  -1 )  vi,aNo  /N  (s  -1 )  Altitude  (km) 
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Fig.  A.l.  Comparison  of  analytical  approximations  with  kinetic  solutions  and  results 
of  swarm  experiments,  (a)  The  discretization  of  altitudes  and  electric  field  values  for  the 
kinetic  solutions;  (b)  The  electron  mobility  corresponding  to  different  models  is  shown  in  the 
same  formal  as  in  Figure  2.5;  (c)  The  ionization  and  attachment  coefficients  corresponding 
to  different  models  are  shown  in  the  same  format  as  in  Figure  2.6;  (d)  The  excitation 
coefficients  of  selected  optical  bands  corresponding  to  different  models  are  shown  in  the 
same  formal  as  in  Figure  2.7;  (e)  The  excitation  coefficient  of  vibrationally  excited  N2 
molecules;  (f)  The  average  energy  of  the  electron  distribution 
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where  oy  and  oa  are  corresponding  cross  sections.  Comparison  of  the  results  of  kinetic  cal¬ 
culations  (A. 2),  (A. 3)  (dots)  and  the  analytical  approximations  (2.6)  (2.7)  (dashed  lines)  for 
v\  and  ua  are  shown  in  the  Figure  A.lc.  For  the  attachment  coefficient  we  also  show  results 
of  experimental  measurements  by  Davies  [1983]  (open  circles).  We  note  once  again  the 
reasonable  agreement  between  the  different  approaches.  For  the  electron  attachment  coef¬ 
ficient  in  the  range  of  electric  fields  EN0/N  >  4  x  106  V/m  we  see  a  significant  difference 
between  the  experimental  data  and  the  results  of  kinetic  solutions  (most  probably  caused  by 
detachment  effects  as  discussed  in  [Davies,  1983]);  however,  in  this  region  of  electric  fields 
the  attachment  coefficient  becomes  much  smaller  than  the  ionization  coefficient  so  that  the 
relative  inaccuracy  of  the  parametric  model  does  not  affect  our  numerical  results. 


A .3  OPTICAL  EMISSIONS 

The  excitation  rate  per  one  electron  of  a  given  energy  state  k  for  a  particular  species  j  has 
the  value: 

J  vf(v)ak(v)d~v,  (A. 4) 

where  rrk  is  the  cross  section  for  excitation  of  state  k.  For  the  sake  of  brevity  in  Figure  A.  Id 
we  show  only  results  corresponding  to  the  excitation  rates  of  the  1  st  positive  band  of  N2 
and  the  Meinel  band  of  N£.  The  comparison  between  kinetic  solutions  (A.4)  (dots)  and  the 
analytical  approximations  (2.10)  (dashed  lines)  demonstrate  that  the  later  can  be  used  with 
negligible  error,  as  we  have  done  in  our  model  calculations. 


A.4  VIBRATIONAL  EXCITATION  OF  MOLECULAR  NITROGEN 


In  Figure  A.le  we  show  the  excitation  rate  (per  one  electron)  of  the  vibrationally  excited 
N2  molecules: 


(A. 5) 
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where  av  is  the  cross  section  corresponding  to  this  process  [Taranenko,  1993].  Although 
not  discussed  in  this  thesis,  production  of  vibrationally  excited  NT2  molecules  may  play  an 
important  role  in  the  chain  of  reactions  leading  to  4.3  //m  infrared  emissions  of  CO2  [e.g., 
Picard  et  al.,  1987].  Results  of  kinetic  solutions  in  this  case  can  also  be  approximated  by  a 
relatively  simple  analytical  formula: 

N  ^ ' 

log  =  2^  a'x'  (A-6) 

?=() 

where  x  =  logC^^1),  and  aQ  =  — 1284.6,  a  1  =  1 154. 1 , 02  =  —408.37,  <23  =  71.513,  «4  = 
—6. 1 882,  05  =  0.2 1158,  which  is  shown  in  Figure  A.  1  e  by  a  dashed  line.  The  approximation 
given  by  (A. 5)  is  valid  for  the  range  of  electric  fields  104  <  ENo/N  <  107V/m. 


A.5  AVERAGE  ELECTRON  ENERGY 

For  illustrative  purposes  in  Figure  A.  If  we  provide  the  average  energy  of  the  electron 
distribution 

IV  =  jj-  f  T^—f(v)d}'v  (A. 7) 

as  a  function  of  the  electric  field.  Calculations  show  that  even  in  cases  of  relatively  strong 
ionization,  as  shown  in  Figure  3.4  for  the  removed  charge  250  C,  the  peak  average  energy 
of  the  electron  distribution  does  not  exceed  ~5  eV,  indicating  that  most  of  the  ionization  is 
produced  by  the  relatively  small  number  of  electrons  in  the  tail  of  the  electron  distribution 
function  which  appear  in  the  region  of  energies  of  ~15  eV  (corresponding  to  ionization 
threshold  of  N2  and  O2). 


B 

Comparison  with  Auroral  Optical 
Emissions 


It  is  interesting  to  compare  the  spectral  content  of  optical  emissions  associated  with  aurorae 
and  those  produced  at  the  same  altitudes  by  electrons  heated  by  electric  fields  resulting  from 
intensive  cloud  to  ground  lightning  discharges. 

The  typical  omnidirectional  flux  of  auroral  electrons  (in  erg/cm2/s/eV)  can  be  approxi¬ 
mated  by  the  expression  [Meier  et  ai,  1989]: 

€'  e 

<t>(£)  =  +  0.4<t>max— e-^.  (B.  1 ) 

where 

(B-2) 

and  <t>max  is  the  maximum  value  of  E  is  energy  in  keV,  and  b  =  O.Sif'  if  S'  <  0.5  keV. 
or  b  =  0.  l<f'+0.35  if  S'  >  0.5  keV,  E  is  the  characteristic  energy,  and  P  is  the  energy  flux. 
The  excitation  rate  of  the  Arth  excitation  process  involving  the  jth  species  can  be  calculated 
as: 

v'k  =  Nj  J  <P(£)nk.(£)d£ .  (B.3) 

where  <r*  is  corresponds  to  fcth  excitation  process  cross  section. 
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(S/A9/SW3/9)  <t> 
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Results  of  calculations  of  14-  using  (B.3)  for  selected  optical  bands  and  range  of  the 
characteristic  energies  S'  0.1-10  keV  for  P  =30  ergs/cm2/s  are  shown  in  Figure  B.lb 
along  with  the  electron  flux  for  three  selected  values  of  S'  (Figure  B.l  a).  The  1st  negative 
and  Meinel  bands  of  N£  are  clearly  dominant,  and  the  auroral  emissions  exhibit  a  rich 
spectral  complexity  in  comparison  with  those  produced  by  quasi-electrostatic  fields  (e.g., 
Figure  3.9b). 


Fig.  B.l.  Comparison  with  auroral  optical  emissions,  (a)  Omnidirectional  flux  of  auroral 
electrons;  (b)  Excitation  rates  of  selected  optical  bands  corresponding  to  (a)  ;  (c)  Typical 
distribution  function  of  heated  electrons  and  cross  sections  of  the  optical  excitation  of  1st 
positive  band  of  N2  and  Meinel  band  of  N,. 


Figure  B.  1  b  shows  the  excitation  rates  rather  than  the  emission  rates.  At  altitudes  of  ~85 
km  used  for  our  calculations,  the  Meinel  Nt  band  is  quenched,  since  the  quenching  height 
for  this  band  is  ~90  km.  The  resulting  emission  rate  is  thus  ~2-3  times  lower  than  the 
excitation  rate  plotted  in  Figure  B.  lb.  The  rest  of  the  bands  shown  have  quenching  heights 
much  lower  than  85  km  so  that  the  corresponding  excitation  rates  can  be  interpreted  as 
emission  rates. 

The  primary  reason  for  the  spectral  differences  between  in  the  auroral  precipitation  and 
quasi-electrostatic  heating  cases  is  the  relative  structure  of  the  electron  distribution  function 
as  a  function  of  energy  and  the  corresponding  excitation  cross  sections.  In  the  case  of 
electron  heating  by  quasi-static  electric  fields,  the  distribution  function  sharply  decreases 
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with  increasing  energy  in  the  region  of  electron  energies  of  10-1 5  eV  (due  to  the  effective  loss 
of  electron  energy  as  a  result  of  inelastic  collisions  [Taranenko  et  al.,  1993])  as  illustrated 
in  Figure  B.lc,  for  h  =85  km  and  an  applied  electric  field  of  45  V/m.  For  illustration,  the 
distribution  of  excitation  cross  sections  for  N2  1  st  positive  and  N2  Meinel  bands  [Taranenko, 
1 993]  are  also  shown  in  Figure  B.lc.  In  the  auroral  case,  the  distribution  function  has  a  wide 
extent  and  covers  a  broad  energy  range  up  to  104-  105  eV  (Figure  B.la).  In  simpler  terms, 
the  secondary  electron  distribution  in  aurorae  are  produced  by  impact  of  higher  energy 
precipitating  electrons  and  thus  has  a  distribution  with  larger  numbers  of  electrons  at  higher 
energies. 


Bibliography 


Armstrong,  W.  C.,  Recent  advances  from  studies  of  the  Trimpi  effect,  Antarctic  J.,  IK  28 1 , 
1983. 

Baginski,  M.  E.,  L.  C.  Hale,  and  J.  J.  Olivero,  Lightning-related  fields  in  the  ionosphere, 
Geophys.  Res.  Lett.,  15,  764,  1988. 

Baginski,  M.  E.,  and  A.  S.  Hodel,  A  case  study  comparing  the  lossy  wave  equation  to  the 
continuity  equation  in  modeling  late-time  fields  associated  with  lightning.  Applied 
Computational  Electromagnetics  Society  Journal,  9,  98,  1994. 

Bell,  T.  F.,  V.  P.  Pasko,  and  U.  S.  Inan,  Runaway  electrons  as  a  source  of  Red  Sprites  in  the 
mesosphere.  Geophys.  Res.  Lett.,  22,  2127,  1995. 

Boccippio,  D.J.,  E.  R.  Williams,  S.  J.  Heckman,  W.  A.  Lyons,  I.  T.  Baker,  and  R.  Boldi, 
Sprites,  ELF  transients,  and  positive  ground  strokes.  Science,  269,  1088,  1995. 

Boeck,  W.  L.,  O.  H.  Vaughan,  Jr.,  R.  Blakeslee,  B.  Vonnegut,  and  M.  Brook,  Lightning 
induced  brightening  in  the  airglow  layer,  Geophys.  Res.  Lett.,  19,  99,  1992. 

Boeck,  W.  L.,  O.  H.  Vaughan,  R.  J.  Blakeslee,  B.  Vonnegut,  M.  Brook,  and  J.  McKune. 
Observations  of  lightning  in  the  stratosphere,  J.  Geophys.  Res.,  100,  1465,  1995. 

Braham,  Jr.,  R.R.  The  water  and  energy  budgets  of  the  thunderstorm  and  their  relation  to 
thunderstorm  development,  J.  Meteor.,  9,  227 ,  1952. 

Brook,  M.,  M.  Nakano,  P.  Krehbiel,  and  T.  Takeuti,  The  electrical  structure  of  the  Hokuriku 
winter  thunderstorms,  J.  Geophys.  Res.,  87,  1207,  1982. 


119 


120 


BIBLIOGRAPHY 


Buneman,  O.,  A  compact  non-iterative  Poisson-solver,  SUIPR  report  No.  294,  Stanford 
University,  1969. 

Chamberlain,  J.  W.,  Theory  of  planetary  atmospheres.  Academic  Press,  New  York,  1978. 

Davies,  D.K.,  Measurements  of  swarm  parameters  in  dry  air,  Theoretical  Notes,  Note  346, 
Westinghouse  R&D  Center,  Pittsburg,  May,  1983. 

Dejnakarintra,  M.,  and  C.  G.  Park,  Lightning-induced  electric  fields  in  the  ionosphere,  J. 
Geophys.  Res.,  79,  1903,  1974. 

Dhali.  S.K.,  and  P.  F.  Williams,  Two-dimensional  studies  of  streamers  in  gases,  J.  Appl. 
Phys.,  62,  4696,  1987. 

Dowden,  R.  L.,  C.  D.  C.  Adams,  J.  B.  Brundell  and  P.  E.  Dowden,  Rapid  onset,  rapid  decay 
(RORD),  phase  and  amplitude  pertubations  of  VLF  subionospheric  transmissions, 
J.  Atmos.  Terr.  Phys.,  56,  1513,  1994. 

Fishman,  G.  J.,  P.  N.  Bhat,  R.  Mallozzi,  J.  M.  Horack,  T.  Koshut,  C.  Kouveliotou.  G.  N. 
Pendleton,  C.  A.  Meegan,  R.  B.  Wilson,  W.  S.  Paciesas,  S.  J.  Goodman,  and  H.  J. 
Christian,  Discovery  of  intense  gamma-ray  flashes  of  atmospheric  origin,  Science, 
264,  1313,  1994. 

Franz,  R.  C.  et  al.,  Television  image  of  a  large  upward  electric  discharge  above  a  thunderstorm 
system,  Science,  249,  48,  1990. 

Fukunishi,  H.,  Y.  Takahashi,  M.  Kubota,  K.  Sakanoi,  U.  S.  Inan,  and  W.  A.  Lyons,  Lightning- 
induced  transient  luminous  events  in  the  lower  ionosphere,  Geophys.  Res.  Lett.,  in 
press,  1996. 

Glukhov,  V.  S.,  V.  P.  Pasko,  and  U.  S.  Inan,  Relaxation  of  transient  lower  ionospheric 
disturbances  caused  by  lightning-whistler-induced  electron  precipitation  bursts,  J. 
Geophys.  Res.,  97,  16971,  1992. 

Grange,  F.,  N.  Soulem,  J.  F.  Loiseau,  and  N.  Spyrou,  Numerical  and  experimental  determi¬ 
nation  of  ionizing  front  velocity  in  a  DC  point-to-plane  corona  discharge,  J.  Phys. 
D:  Appl.  Phys.,  28,  1619,  1995. 

Greifinger,  C.,  and  P.  Greifinger,  Transient  ULF  electric  and  magnetic  fields  following  a 
lightning  discharge,  J.  Geophys.  Res.,  81,  2237,  1976. 


BIBLIOGRAPHY 


121 


Griffiths,  R.  F.,  and  C.  T.  Phelps,  A  model  of  lightning  initiation  arising  from  positive  corona 
streamer  development,  J.  Geophys.  Res.,  31,  3671,  1976. 

Grigor’ev,  A.  I.,  and  S.  O.  Shiryaeva,  Mechanism  for  development  of  the  stepped  leader  and 
intracloud  branching  of  streak  lightning,  Sov.  Phys.  Tech.  Phys.,  34,  502,  1989. 

Gurevich,  A.  V.,  Nonlinear  phenomena  in  the  ionosphere,  Berlin:  Spriger-Verlag,  1978. 

Gurevich,  A.  V.,  G.  M.  Milikh,  and  R.  A.  Roussel-Dupre,  Runaway  electron  mechanism  of 
air  breakdown  and  preconditioning  during  a  thunderstorm,  Phys.  Lett.  A,  165,  463, 
1992. 

Gurevich,  A.  V.,  G.  M.  Milikh,  and  R.  A.  Roussel-Dupre,  Nonuniform  runaway  air-  break¬ 
down,  Phys.  Lett.  A,  187,  197,  1994. 

Hale,  L.  C.,  and  M.  E.  Baginski.  Current  to  the  ionosphere  following  a  lightning  stroke, 
Nature,  329,  814,  1987. 

Hale,  L.  C.,  Coupling  of  ELF/ULF  energy  from  lightning  and  MeV  particles  to  the  middle 
atmosphere,  ionosphere,  and  global  circuit,  J.  Geophys.  Res.,  99,  21089,  1994. 

Hampton,  D.  L.,  M.  J.  Heavner,  E.  M.  Wescott,  and  D.  D.  Sentman,  Optical  spectral  char¬ 
acteristics  of  sprites,  Geophys.  Res.  Lett.,  22,  89,1996. 

Hegerberg,  R.,  and  I.  D.  Reid,  Electron  drift  velocities  in  air,  Aust.  J.  Phys.,  33,  221,  1980. 

Holden.  D.  N.,  C.  P.  Munson,  and  J.  C.  Devenport,  Satellite  observations  of  transionospheric 
pulse  pairs,  Geophys.  Res.  Lett.,  22,  889,  1995. 

Holzer,  R.  E.,  D.  S.  Saxon,  Distribution  of  electrical  conduction  currents  in  the  vicinity  of 
thunderstorms,  J.  Geophys.  Res.,  57,  207,  1952. 

Holzworth,  R.  H.,  M.  C.  Kelley,  C.  L.  Siefring,  L.  C.  Hale  and  J.  T.  Mitchell,  Electrical 
measurements  in  the  atmosphere  and  the  ionosphere  over  an  active  thunderstorm.  2. 
Direct  current  electric  fields  and  conductivity,  J.  Geophys.  Res.,  90,  9824,  1985. 

Idone,  V.P.,  A.B.  Saljoughy,  R.  W.  Henderson,  P.  K.  Moore,  and  R.  B.  Pyle,  A  reexamina¬ 
tion  of  the  peak  current  calibration  of  the  National  Lightning  Detection  Network, 
Geophys.  Res.  Lett.,  98,  18323,  1993. 

Illingworth,  A.  J.,  Electric  field  recovery  after  lightning  as  the  response  of  the  conducting 
atmosphere  to  a  field  change,  Quart.  J.  R.  Met.  Soc.,  98,  604,  1972. 


122 


BIBLIOGRAPHY 


Inan,  U.  S.,  D.  C.  Shafer,  W.  Y.  Yip,  and  R.  E.  Orville,  Subionospheric  VLF  signatures  of 
nighttime  D-region  perturbations  in  the  vicinity  of  lightning  discharges,  J.  Geophys. 
Res.,  93,  11455,  1988. 

Inan,  U.  S.,  T.  F.  Bell,  and  J.  V.  Rodriguez,  Heating  and  ionization  of  the  lower  ionosphere 
by  lightning,  Geophys.  Res.  Lett.,  18,  705,  1991. 

Inan,  U.  S.,  J.  V.  Rodriguez,  and  V.  R  Idone,  VLF  signatures  of  lightning-induced  heating 
and  ionization  of  the  nighttime  D-region,  Geophys.  Res.  Lett.,  20,  2355,  1993. 

Inan.  U.  S.,  T.  F.  Bell,  V.  P.  Pasko,  D.  D.  Sentman,  and  E.  M.  Wescott,  and  W.  A.  Lyons, 
VLF  signatures  of  ionospheric  disturbances  associated  with  sprites,  Geophys.  Res. 
Lett.,  22,  3461,  1995. 

Inan,  U.  S.,  A.  Slingeland,  V.  P.  Pasko,  and  J.  Rodriguez,  VLF  signatures  of  mesospheric/ 
lower  ionospheric  response  to  lightning  discharges,  J.  Geophys.  Res.,  101,  5219, 
1996a. 

Inan,  U.  S.,  W.  A.  Sampson,  and  Y.  N.  Taranenko,  Space-time  structure  of  lower  ionospheric 
optical  flashes  and  ionization  changes  produced  by  lightning  EMP,  Geophys.  Res. 
Lett.,  23,  133,  1996b. 

Inan,  U.  S.,  S.  C.  Reising,  G.  J.  Fishman,  and  J.  M.  Horack,  On  the  association  of  terrestrial 
gamma-ray  bursts  with  lightning  discharges  and  sprites,  Geophys.  Res.  Lett.,  23, 
1017,  1996c. 

Inan,  U.  S„  V.  P.  Pasko,  and  T.  F.  Bell,  ’Early/Fast  VLF  events  as  evidence  of  sustained 
heating  of  the  ionosphere  above  thunderclouds,  Geophys.  Res.  Lett. ,23,  1067,  1996d. 

Lagarkov,  A.  N.,  and  I.  M.  Rutkevich,  Ionization  waves  in  electrical  breakdown  in  gases, 
Springer- Verlag,  New  York,  1994. 

Lowke,  J .  J .,  Theory  of  electrical  breakdown  in  air  -  the  role  of  metastable  oxygen  molecules, 
J.  Phys.  D:  Appl.  Phys.,  25,  202,  1992. 

Lyons,  W.  A.,  Characteristics  of  luminous  structures  in  the  stratosphere  above  thunderstorms 
as  imaged  by  low-light  video,  Geophys.  Res.  Lett.,  21,  875,  1994. 

Lyons,  W.  A.,  Low-light  video  observations  of  frequent  luminous  structures  in  the  strato¬ 
sphere  above  thunderstorms,  Monthly  Weather  Review,  122,  1940,  1995. 


BIBLIOGRAPHY 


123 


Lyons,  W.  A.,  Sprite  observations  above  the  U.S.  high  plains  in  relation  to  their  parent 
thunderstorm  systems,  J.  Geophys.  Res.,  in  review,  1996. 

Malan,  D.  J.,  Physics  of  lightning.  The  English  Universities  Press  Ltd,  London,  1963. 

Marshall,  T.  C.,  and  W.  D.  Rust,  Two  types  of  vertical  electrical  structures  in  stratiform 
precipitation  regions  of  mesoscale  convective  systems,  Bull.  Amer.  Meteor.  Soc.,  74, 
2159,  1993. 

Marshall,  T.  C.,  M.  Stolzenburg,  W.  D.  Rust,  Electric  field  measurements  above  mesoscale 
convective  systems,  J.  Geophys.  Res.,  1 01,  6979,  1996. 

McCarthy,  M.  P.,  and  G.  K.  Parks,  On  the  modulation  of  X  Ray  fluxes  in  thunderstorms,  J. 
Geophys.  Res.,  97,  5857,  1992. 

Mende,  S.  B.,  R.  L.  Rairden,  and  G.  R.  Swenson,  and  W.  A.  Lyons,  Sprite  spectra:  N2  1  PG 
band  identification,  Geophys.  Res.  Lett.,  22,  2633,  1995. 

Meier,  R.  R.,  D.  J.  Strickland,  J.  H.  Hecht,  and  A.  B.  Christensen,  Deducing  composition 
and  incident  electron  spectra  from  ground-based  auroral  optical  measurements:  A 
study  of  auroral  red  line  processes,  J.  Geophys.  Res.,  94,  13541,  1989. 

Morrill,  J.,  B.  A.  Garragher,  and  W.  Benesh,  Population  development  of  auroral  molecular 
nitrogen  species  in  a  pulsed  electric  discharge,  J.  Geophys.  Res.,  93,  963,  1988. 

Milikh,  G.  M.,  K.  Papadopoulos,  C.  L.  Chang,  On  the  physics  of  high  altitude  lightning, 
Geophys.  Res.  Lett.,  22,  85,  1995. 

Ogawa,  T.,  Fair-weather  electricity,  J.  Geophys.  Res.,  90,  5951,  1985. 

Orville,  R.  E.,  Calibration  of  a  magnetic  direction  finding  network  using  measured  triggered 
lightning  return  stroke  peak  currents,  J.  Geophys.  Res.,  96,  17135,  1991. 

Papadopoulos,  K.,  G.  Milikh,  A.  Gurevich,  A.  Drobot,  and  R.  Shanny,  Ionization  rates  for 
atmospheric  and  ionospheric  breakdown,  J.  Geophys.  Res.,  98,  17593,  1993. 

Pasko,  V.  P.,  and  U.  S.  Inan.  Recovery  signatures  of  lightning-associated  VLF  perturbations 
as  a  measure  of  the  lower  ionosphere,  J.  Geophys.  Res.,  99,  17,523,  1994. 

Pasko,  V.  P.,  U.  S.  Inan,  Y.  N.  Taranenko,  and  T.  F.  Bell,  Heating,  ionization  and  upward 
discharges  in  the  mesosphere  due  to  intense  quasi-electrostatic  thundercloud  fields, 
Geophys.  Res.  Lett.,  22,  365,  1995. 


124 


BIBLIOGRAPHY 


Pasko,  V.  P,  U.  S.  Inan,  and  T.  F.  Bell,  Sprites  as  luminous  columns  of  ionization  produced 
by  quasi-electrostatic  thundercloud  fields,  Geophys.  Res.  Lett.,  23,  649,  1996a. 

Pasko,  V.  P,  U.  S.  Inan,  andT.  F.  Bell,  Blue  jets  produced  by  quasi-electrostatic  pre-discharge 
thundercloud  fields,  Geophys.  Res.  Lett.,  23,  301,1996b. 

Pasko,  V.  P,  U.  S.  Inan,  T.  F.  Bell,  and  Y.  N.  Taranenko,  Sprites  produced  by  quasi¬ 
electrostatic  heating  and  ionization  in  the  lower  ionosphere,  J.  Geophys.  Res.,  in 
review,  1996c 

Pasko,  V.  P,  U.  S.  Inan,  and  T.  F.  Bell  Ionospheric  Effects  due  to  Electrostatic  Thundercloud 
Fields,  Geophys.  Res.  Lett.,  in  review,  1996d 

Petrov,  N.  I.,  V.  R.  Avanskii,  and  N.  V.  Bombenkova,  Measurements  of  the  electric  field  in 
the  streamer  zone  and  in  the  sheath  of  the  channel  of  a  leader  discharge,  Tech.  Phys., 
39,  546,  1994. 

Potter,  D.,  Computational  physics,  John  Wiley  &  Sons,  London,  1973. 

Picard,  R.  H.,  J.  R.  Winick,  R.  D.  Sharma,  A.  S.  Zachor,  P.  J.  Espy,  and  C.  R.  Harris, 
Interpretation  of  infrared  measurements  of  the  high-latitude  thermosphere  from  a 
rocket-borne  interferometer,  A dv.  Space  Res.,  7,  (10)23,  1987. 

Rairden,  R.  L.,  and  S.  B.  Mende,  Time  resolved  sprite  imagery,  Geophys.  Res.  Lett.,  22, 
3465,  1995. 

Reid,  G.  C.,  The  middle  atmosphere,  NASA  CP-2090,  National  Aeronautics  and  Space 
Administration,  27-42,  1979. 

Roussel-Dupre,  R.  A.,  A.  V.  Gurevich,  T.  Tunnell,  and  G.  M.  Milikh,  Kinetic  theory  of 
runaway  air  breakdown,  Phys.  Rev.  E,  49,  2257,  1994. 

Roussel-Dupre,  R.  A.,  and  A.  V.  Gurevich,  On  runaway  breakdown  and  upward  propagating 
discharges,  /  Geophys.  Res.,  101,  2297,  1996. 

Rowland,  H.  L.,  R.  F.  Femsler,  J.  D.  Huba,  and  P.  A.  Bernhardt,  Lightning  driven  EMP  in 
the  upper  atmosphere,  Geophys.  Res.  Lett.,  22,  361,  1995. 

Rust,  W.  D.,  T.  C.  Marshall,  M.  Stolzenburg,  M.  G.  Boteman,  T.  R.  Shepherd,  B.  F.  Smull, 
Observations  of  electrical  structure  of  mudlatitude  mesoscale  convective  systems: 
Review  and  update,  in  Proceedings  of  10th  International  Conference  on  Atmospheric 
Electricity',  p.  108,  Osaka,  Japan,  June  10-14,  1996. 


BIBLIOGRAPHY 


125 


Sentman,  D.  D.,  and  E.  M.  Wescott,  Observations  of  upper  atmospheric  optical  flashes 
recorded  from  an  aircraft,  Geophys.  Res.  Lett.,  20,  2857,  1993. 

Sentman,  D.  D.,  and  E.  M.  Wescott,  Red  sprites  and  blue  jets.  Geophysical  Institute  Video 
Production,  University  of  Alaska,  Fairbanks,  9  July,  1994. 

Sentman,  D.  D.,  and  E.  M.  Wescott,  Red  sprites  and  blue  jets:  Thunderstorm-excited  optical 
emissions  in  the  stratosphere,  mesosphere,  and  ionosphere,  Phys.  Plasmas,  2,  2514, 
1995. 

Sentman,  D.  D.,  E.  M.  Wescott,  D.  L.  Osborne,  D.  L.  Hampton,  and  M.  J.  Heavner,  Pre¬ 
liminary  results  from  the  Sprites94  campaign:  Red  Sprites,  Geophys.  Res.  Lett.,  22, 
1205,  1995. 

Sipler,  D.  P,  and  M.  A.  Biondi,  Measurements  of  O('D)  quenching  rates  in  the  F  region,  J. 
Geophys.  Res.,  17,  6202,  1972. 

Stergis,  C.G.,  G.C.  Rein,  and  T.  Kangas,  Electric  field  measurements  in  the  stratosphere  J. 
Atmos.  Terr.  Phys.,  11,  77,  1957. 

Swarztrauber,  P.,  and  R.  Sweet,  Efficient  FORTRAN  subprograms  for  the  solution  of  elliptic 
partial  differential  equations,  NCAR  Technical  Note,  NCAR/TN-109+IA,  July,  1975. 

Taranenko,  Y.  N.,  U.  S.  Inan,  and  T.  F.  Bell,  Interaction  with  the  lower  ionosphere  of 
electromagnetic  pulses  from  lightning:  heating,  attachment,  and  ionization,  Geophys. 
Res.  Lett.,  20,  1539,1993a. 

Taranenko,  Y.  N.,  U.  S.  Inan,  and  T.  F.  Bell,  The  interaction  with  the  lower  ionosphere 
of  electromagnetic  pulses  form  lightning:  excitation  of  optical  emissions,  Geophys. 
Res.  Lett.,  20,  2675,  1993b. 

Taranenko,  Y.  N.,  Interaction  with  the  lower  ionosphere  of  electromagnetic  pulses  from 
lightning:  heating,  attachment,  ionization,  and  optical  emisssions,  Ph.D.  dissertation. 
Stanford.  Univ.,  Calif.,  September,  1993. 

Taranenko,  Y.  N.,  and  R.  A.  Roussel-Dupre,  High  altitude  discharges  and  gamma-ray  flashes: 
a  manisfestation  of  runaway  air  breakdown,  Geophys.  Res.  Lett.,  23,  571 ,  1996. 

Taylor,  L.  L.,  Mesospheric  heating  due  to  intense  tropospheric  convection,  NASA  Contractor 
Report  3132,  Washington,  D.  C.,  1979. 


126 


BIBLIOGRAPHY 


Tzur,  I.  and  R.  G.  Roble,  The  interaction  of  a  polar  thunderstorm  with  its  global  electrical 
environment,  J.  Geophys.  Res.,  90,  5989,  1985. 

Valence-Jones,  A.,  Aurora,  D.  Reidel  publishing  Co.,  Dordrecht,  1974. 

Vaughan,  O.  H.,  Jr.,  and  B.  Vonnegut,  Recent  observations  of  lightning  discharges  from  the 
top  of  a  thunderstorm  into  the  clear  air  above,  J.  Geophys.  Res.,  94,  13179,  1989. 

Vaughan,  O.  H.,  Jr.,  R.  J.  Blakeslee,  W.  L.  Boeck,  B.  Vonnegut,  M.  Brook,  and  J.  McKune 
Jr.,  A  cloud-to-space  lightning  as  recorded  by  the  space  shuttle  payload-bay  TV 
cameras,  Mon.  Weather  Rev.,  120,  1459,  1992. 

Volland  H.,  Atmospheric  Electrodynamics,  Springer-Verlag,  New  York,  1984. 

Vitello,  R  A.,  B.  M.  Penetrante,  and  J.  N.  Bardsley,  Simulation  of  negative-streamer  dynam¬ 
ics  in  nitrogen,  Phys.  Rev.  E,  49,  5574,  1994. 

Wen,  C.,  and  J.  M.  Wetzer,  Electron  avalanches  influenced  by  detachment  and  conversion 
processes,  IEEE  Trans.  Electr.  Instil.,  23,  999,  1988. 

Wescott,  E.  M,  D.  Sentman,  D.  Osborne,  D.  Hampton,  and  M.  Heavner,  Preliminary  results 
from  the  Sprites94  aircraft  campaign:  2.  Blue  jets,  Geophys.  Res.  Lett.,  22,  1209, 
1995a. 

Wescott,  E.M.,  D.  D.  Sentman,  M.  J.  Heavner,  and  D.  L.  Hampton,  Blue  starters,  discharges 
above  an  intense  thunderstorm  over  Arkansas,  July  1,  1994,  EOS  Trans.  AGU,  1995 
Fall  Meeting,  76,  FI 04,  1995b. 

Wescott,  E.  M.,  D.  D.  Sentman,  M.  J.  Heavner,  D.  L.  Hampton,  D.  L.  Osborne,  and  O.  H. 
Vaughan  Jr.,  Blue  Starters:  Brief  upward  discharges  from  intense  arkansas  thunder¬ 
storm,  Gephys.  Res.  Lett.,  23,  in  press,  1996. 

Wilson,  C.  T.  R.,  The  electric  field  of  a  thundercloud  and  some  of  its  effects,  Proc.  Phys. 
Soc.  London,  37,  32D,  1925. 

Wilson,  C.  T.  R.,  A  theory  of  thundercloud  electricity,  Proc.  R.  Soc.  London  Ser  A,  236, 
297,  1956. 

Winckler,  J.  R.,  Further  observations  of  cloud-ionosphere  electrical  discharges  above  thun¬ 
derstorms,  J.  Geophys.  Res.,  100,  14335,  1995. 

Winckler,  J.  R.,  W.  A.  Lyons,  T.  Nelson,  and  R.  J.  Nemzek,  New  high  resolution  ground- 
based  studies  of  cloud-ionosphere  discharges  over  thunderstorms  (Cl  or  Sprites),  J. 
Geophys.  Res.,  101,  6997,  1996. 


BIBLIOGRAPHY 


127 


Uman,  M.  A.,  Lightning,  Dover  Publications,  1984. 

Uman.  M.A.,  The  lightning  discharge,  Academic  Press,  Orlando,  1987. 

Uman,  M.A.,  The  earth  and  its  atmosphere  as  a  leaky  spherical  capacitor,  Am.  J.  Phys.,  42, 
1033,  1974. 


